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ENERGY  TRANSFER  AND  DISSOCIATION  IN  HYPERTHERMAL 
ATOM-DIATOM  COLLISIONS 

By 

Lyntis  H.  Beard,  Jr. 
March  1979 

Chairman:    David  A.  Micha 
Major  Department:  Chemistry 

Energy  transfer  and  dissociation  processes  in  hyperthermal  atom- 
diatom  collisions  are  investigated  using  approximations  based  on  the 
multiple  collision  expansion  of  the  Faddeev  equations.    Within  this  many- 
body  approach,  the  collision  process  is  treated  as  a  sequence  of  atom- 
atom  encounters.    This  allows  one  to  obtain  information  on  a  three-atom 
system  using  only  two-atom  transition  operator  matrix  elements. 

A  new  numerical  procedure,  the  "variable-phase  and  -amplitude" 
method,  is  developed  to  compute  the  required  off-energy-shel 1  two-body 
t-matrix  elements  for  an  arbitrary  radial  diatomic  potential.    An  alter- 
nate approach  based  on  the  comparison  potential  method  is  also  investi- 
gated, leading  to  an  improved  computational  algorithm.    Whereas  these 
methods  are  based  on  a  partial  wave  decomposition  of  the  two-body  t-opera- 
tor,  the  need  for  non-partial  wave  techniques  is  pointed  out  at  high 
energies  and  large  momentum  transfers.    In  this  connection  the  Bateman 
method  is  investigated. 

In  the  present  study,  the  single  coll ision  approximation  is  considered 
for  both  collinear  and  three-dimensional  motions.    In  the  collinear  studies, 
it  is  shown  that  multiple  collision  effects  are  important  in  the  case  of 
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inelastic  scattering  near  threshold,  and  that  for  the  case  of  dissocia- 
tive scattering  they  become  crucial.    A  study  is  also  made  of  the  valid- 
ity of  the  single  collision  peaking  approximation,  and  it  is  shown  that 
this  approximation  breaks  down  when  soft  two-atom  potentials  are  involved. 
In  order  to  compare  to  experiment,  three-dimensional  results  are  presented 
for  the  scattering  of  Li    with  N?  and  CO  at  hyperthermal  energies. 
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INTRODUCTION 

With  the  advent  of  new  experimental  techniques,  the  modern-day 
chemist  is  now  more  than  ever  able  to  discern  the  specificity  and  selec- 
tivity shown  by  the  reactants  and  products  in  a  chemical  system.  The 
goal  of  a  chemist  is  to  properly  interpret  and  understand  in  a  general 
way  the  events  that  occur  in  a  chemical  reaction.    The  task  at  nand, 
however,  is  quite  formidable  because  of  the  complexity  of  the  particles 
involved.    This  complexity  implies  that  the  investigation  of  the  simples 
chemical  systems  is  necessary  to  lay  the  foundation  of  a  deeper  and  more 
conceptual  understanding  of  chemical  phenomena. 

In  this  dissertation  we  will  be  concerned  with  the  investigation 
of  energy  transfer  processes  in  atom-diatom  collisions  at  hyperthermal 
energies.    In  particular,  we  shall  be  concerned  with  processes  of  the 
type 

A  +  BC  ■+  A  +  BC 
+  B  +  AC 
C  +  AB 
+  A  +  B  +  C  , 

where  A,  B  and  C  represent  the  atoms  in  the  atom-diatom  system.  Theoret 
ically,  atom-diatom  systems  are  of  interest  because  they  exhibit  the  en- 
tire spectrum  of  possible  modes  of  energy  transfer,  i.e.,  they  allow  for 
vibrational  and  rotational  energy  transfer,  as  well  as  electronic  and 
translational  energy  transfer.    Therefore,  we  see  that  the  elucidation 
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of  the  dynamics  of  atom-diatom  collisions  would  provide  an  impetus  to 
the  eventual  understanding  of  more  complicated  systems.    Our  goal  in 
this  chapter  will  be  to  clarify  and  discuss  the  treatment  of  electronic 
structure  in  atom-diatom  scattering  processes.    The  theoretical  aspects 
of  the  scattering  problem  will  be  the  subject  of  the  remaining  chapters 
in  this  dissertation. 

1.    Theoretical  Treatment  of  Nuclear  and  Electronic  Motions 
A  formal  treatment  of  an  atom-diatom  system  requires  that  one  con- 
sider both  the  nuclear  and  electronic  degrees  of  freedom  of  the  system. 
In  other  words,  one  would  have  to  solve  a  many-particle  problem.  From 
a  practical  standpoint,  such  a  problem  would  not  be  tractable  unless 
approximations  were  introduced.    In  this  section  therefore,  we  will 
sketch  various  approaches  that  can  be  taken  to  treat  the  atom-diatom 
system  as  an  effective  three-body  problem. 

The  first  approach  we  shall  consider  is  that  of  introducing  approx- 
imations which  allow  for  the  separation  of  nuclear  and  electronic  motions. 
We  begin  then  by  considering  the  Schrb'dinger  equation 

(H-E)y  =  0  (1.1) 

where  E  is  the  energy  of  the  system,  ¥  is  the  wave  function  and  the  Hamil- 
tonian  H  can  be  written  as  (in  atomic  units) 

H=-^  2^Va"i.V?  +  Vcoul  +HS0  •  t1-1) 
a= 1     a  1=1 

in  the  center-of-mass  system  (Fr62).  The  Greek  indices  in  Eq.  (1.2)  refer 
to  the  nuclei  and  the  Latin  indices  refer  to  the  N  electrons  in  the  atom- 
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diatom  system.    We  note  that  the  Coulomb  interactions  of  the  various 
particles  involved  are  contained  in  Vcou-j  and  the  spin-orbit  interactions 
are  contained  in  HSQ.    The  electronic  and  nuclear  mass  polarization 
terms  have  been  neglected.    The  reduced  mass       is  defined  by  the  rela- 


tion 


_L  =  _L  _  4  (1.3) 

u        m        M  u  ' 

Ma  a 

where  m    is  the  mass  of  nucleus  a  and 
a 

3 

M  =    Z   m0  (1.4) 
6=1  B 

is  the  total  mass  of  the  nuclei.    In  order  to  separate  nuclear  and  elec- 
tronic degrees  of  freedom,  the  wave  function  ¥  is  generally  expanded  in 
terms  of  a  set  of  electronic  functions  {<f>-|(ra;  r-j)>  which  depend  on  the 
electronic  coordinates  {jr;}  and  only  parametrically  on  the  nuclear  coor- 
dinates {r  }.    Thus  one  has 
—a 

*(v  Hi> =  { <Mv  ziWzJ  '  (1-5) 

where  Xid^)  can  be  considered  as  the  nuclear  wave  function  describing 
the  motion  of  the  nuclei  on  the  potential  energy  surface  which  charac- 
terizes the  Tth  electronic  state  [(Fr62),  (Hi67)].    Equations  (1.5)  and 
(1.1)  lead  to  the  following  set  of  coupled  equations: 

*VTiV  uirE]xi  =  Y^iv  +TiV  +  uir]*v  ■  (1'6) 

where  we  have  introduced  the  definitions 

'.--ij^i'  (1-7) 
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N 

H  ,  =  -  I  I  7?  +  V„rtl  ,  ,  (1.8) 
el        2         i  coul 


=1 


U„.  •  U{«]>  +  uft?>  .  (1.9) 
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x'  ,  =  Z  d{?!  '       ,  (1.12) 


and 


IT    ^ZP-^T  '3 


TiV  =  z  srD{r  (1-13) 

1 1       3=1  B 


djf)  -  <4»1|^1,>  (1.14) 


D1V  =  ^ll^T3,  '     '  (1'15) 


These  equations  are  well  known  in  the  literature  [(Fr62),  (Hi 67),  (Tu76a)], 
and  in  general  form  the  basis  of  most  treatments  of  atomic  and  molecular 
scattering  problems.    Our  notation  and  discussion  follow  that  of  Tully 
(Tu76a)].The  diagonal  matrix  elements  U^r^)  are  the  effective  potential 
energy  surfaces  for  the  atom-diatom  system  and  the  nondiagonal  matrix 
elements  x-| -j , ,       ,  and       ,  are  the  coupling  elements  responsible  for 
the  promotion  of  electronic  transitions.    The  diagonal  term        is  a  non- 
adiabatic  correction  to  U^.    The  diagonal  term  x-j.  can  be  shown  to  be 
exactly  zero  if  the  electronic  functions  {<J>, }  are  taken  to  be  real  (Hi 67 ) . 


The  neglect,  of  all  the  momentum  dependent  coupling  terms  T-j-j,  and  , 
(all  1  and  T)  is  known  as  the  Born-Oppenheimer  approximation  [(Hi67), 
(Tu76a)]and  leads  to  a  new  set  of  coupled  equations 

K  +  Uii  "  E]*l  =  "r2=1  U1TXT  (1'16) 

in  which  electronic  transitions  may  occur  only  through  the  nondiagonal 
potential  terms  U,-. , . 

The  effectiveness  of  Eq.  (1.6)  in  properly  describing  collision 
processes  involving  electronic  transitions  depends  crucially  on  the  choice 
of  the  functions  For  example,  if  {^}  is  chosen  to  belong  to  the 

adiabatic  representation,  i.e.,  that  representation  which  diagonal izes 
the  Hamilton1! an 

H'  =  Hel  +  HSQ  ,  (1.17) 
one  obtains  the  set  of  uncoupled  equations 

[xa  +  V1  -  E]Xl  *  0  ,  (1.18) 


where 


V,(r  )  =  U,,(r  )  .  (1.19) 
1  -a         1 l x-a' 


In  this  representation,  transitions  between  the  various  electronic  states 
are  not  allowed,  unless  one  considers  the  momentum  dependent  coupling 
terms.    The  term  adiabatic  representation  is  also  used  in  the  literature 
for  that  representation  in  which  the  electronic  Hamiltonian  H  -j  is  diag- 
onal [ (Hi67) ,  (Tu76a)].  There  are,  however,  an  infinite  number  of  non- 
adiabatic  representations.    One  set  of  representations  encountered  in 
the  literature  are  the  so  called  diabatic    representations  [(0m71), 
(Sm69)].    These  are  representations  in  which  the  electronic  functions 


approach  well  defined  atoric  states  asymptotically  (0m71).    There  is  no 
unique  recipe  for  constructing  such  representations.    In  general  they 
are  chosen  such  that  the  momentum-dependent  coupling  terms  t^,  are 
small  (Sm69).    One  should  note  however,  as  Tully  (Tu76a)  points  out.  that 
diabatic  representations  are  useful  in  atom-atom  scattering  because  in 
those  situations  avoided  crossings  in  the  adiabatic  curves  are  localized 
to  a  point.    In  other  situations  where  there  are  more  nuclear  degrees 
of  freedom,  as  in  atom-diatom  scattering,  one  is  dealing  with  avoided 
surface  crossings  at  which  large  variations  in  the  nonadiabatic  couplings 
can  occur.    This  indicates  that  one  might  not  be  able  to  find  a  diabatic 
representation  which  is  suitable  in  a  global  sense. 

2.    Adiabatic  Collision  Processes 
Because  of  the  complexity  of  atomic  and  molecular  collisions,  most 
work  in  the  literature  has  focused  on  the  theory  of  adiabatic  collision 
processes.    The  theoretical  investigation  of  such  processes  is,  however, 
by  no  means  trivial.    For  example,  in  the  case  of  inelastic  atom-diatom 
scattering,  Eq.  (1.18)  would  generally  involve  solving  a  large  set  of 
coupled  differential  equations.    If  a  large  number  of  open  channels  for 
vibrational  and  rotational  excitation  exist,  an  accurate  and  economical 
solution  to  Eq.  (1.18)  may  be  difficult  to  obtain.    In  order  to  solve 
the  problem  for  such  collision  processes,  further  approximations  have 
been  sought  which  reduce  the  dimensionality  of  such  coupled  differential 
equations.  A  review  of  these  techniques  is  given  by  Micha  (Mi74),  Lester 
(Le76)  and  Rabitz  (Ra76).    In  the  case  of  reactive  scattering  the  situa- 
tion is  even  more  complicated.    For  an  atom-diatom  system,  there  are 
four  arrangement  channels  [A  +  BC,  B  +  AC,  C  +  AB  and  A  +  B  +  C]  and  as 


3  consequence  *ree  motion  in  the  asymptotic  region  is  governed  by  differ- 
ent channel  Hamiltonians.    This  aspect  of  reactive  scattering  introduces 
compli cations,  because  the  boundary  conditions  become  difficult  to  apply 
and  in  addition  no  unique  and  simple  set  of  coordinates  is  suited  for 
propagating  the  solution  in  all  regions  of  configuration  space.    We  note 
that  in  each  rearrangement  channel,  a  natural  choice  of  coordinates  would 
be  that  of  the  Jacobi  coordinates  t.  and  r . ,  where  t.  describes  the  rela- 
tive motion  of  the  "free"  atom  with  respect  to  the  center  of  mass  of  the 
diatom  and  r.  the  relative  motion  of  the  diatom.    The  label  "i"  is  a 
channel  index  (see  Appendix  IV).    In  order  to  overcome  the  difficulty 
with  the  choice  of  coordinates,  Marcus  (Ma66)  introduced  the  concept  of 
natural  collision  coordinates,  e.g.,  a  set  of  coordinates  which  goes 
smoothly  from  one  set  of  Jacobi  coordinates  to  the  other.  However, 
these  coordinates  complicate  the  structure  of  the  kinetic  energy  opera- 
tor and  may  in  some  instances  become  multivalued  according  to  how  they 
are  constructed  [(Ma66),  (W177)].    A  clarification  of  many  of  the  problems 
of  reactive  scattering  may  be  obtained  by  considering  the  structure  and 
many-body  aspects  of  the  wave  function  [(Mi72b),  (Re77)].    A  general 
review  of  reactive  scattering  may  be  found  in  the  articles  by  George  and 
Ross  (Ge71),  Kouri  (Ko73)  and  Micha  (Mi75a). 

3.    Nonadiabatic  Collision  Processes 
Nonadiabatic  collision  processes  have  not  been  studied  in  as  much 
detail  as  adiabatic  processes  for  two  general  reason.    One  reason  is  the 
added  complexity  of  dealing  with  a  coupled  set  of  equations  [Equations 
(1.6)  and  (1.16)];  the  second  reason  is  the  lack  of  information  on  the 
potential  energy  surfaces  and  on  the  various  couplings  involved.    It  has 
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been  found,  however,  that  for  some  problems  it  is  possible  to  overcome 
the  first  of  these  obstacles  by  treating  some  or  all  of  the  nuclear  de- 
grees of  freedom  either  classically  or  semiclassically.    As  was  pointed 
out  earlier,  in  some  instances  nonadiabatic  effects  may  be  localized  to 
certain  regions  of  the  potential  energy  surfaces.    This,  for  example, 
would  allow  one  to  use  classical  trajectory  methods  in  regions  where 
motion  along  the  adiabatic  potential  surface  is  allowed.    In  such  an 
approach  the  nonadiabatic  regions  serve  only  to  determine  branch  points 
at  which  a  jump  to  another  surface  can  occur.    This  approach,  which  is 
known  as  the  trajectory  surface  hopping  model  (TSH),  was  proposed  by 
Tully  and  Preston  (Tu71)  to  study  the  systems  [H+,D2]  (Tu71),  [H,D*] 
(Pr73)  and  [Ar+,H2]  (Ch74).    In  order  to  determine  whether  a  trajectory 
remains  in  the  same  potential  energy  surface  or  branches  to  another  sur- 
face, the  probability  for  hopping  must  be  computed.    A  discussion  of 
how  this  probability  is  obtained  is  given  in  the  above  references.  We 
shall  only  mention  one  method;  that  of  the  well  known  Landau-Zener- 
Stuckelbert  (LZS)  approximation  [see  for  example  Bransden  (Br70)]. 
Strickly  this  approximation  is  valid  only  for  the  one-dimensional  case. 
Here,  however,  one  applies  the  LZS  result  to  a  cut  perpendicular  to  the 
nonadiabatic  seam  (Tu76a)  in  the  potential  energy  surfaces.    Since  trajec- 
tories may  encounter  various  branch  points,  we  see  that  in  many  physical 
applications  the  TSH  method  could  become  computationally  unmanageable. 

Another  approach,  which  is  useful  in  treating  nonadiabatic  processes 
such  as 

A  +  BCCn^j)  -  A+  +  BC'Cn^j^  (1.20) 
+  A    +  BC  (nj,jj)  (1.21) 
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which  involve  charge  transfer  and  vibrational-rotational  excitation,  is 
that  given  by  the  multiple-curve  crossing  model.    In  order  to  define  this 
model,  let  us  introduce  the  Jacobi  coordinates  ^  and      which  are 
appropriate  in  describing  the  arrangement  channel  specified  by  Eq.  (1.21). 
Furthermore,  let  us  expand  the  nuclear  wave  function  x-j^,^)  in  terms 
of  a  set  of  vibration-rotation  functions  (u^^i^)}  as  follows, 

X  (*  ?  )=E       (ftjfye^Cfy.  (1.22) 
v 

where  v  labels  the  vibration-rotation  states  and  ^(Rj)  acts  as  a  nu- 
clear wave  function  for  the  relative  motion  of  atom  A  with  respect  to 
the  center  of  mass  of  the  diatom  BC.    Note  that  the  {ulv}  depend  only 
parametrically  on  %.    If  we  now  use  Eq.  (1.6)  expressed  in  Jacobi  coor- 
dinates  and  Eq.  (1.22),  we  obtain  a  set  of  coupled  differential  equations 
given  by 

[*E  +  <ulvlTll'ulv>  +  Wlvlv-  E^lv(Rl}  = 

-  ry<UlvlTll'  +  T11'IU1V>  +  WlvlV]bW(1*l>'  (1-23) 


v' 


where 


E-       1    V2       ,  (1.24) 


T  = 


m.(mR+mr) 
1  WmC 


Tn-  = "  iq^ilfyV  '  %  •  (1'26) 
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T 


E"  1    ,±  ,„2 


jr^iK  M  '  (1-27) 


ir     2Ml    ,  ^ 


wiviv(]V  5<hv,l6irTSS'^ii**<*i,Tl+H,,*i,>,uiv>  •  (1*28) 


ml     ~  (mB+mc)  ' 


(1.30) 


Tir--^<*i^p1i*i,>^pr  {1'31) 


I"         1  i„2 

T 


IT    "  2mx      K^T*    '  (1'32) 


and  the  quantities  m^,  mB  and  mc  are  the  masses  of  the  atoms  composing 
the  system.    Equation  (1.23)  is  the  basis  for  treating  an  inelastic  atom- 
diatom  scattering  problem.    The  diagonal  matrix  elements  W]v-jv(^j_) 
correspond  to  effective  potential  energy  curves,  each  of  which  is  asso- 
ciated to  a  particular  electronic-vibrational -rotational  state  (Tu76). 
These  potential  energy  curves  then  serve  as  a  basis  for  treating  charge 
transfer  and  inelastic  scattering  within  the  multiple-curve  crossing 
model  [(Ba69b),  (Fi71),  (Gi75),  (Tu76a)].    It  is  obvious,  however,  that 
due  to  the  enormous  number  of  these  curves,  various  other  approximations 
must  be  introduced.    Generally  one  averages  over  the  rotational  states 
and  assumes  that  nonadiabatic  effects  are  confined  to  the  avoided-cross- 
ing  regions  where  the  LZS  approximation  can  be  used.    A  discussion  of 
the  approximations  and  implementation  of  this  model  can  be  found  in  the 


work  of  Fisher  et  al .  [(Fi71),  (Fi72),  (Ba69b)],  Gislason  and  Sachs 
(Gi75),  and  in  the  review  article  by  fully  (Tu76a).    We  note  that  Eq. 
(1.23)  is  not  suitable  for  treating  rearrangement  collision  processes. 
However,  using  the  many-body  formalism  of  the  Faddeev  equations  Yuan  and 
Micha  (Yu76)  have  implemented  multiple  curve  crossing  information  to 
study  rearrangement  scattering  of  the  atom-diatom  system  [K,IBr]. 

4.    Many-Body  Theory  and  the  Treatment  of  Nuclear  and  Electronic  Motions 

All  the  methods  discussed  in  the  previous  sections  involved  the 
treatment  of  electronic  and  nuclear  motions  essentially  as  two  separate 
problems.    This  approach  has  resulted  in  the  introduction  of  the  concept 
of  a  potential  energy  surface,  i.e.,  that  quantity  which  characterizes 
the  electronic  state  of  the  system  and  serves  as  the  effective  potential 
governing  nuclear  motion.    An  alternative  to  this  picture  would  be  the 
simultaneous  treatment  of  nuclear  and  electronic  motions.    We  will 
briefly  sketch  how  this  can  be  accomplished  by  the  use  of  many-body 
theories.    This  will  be  done  for  three  reasons.    First,  it  will  serve 
to  contrast  the  more  traditional  methods  of  scattering  theory  based  on 
Eq.  (1.6).    Secondly,  it  will  show  how  one  can  use  known  two-body  data 
as  input  to  the  solution  of  an  inherent  many-body  problem.    Lastly,  it 
will  serve  as  a  reference  with  which  to  compare  the  approximate  treat- 
ment of  atom-diatom  collisions  given  in  later  chapters  of  the  present 
work. 

The  many-body  formalisms  that  will  be  discussed  are  based  on  the 
treatment  of  the  three-body  problem  pioneered  by  Faddeev  (Fa61).  Earlier 
it  was  pointed  out  that  for  a  three-particle  system,  there  are  four 
arrangement  channels.    This  implies  that  free  motion  in  each  asymptotic 


region  is  governed  by  a  different  channel  Hamilton! an.    In  order  to 
clarify  this,  consider  the  Hamiltonian  for  three  interacting  atoms  A,  B 
and  C  (in  laboratory  frame), 


H=ZK   +  Z  HT  +  I  V 

OC  T  I  I 

a  I  J 


JK 


(1.33) 


where  K    is  the  kinetic  energy  of  nucleus  a  =  a,  b,  or  c  (associated  with 

a 

A,  B  and  C  respectively);  Hj  is  the  electronic  Hamiltonian  of  atom  I, 
and  V,,,  is  the  Coulomb  interaction  among  the  charges  of  atoms  J  and  K 
(Mi 77b).    The  channel  Hamiltonian  for  a  given  rearrangement  channel 
i[l  for  (A+BC),  2  for  (B+AC),  3  for  (C+AB)]  is  defined  by  the  relation 

H.  =  lim   H    ,  (1.34) 
1  R.-*» 

and  as  can  be  seen  using  Eq.  (1.33) 


In  the  breakup  arrangement  channel  one  can  analogously  define  the  free 
particle  Hamiltonian 


H.  =  H    +  E  HT  +  V 
l  ex      T  I 


JK  ' 


(1.35) 


(1.36) 


which  can  be  used  in  Eq.  (1.35)  to  obtain 


and 


(1.37) 
(1.38) 


where  V.  is  equal  to  Vlt/.    The  Faddeev  formalism,  then,  is  merely  a  way 

1  J  Is. 

of  decomposing  the  total  wave  function  of  the  system  into  channel  com- 
ponents which  can  be  propagated  by  their  respective  channel  Hamiltonians. 
The  result  is  given  by  the  relation  (Fa61) 

3 

(E  -  H.)y(j)  =  V.    I  (1  -  6., )¥(k),  (j  =  1,  2,  3)  (1.39) 
J  J  k=1  JK 

where  the  total  wave  function  ¥  has  been  decomposed  into  the  channel  com- 
ponents ¥    ^ ,  i .e 


4»  =    I  .  (1.40) 

j=l 

A  detailed  derivation  of  the  Faddeev  equations  will  be  given  in 
Chapter  IV  of  the  present  work.    In  the  following  discussion,  our  con- 
cern will  be  the  use  of  Eq.  (1.39)  solely  as  the  basis  for  the  simul- 
taneous treatment  of  nuclear  and  electronic  motions.    Micha  (Mi 77d)  has 
shown  that  this  can  easily  be  accomplished  within  the  framework  of  the 
method  of  diatomics-in-molecules  (DIM)  (E163).    We  begin  by  considering 
the  total  wave  function  for  the  system,  keeping  in  mind  that  for  struc- 
tured particles  it  must  have  the  form 

¥  =  Cj\X  ,  (1.41) 

where  C  is  a  normalization  constant,  J{  is  the  N  electron  antisymmetrizer 
and  X  is  an  unsymmetrized  total  wave  function  for  both  nuclear  and  elec- 
tronic motions.    Since  H  and  J\  commute,  one  can  cast  the  Schrb'dinger 
equation  into  the  form 


[EjA  -  |  (Hj\         H)]X  =  0  ,  (1.42) 


or  in  the  Faddeev  form 

3  / 

(M  -  fr*yJ  Vx(j)  -  hyy*  h  k^(1-6jk)x     •  (1-43) 

where 

[Hj.jl  ]+  s  +  ^  H.    .  (1.44) 

Equation  (1.39)  may  now  be  replaced  by  Eq.  (1.43),  which  is  more  suitable 
for  introducing  DIM  approximations.    A  symmetrized  form  of  the  Hamil- 
tonian  such  as 

i  HjA    -  jljl  »H]+  (1.45) 

has  been  introduced  in  order  to  avoid  nonhermitian  matrices  (Tu76a)  which 
occur  whenever  the  x'^'  are  expanded  in  terms  of  a  set  of  channel  elec- 
tronic states  |£    >  and  used  in  Eq.  (1.43).    The  channel  electronic 
states  |£^>  correspond  to  a  collection  of  m'^'  primitive  channel  wave 
functions  constructed  in  a  step  by  step  procedure  from  a  set  of  atomic 
wave  functions  g.j,  cBm  and  ?Cm,  as  is  usually  done  in  the  DIM  method 
[ ( El 63) ,  (Mi 77d) ] .    When  these  channel  functions  are  constructed,  how- 
ever, the  arrangement  structure  of  the  channel  must  be  considered. 
Therefore,  one  begins  by  constructing  a  collection  of  Mq  functions 

2 

{ l5A1CBn5Cm>}.    Each  of  the  atomic  functions  Cj-j  are  eigenfunctions  Sj 
and  Sj^,  and  are  composed  of  nonsymmetrized  combinations  of  products  of 
atomic  orbitals  centered  on  the  nucleus  of  atom  I.    This  set  of  func- 
tions we  designate  |£q>.    A  set  of  M^  un symmetrized  or  primitive 

channel  wave  functions  |£    >  can  be  obtained  from  |£^>  by  first  com- 

2 

bining  the  products  £jm?Kn  to  form  eigenfunctions.  of      and  S^,  and 


then  combining  the  products  of  these  and  £•,,  to  form  eigenf unctions  of 

0 

S    and  Sz  (Mi77d).    This  leads  to  the  relations 

<ktee =  k  {lA7) 

and 

<^(j)|^(j)>  =  jU)  (1>48) 

where  J_q  and         are  unit  matrices  and         defines  a  channel  coupling 
matrix  between  the  breakup  channel  (channel  0)  and  the  arrangement 
channel  i.    The  unsymmetrized  or  primitive  channel  wave  functions  x'J^ 
can  now  be  expressed  as 

x(j)  s  |r(j)>x(j)  ,  (1.49) 

where         denotes  a         x  1  column  of  scattering  functions.    In  ob- 
taining Eq.  (1.49),  one  has  assumed  the  approximate  completeness  rela- 
tion 

I^W^I  =  I    ,  (1.50) 

u 

where  1.  is  the  unit  operator  in  the  space  of  many  electron  states  in 
channel  j  (Mi77d).    Using  equations  (1.43)  and  (1.49)  the  end  result  of 
the    DIM-Faddeev  equations  is  given  by 

{A(j)E  -  E  K    -  ^-[A(j),  E  Ej  +  V.]+}  x(j)  = 
a  I  J 

|[A(j),  Vj]+  E  (1  -  5jk)  C(jk)  X(k)  .  (1.51) 

where 

A(j)  =  <i(j)|j^  |?(j)>  ,  (1.52) 


h  -  £ik  -  k  -  h  •  (i-54' 

£j  is  a  diagonal  matrix  containing  atomic  energies  and  E^.,  is  a  matrix 
containing  diatomic  energies.    The  energies  Ej  are  associated  with  the 
unbound  atom  in  channel  j  and  the  E^  energies  of  the  corresponding  dia- 
tomic in  the  same  channel.    Equation  (1.51)  is  a  result  of  the  Born- 
Oppenheimer-type  approximation 

<I(j)|jA  Kall(j)>  *  A(j)  Ka  (1.55) 
and  the  DIM  approximations 

<C(j)|j\  Hj|C(j)>  -  A(j)  Ej  (1.56) 


and 


Huli(j)>  "  A(j)  Iik  •  (1'57) 


( i  k ) 

The  quantities  C/J '  in  Eq.  (1.51)  act  as  recoup! ing  coefficients  among 
the  various  channels,  allowing  for  inelastic,  rearrangement  and  dissoci- 
ative processes  to  occur.    Furthermore,  as  long  as  the  proper  atomic 
and  diatomic  energies  are  included,  Eq.  (1.51)  allows  for  electronic  ex- 
citation.   A  detailed  application  of  this  formalism  to  the  [K  +  Hg]  sys- 
tem is  given  in  the  literature  [ (Mi 72b) ,  (Mi77d)].    In  these  references 
it  was  shown  that  one  can  construct  a  suitable  expression  of  the  ground 
state  potential  energy  surface  of  the  system  in  terms  of  two-body  inter- 
actions provided  one  uses  spin-dependent  potentials. 

Atom-diatom  collisions  involving  charge  transfer  cannot  be  treated 
properly  using  Eq.  (1.51)  because  this  equation  assumes  a  unique  parti- 


tion  of  the  electrons.    There  is,  however,  an  alternative  approach  based 
on  the  coupled-channel  many-body  formalism  [(Ko75),  (To74),  (Ra77)], 
which  is  suitable  for  treating  charge  transfer.    In  this  approach  one 
simply  considers  various  partitions  of  the  Hamiltonian  consistent  with 
the  atomic  and  molecular  fragments  involved  in  the  scattering  process 
of  interest.    A  set  of  coupled  equations  may  then  be  constructed  that 
are  analogous  to  the  Faddeev  equations  (Mi 76b). 

5.  Analysis  of  Potential  Energy  Surface  Information 
In  the  present  work,  we  shall  be  concerned  with  ion-molecule  reac- 
tions of  the  type  (Li+,  CO)  and  (Li  +  ,  Np).  Since  CO  and  N2  have  ident- 
ical masses  and  very  similar  molecular  properties,  these  systems  are 
very  interesting  for  comparison  purposes.  Such  a  study  is  also  attrac- 
tive because  a  comparison  between  experiment  and  theory  is  now  possible 
in  light  of  recent  molecular  beam  experiments  [(Bo76),  (Ea78)]. 

Our  investigation  of  these  systems  will  be  based  on  an  analysis  of 
the  various  interaction  potentials  involved  and  upon  the  potential  energy 
surface  data  available  in  the  literature  [(St75),  (St76),  (Th78a)]. 
This  analysis  though  qualitative  in  nature,  will  serve  as  a  guide  for 
the  development  of  an  impulsive  model  which  will  be  specified  in  the 
next  section. 

The  only  available  potential  energy  surface  Information  on  the 
(Li+,  CO)  and  (Li+,  N2)  systems  corresponds  to  a  restricted  number  of 
geometries  for  the  ground  state  of  the  system.    Because  of  the  relatively 
large  number  of  electrons  in  these  systems,  accurate  boundstate  calcula- 
tions are  very  expensive.    This  is  further  complicated  by  the  large 
number  of  points  required  to  adequately  characterize  the  potential 
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energy  surface.    The  available  calculations  [(St75),  (St76),  (Th78a)] 
were  performed  within  the  SCF  approximation  and  corrected  to  some  degree 
for  correlation  effects.    SCF  calculations  alone  cannot  sufficiently 
describe  the  potential  energy  surface  because  such  an  approximation 
fails  to  properly  describe  the  dissociation  of  closed  shell  systems  such 
as  No  and  CO  into  open  shell  atoms  (St76).    SCF  methods  in  general  fail 
to  properly  characterize  the  large  distance  multipole  type  interactions 
between  atomic  and  molecular  fragments.    Staemmler  [ (St75 ) ,  (St76)]  and 
Thomas  et  al.  (Th78a)  corrected  the  SCF  results  by  making  use  of  CI 
(configuration  interaction)  methods. 

Since  the  information  obtained  in  the  above  calculations  is  re- 
stricted to  the  ground  states  of  the  (Li+,  N«)  and  (Li+,  CO)  systems, 
the  possibility  of  nearby  excited  states  must  be  ascertained  by  consider- 
ing the  possible  electronic  states  of  the  various  atom  and  diatom  frag- 
ments of  these  systems  (Ma75).    For  example,  Fig.  (1-1)  considers  some 
of  the  electronic  states  of  the  (Li+,  N«)  system  when  Li+  and  N9  are 
infinitely  separated.    The  resulting  potential  energy  curves  were  con- 
structed using  the  information  given  by  Bond  et  al.  (Bo65).    We  note 
that  the  ground  state  of  N2(1E^  state)  correlates  with  the  two  nitrogen 
atoms  in  the  ^S  state.    Using  the  Wigner-Witmer  rules  (He50),  one  can 
predict  that  the  other  states  that  correlate  with  the  [N(  S)  +  N(^S)] 
asymptotic  state  are  the  and  7Z~  states.    In  Fig.  (1-1)  the 

u         y  u 

possibility  of  charge  transfer  (Li,  N«)  and  of  electronic  excitation  of 

nitrogen  (Li  ,  N0)  is  also  considered.  Asymptotically  we  see  that  the 

ground  state  of  n!(22^)  correlates  with  [(N(4S)  +  N+(3P)].    Of  more 

interest  is  the  excited  electronic  state  which  asymptotically  leads  to 

A '  2  2 

[N(  S)  +  N(  D)],  where  the    D  state  corresponds  to  the  first  excited 


Figure  (1-1)    Potential  energy  curves  for  the  diatoms 
in  the  reactant  and  product  regions  of 
the  (Li+,  N2)  system  [(Bo65).  (Ma75)] . 
The  Li  +  (Li)  atom  and  the  N2(N;>)  molecule 
are  at  infinite  separation. 


state  of  nitrogen.    Mote  that  in  constructing  the  curves  in  Fig.  (1-1), 
+  1  2 

the  energies  of  Li  (  S)  and  Li(  S)  must  be  included.    The  possibility  of 

reactive  channels  is  ruled  out  by  the  closed  shell  nature  of  the 

+  + 
(Li  ,        system.    Furthermore,  the  electronic  excitation  of  Li 

[Li  (  S)  +  N^(  Zg)]  is  not  considered  since  it  would  lie  much  higher  in 

energy  than  those  states  already  given  in  Fig.  (1-1).    It  is  apparent 

then  from  Fig.  (1-1),  that  there  is  a  possibility  of  electronic  excita- 

tion,  since  the    Zu  state  of  Ng  lies  only  6.1  eV  (Bo65)  above  that  of 

the  ground    Eg  state.    Charge  transfer  on  the  other  hand  would  seem 

unl ikely. 

As  the  Li+  is  brought  near  N«,  the  relative  position  of  the  various 

electronic  states  would  change.    If  the  spacing  between  the  ground  and 

excitated  states  grows  closer,  nonadiabatic  effects  will  be  important 

for  the  given  system.    Figure  (1-2)  illustrates  what  happens  as  the  Li  + 

atom  approaches  the  N2  molecule.    This  figure  shows  three  planes. 

Plane  III  corresponds  to  the  same  cut  of  the  potential  energy  surface 

given  in  Fig.  (1-1)  in  which  R  (the  distance  between  Li+  and  the  center- 

of-mass  of  N2)  is  infinite.    Planes  I  and  II  are  perpendicular  to  plane 

III  and  intersect  that  plane  at  r  equal  to  zero  and  at  r  equal  to  the 

equilibrium  internuclear  separation  (rg)  of  Ng.    It  is  shown  then  that 

4- 

the  ground  state  potential  energy  of  the  system  is  decreased  as  Li 
approaches  and  increases  rapidly  as  R  goes  to  zero.    The  reason  for  this 
initial  decrease  of  the  potential  as  the  Li+  approaches  N2  is  due  to 
attractive  polarization  forces.    This  will  be  explained  later  when  the 
long  range  interactions  of  the  system  are  discussed.    On  the  other  hand 
at  small  values  of  R,  strong  repulsion  forces  exist  because  of  the  closed 
shell  interactions.    A  similar  qualitative  behavior  is  expected  for  the 


Figure  (1-2)    Various  cuts  through  the  potential  energy  surface 
of  en  a torn- diatom  system  are  illustrated.  The 
parameter  R  measures  the  distance  of  the  projectile 
atom  to  the  center-of-mass  of  the  target  diatom; 
whereas,  r  is  the  internuclear  separation  of  the 
target  diatom.    Planes  I  and  II  are  perpendicular 
to  plane  III  at  r  equal  zero  and  r  ,  respectively. 
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excited  states.    It  should  be  mentioned,  however,  that  the  reaction 

would  be  spin-forbidden  (Ge71).    This  conclusion  is  only  valid  if  one 
can  neglect  effects  such  as  spin-orbit  interaction.    A  discussion  of 
the  applicability  of  symmetry  rules  to  collisions  is  given  in  the  work 
by  George  and  Ross  (Ge71)  and  that  by  Shuler  (Sh53).    From  an  experimenta' 
viewpoint,  electronic  excitation  has  not  been  observed  at  the  energies 
[4  to  8  eV  in  the  center-of-mass  system]  and  scattering  angles  studied 
in  the  present  work  (B676).    A  similar  analysis  could  be  done  for 
(Li+,  CO);  however,  this  system  is  isoelectronic  with  (Li+,  1^)  and  is 
expected  to  behave  similarly. 

We  have  not  considered  changes  of  the  equilibrium  internuclear 
separation  of  N?  as  Li+  approaches.    In  order  to  do  so,  we  now  discuss 
the  potential  energy  surface  in  more  detail.    The  long  range  interaction 
between  a  charged  atom  and  a  diatom  with  a  nondegenerate  electronic 
state  is  dominated  by  electrostatic  and  inductive  terms  (Bu67),  i.e. 


V(R,  r,  0)  -  VC0UL  +  VIND    ,  (1.58) 


wher 


VC0UL  *  j£  +  J  P2(COS9)  *  jj?  P3(a,50) 


nH 

+  an.  p.(cose)  +  .  .  .  .  (1.59) 

R  ' 


VIND  =  "  I  q2  "V  [a  +  I  (a||  "  a_L  ^cosQ)  +•••]> 
R 


(1.60) 
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a  =     (aj(  +  2a±  )  ,  (1.61) 

where  q  corresponds  to  the  charge  of  the  atom,  the  molecular  parameters 
Q,  Q  and  H  correspond  to  the  dipole,  quadrupole,  octapole  and  hexa- 
decapole  moments  of  the  diatom,  and  the  parameters        and  correspond 
respectively  to  the  parallel  and  perpendicular  components  of  the  diatomic 
polarizabi lity  relative  to  the  internuclear  axis  of  the  molecule.  The 
functions  P-|(cosO)  are  Legendre  polynomials  and  0  is  the  angle  between 
the  vector  ft  from  the  center-of-mass  of  the  diatom  to  the  projectile 
atom  and  the  internuclear  axis  of  the  diatom.    Beth  the  ^  and  CO  mole- 
cules have  very  similar  values  for  the  various  parameters  used  in  equa- 
tions (1.59)  and  (1.60),  with  the  exception  that  N~  does  not  have  a 
permanent  dipole  moment.    This  difference,  however,  is  quite  important 
since  the  charge-dipole  interaction  determines  the  R     term  in  Eq.  (1.59), 
and  as  such  dominates  the  long  range  portion  of  the  (Li+,  CO)  potential. 
It  is  precisely  these  parameters  [y,  Q,  etc.]  which  are  not  well  speci- 
fied in  SCF  calculations  [(St75),  (St76),  (Th78a)].    The  value  of  the 
dipole  moment  of  CO  in  fact  has  been  the  subject  of  wide  controvercy 
(Ne64).    Experimentally  this  value  is  small  and  negative  [u^q  =  -0.049 
a.u.]  (St76)  and  implies  the  polarity  of  the  CO  molecule  is  given  by 
C~0  .    More  recent  calculations  (Ch76)  have  now  shown  agreement  with  the 
experimental  results.    Unlike  the  situation  for  N^,  in  a  colli  near  reac- 
tion the  potential  between  Li+  and  CO  is  not  symmetric.    If  one  examines 
the  behavior  of  the  equilibrium  internuclear  separation  of  CO  as  a  func- 
tion of  R,one  would  find  that  it  shows  much  more  deviation  from  its 
value  at  R  equal  to  infinity  than  that  for         In  the  scattering  exper- 
iments on  these  systems  (B576),  it  is  found  that  (Li+,  CO)  shows  more 
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vibrational  excitation  than  the  (Li+,  No)  systems.    Various  investi- 
gators [(B676),  (St76),  (Mi78b)]  have  attributed  this  observed  result 
to  the  charge-dipole  interaction  in  the  (Li+,  CO)  system.    The  relative 
insensitivity  of  r   as  a  function  of  R  for  the  (Li+,  Ng)  system  is 
shown  in  Fig.  (1-3). 

At  small  values  of  R  the  interaction  of  Li+  with  N2  and  CO  is 
highly  repulsive  and  can  be  easily  characterized  by  an  exponential  poten- 
tial [St76),  (Th78a)].    Thomas  et  al .  (Th78a)  use  the  potential 

V(R,  r,  0)  =  Vj(r)  +  I  vx(r,  R)PA(cosG)  ,  (1.62) 

where  V,(r)  is  the  potential  for  the  free  diatom  and  vx(r,  R)  is  deter- 
mined by  using  various  parametrized  functions  [sums  over  exponentials 
and  R"n  terms]  to  yield  the  best  curve  fit  to  the  potential.  Equation 
(1.62)  behaves  asymptotically  as  Eq.  (1.58).    Staemmler  (St76)  points 
out  that  the  Legendre  expansion  of  the  potential  in  Eq.  (1.62)  loses 
its  significance  for  small  R[R  <  3.5  a.u.]  and  may  in  fact  not  converge. 

6.    Proposed  Application  of  the  Faddeev  Formalism 
Our  study  of  the  (Li+,  CO)  and  (Li+,  N2)  systems  will  be  based  on 
the  many  body  formalism  of  the  Faddeev  equations.    Within  this  approach 
we  shall  consider  an  iteration  of  the  Faddeev  equations  (in  transition 
operator  form)  which  yields  a  miltiple-collision  expansion  (Mi75). 
This  allows  us  to  describe  an  atom-diatom  collision  process  as  one 
which  takes  place  as  a  sequence  of  two-body  encounters.    Since  the  colli- 
sion energies  are  large  and  the  scattering  angles  are  beyond  the  rainbow 
angle,  it  is  expected  that  the  scattering  cross  section  will  be  largely 


Figure  (1-3)    Contour  plots  of  the  ground  state  potential 
energy  surface  for  the  (Li+,  N?)  system  as  a 
function  of  the  N-N  internuclear  separation 
(St75).    Each  curve  corresponds  to  a  fixed 
value  of  R,  the  distance  between  Li+  and 
the  center-of-mass  of  the  N2  molecule. 
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characterized  by  the  repulsive  forces  of  the  system.    Thus,  an  impulsive 
model  based  on  the  truncation  of  the  multiple-collision  expansion  would 
be  expected  to  yield  useful  information  on  these  systems.    An  advantage 
to  this  approach,  as  will  be  shown  in  Chapter  V,  is  that  one  is  able  to 
obtain  state-to-state  information  directly  from  the  initial  and  final 
momentum  distributions  of  the  target  and  product  diatom,  and  from  the 
various  two-body  potentials  involved.    Due  to  the  small  vibrational- 
rotational  spacings  of  CO  and  N2,  other  methods  such  as  the  coupled- 
channel  method,  based  on  Eq.  (1.23)  would  be  computationally  prohibitive 
(Th78b). 

We  shall  assume  in  our  implementation  of  the  Faddeev  equations  that 
electronic  motions  have  already  been  eliminated  and  that  only  information 
from  the  ground  adiabatic  potential  energy  surface  of  the  system  is 
needed.    As  mentioned  previously,  the  scattering  process  at  hyperthermal 
energies  probes  the  short  range  part  of  the  potential.    Therefore,  we 
shall  develop  a  many  body  theory  based  on  the  decomposition  of  the  adia- 
batic potential  into  the  form  (Mi 79 ) 

V(3,  r)  =  Vsr  +  Vlr  (1.63) 

where 

3 

V.    =    I    V.(r.)  (1.64) 
sr     i=l    1  1 

and 

Vlr  =  I  V>(R,  r)Yx(R,  r)  .  (1.65) 
A 

The  V-jr  term  denotes  the  long  range  potential  and  is  analogous  to  the 
Legendre  expansion  used  in  Eq.  (1.62).    Here,  however,  Y. (R,  r)  is  a 
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general  function  depending  only  on  the  relative  orientation  of  the  atom 
and  diatom.    The  Vgr  term,  on  the  other  hand,  corresponds  to  the  short 
range  potential  and  is  characterized  by  a  sum  of  many-body  terms.  In 
Eq.  (1.54),  V     was  assumed  to  be  a  sum  of  radial  two-body  potentials; 
one  potential  for  each  two-atom  fragment  in  the  atom-diatom  system. 
This  breakdown  of  the  potential  circumvents  the  difficulties,  mentioned 
in  the  last  section,  in  representing  the  short  range  behavior  of  V  by 
Eq.  (1.62).    Having  introduced  Eq.  (1.63),  we  propose  an  impulse  model 
in  which  only  the  short  range  part  of  the  potential  is  kept.    The  re- 
sulting problem  is  then  ideally  suited  for  the  many-body  formalism  of 
the  Faddeev  equations. 

We  note  that,  at  low  collision  energies,  one  expects  that  the 
approach  of  Li+  to  CO  or  ^  would  heavily  influence  the  orientation  of 
the  diatom.       This  would  be  particularly  true  for  a  molecule  with  a 
dipole  moment  like  CO.    The  net  result  in  our  impulsive  model  would  be 
to  overestimate  rotational  excitation.    At  the  high  energies  considered 
in  this  study,  however,  this  effect  is  expected  to  be  minimal. 

7.    Plan  of  Dissertation 
The  present  investigation  of  atom-diatom  collisions  will  be  given 
in  three  major  parts.    The  first  part  presents  an  analysis  of  the  two- 
body  problem  and  discusses  various  computational  schemes  for  the  deter- 
mination of  off-shell  matrix  elements  of  the  two-body  transition  opera- 
tor (Chapters  II  and  III).    In  the  next  part  (Chapters  IV  and  V),  the 
formalism  of  the  three-body  problem  is  presented  and  the  various  approx- 
imations obtained  from  the  multiple  collision  expansion  are  discussed. 
Finally,  in  the  third  part  (Chapters  VI  and  VII),  computational  results 


based  on  the  single-collision  approximation  are  presented  for  both 
collinear  and  three-dimensional  atom-diatom  scattering  processes. 


CHAPTER  II 
THE  TWO- BODY  t-MATRIX 

Scattering  processes  in  few-  or  many-body  systems  require  the 
introduction  of  two-body  transition  operators,  whose  matrix  elements 
must  be  known  for  arbitrary  initial  and  final  relative  momenta  of  the 
two  bodies  and  for  arbitrary  energies,  i.e.,  they  must  be  known  in 
general  "off-the-energy-shell ."    In  other  words,  the  operator  quantities 
of  interest  are  those  describing  the  interaction  of  two  particles 
embedded  in  a  many  particle  system,  and  as  a  consequence  energy  is  not 
conserved  locally,  thus  explaining  the  term  "off-shell".    Although  these 
two-body  transition  operators  stand  for  basic  physical  concepts,  their 
properties  and  values  are  little  known  for  molecular  systems,  where  the 
bodies  are  atoms  or  ions. 

The  computation  of  two-body  t-matrix  elements  has  long  been  a 
problem  of  intense  interest  in  the  field  of  nuclear  physics.  Conse- 
quently, there  exists  an  extensive  literature  devoted  to  this  subject 
as  it  pertains  to  potentials  describing  nucleon-nucleon  interactions. 
A  very  nice  and  extensive  review  in  this  regard  is  the  work  of  Srivastava 
and  Sprung  (Sr75).    In  the  chemical  literature,  work  in  this  area  has 
been  sparse.    However,  there  have  now  appeared  various  papers  which  deal 
with  the  computation  of  t-matrix  elements  for  potentials  of  chemical 
interest,  for  example:    van  Leeuwen  and  Reiner  (vL61)  proposed  a  numer- 
ical procedure  based  on  the  comparison  potential  approach,  which  has  now 
been  implemented  in  the  work  of  Brumer  and  Shapiro  (Br75);  Kuruoglu  and 
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Micha  (Ku78)  have  developed  a  variational  approach  leading  to  a 
separable  approximation  to  the  t-matrix;  and  along  more  approximate 
methods,  Korsch  and  MShlenkamp  (Ko77)  have  developed  a  semiclassical 
approach  based  on  the  JWKB  approximation  which  is  applicable  to  repulsive 
potentials. 

The  principal  aim  of  this  chapter  will  be  to  present  a  new  computa- 
tional method  for  obtaining  matrix  elements  of  two-atom  transition  oper- 
ators between  arbitrary  momentum  states,  of  high  accuracy,  and  applicable 
to  any  radial  diatomic  potential.    The  approach  we  have  taken  is  based 
on  equations  presented  some  time  ago  (vL61),  which  reduce  the  calculation 
of  transition  operator  matrix  elements  to  integration  of  inhomogeneous 
radial  differential  equations  followed  by  numerical  quadratures.  Our 
procedure  is  based  on  the  well-known  variable-phase  approach  [(Ca63), 
(Ca67)],  but  extended  to  the  present  inhomogeneous  equations.  This 
extension  shows  that,  unlike  the  situation  for  standard  radial  scatter- 
ing equations,  the  variable  phase  is  coupled  at  every  distance  with  a 
variable  amplitude.    To  emphasize  this  point  v/e  refer  here  to  a 
"Variable-Phase  and  -Amplitude"  approach.    Another  feature  of  our 
approach  is  the  way  in  which  one  can  treat  propagation  through  classi- 
cally forbidden  regions.    Along  with  the  development  of  the  VPA  method 
we  shall  present  some  suggested  improvements  to  the  comparison  potential 
method. 

All  the  methods  we  have  mentioned  above    make  use  of  a  partial 
wave  expansion  to  construct  the  two-body  t-matrix  elements.    When  we 
discuss  the  single-collision  approximation  in  Chapter  V,  it  will  become 
apparent  that  it  would  be  desirable  to  have  a  procedure  for  the  direct 
computation  of  these  matrix  elements.    This  will  even  become  crucial 
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when  one  deals  with  physical  systems  at  energies  for  which  a  large  number 
of  partial  waves  are  required  and  the  computational  expense  becomes  im- 
portant.   Consequently,  in  light  of  this  problem,  we  shall  also  investi- 
gate a  non-partial  wave  approach  that  has  been  recently  developed  by 
Belyaev  and  coworkers  (Be74),  which  is  based  on  the  Bateman  method 
[(Ba22),  (Ka58)]  for  solving  integral  equations. 

1.    General  Two-Body  Scattering  Theory 
The  general  computation  of  off-shell  matrix  elements  of  the  two-body 
t-operator  for  simple  non-singular  potentials  can  in  general  be  solved 
without  much  difficulty.    This  may  be  done  by  direct  consideration  of 
the  Lippmann-Schwinger  equation  [(Si71),  (Sr68)] 

t(z)  =  V  +  VGo(z)t(z),  '  (2.1) 

where 

GQ(z)  =  (z-H^"1.  (2.2) 
In  these  relations  HQ  is  the  kinetic  energy  operator  for  relative  motion, 
z  is  a  complex  energy  parameter  (z  =  E+ie)  and  V  Is  the  interaction 
potential  between  two  particles.    V  will  always  foe  assumed  local  and 
spherically  symmetric.    From  equation  (2.1)  it  is  apparent  that  the 
problem  must  be  reformulated  for  singular  potentials,  those  which  contain 
a  hard  core  or  have  a  singularity  stronger  that  /rc  at  the  origin 
(Ta72).    This  can  be  accomplished  by  introducing  the  wave  operator  W(z) 
via 

t(z)  =  VW(z),  (2.3) 

where 

W(z)  =  1  +  G0(z)VW(z),  (2.4) 
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and  by  relating  the  matrix  elements  of  the  t-operator  to  those  of  the 
wave  operator  in  a  mixed  coordinate  momentum  representation  (vL61). 
Multiplying  (2.4)  by  GqX(z)  on  the  left,  one  finds  that  W(z)  also 
satisfies 

(z  -  HQ  -  V)W(z)  =  (z  -  HQ).  (2.5) 

In  a  mixed  coordinate-momentum  representation  one  then  obtains  the 
inhomogeneous  Schrbdinger  equation  below  (ti  =  1): 

[z  +  (2m)"1V2  -  V(r)]<r|o^(z)>  =  [z  -  (2m)-1q2]<r|q>,  (2.6) 

where  m,  f  and  q  are  respectively,  the  reduced  mass,  the  relative 
coordinate  and  the  momentum  of  the  two  particles  being  considered,  and 
!<jj^(z)>  is  the  off-shell  state  (since  q2/2m  f  E)  defined  as 

|u*(z)>  =  W(z)|q>  .  (2.7) 

In  the  nuclear  physics  literature,  (2.6)  is  known  as  the  Bethe-Goldstone 
equation  (Sr68).    A  reduced  radial  equation  may  now  be  obtained  by 
introducing  into  equation  (2.6)  the  following  partial  wave  decompositions 

<i*|q%  =  (2tt)"3/2  e^'r 

=  (2tt)"3/2    E   i1  (21  +  l)P1(5-r)j1(qr)/(qr)  (2.8) 

and 

<r|c^(z)>  =  (2tt)"3/2    I  i1  (21  +  l)P1(5-r)Wl(q,kE;r)/(qr),  (2.9) 

where  j-,  (qr)  is  the  regular  Riccati-Bessel  function  and  k£  =  \2mt\1^2. 
The  resulting  reduced  radial  equation  is 
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t-2-  +  o.k2  -  1*1^1  -  v(r)JuWq,k_;r)  ■  (a  k*  -  q2)j_(qr),  (2.10) 
dr        E  E  r£  If  EE  1 

where      is  +1  if  E  is  positive  or  -1  otherwise,  and  v(r)  =  2mV(r). 
The  matrix  elements  of  the  t-operator  can  now  be  obtained  from  those  of 
the  wave  operator  by  using  equation  (2.3).    Thus  one  has  the  relation 

t(q',  q;  z)  =     /  <q'  | r>  V(r)  <?|<^(z)>dr\  (2.11) 

which  may  also  be  decomposed  into  partial  waves  with  the  result 

t(q\  q;  z)  =  (47T)"1     I    (21  +  ljP^q' -q)^  (q1 ,  q;  E),  (2.12) 

where 

tl(q''  q;  E)  =   f  q^     '    jVq,r)V(r)aJl(q'  kE;  r)  dr-  (2*13) 

This  definition  of  the  partial  wave  matrix  elements  t-j(q',  q;  E)  follows 
directly  from  equations  (2.8),  (2.9)  and  (2.12). 

2.    Boundary  Conditions  For  The  Off-Shell  Wave  Function 
In  order  to  obtain  the  off-shell  wave  function  w^(q,  k^;  r),  a 
solution  to  equation  (2.10)  must  be  found  subject  to  two  boundary  condi- 
tions (vL61): 

Ul(q,  kE;  r)  p  ~  q  o  (2.14) 

and 

^(q,  k  ;  r)     ~     j  (qr)  +  C  h(+)(k  r),  (2.15) 

y  -*■  oo  lit 

where      is  a  constant  and  hj+)  is  a  Riccati-Hankel  function  (Ca67). 
Thus,  u-|(q,  k£;  r)  is  a  complex  wave  function  which  asymptotically 
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behaves  as  the  sum  of  a  component  from  a  plane  wave  and  that  of  an  out- 
going spherical  wave,  which,  unlike  on-shell  scattering  expressions, 
depends  on  different  momenta  (Ta72).    For  potentials  containing  a  hard 
core  (2.14)  must  be  replaced  by 

W,(q,  k  ;  r)     ~       0  (2.16) 
r  ->  r 

c 


where  r   is  the  hard  core  radius.    Once  the  off-shell  wave  function  is 
c 

obtained,  the  partial  wave  matrix  elements  t^(q',  q;  E)  follow  directly 
from  (2.13).    However,  if  the  potential  contains  a  hard  core,  care  must 
be  taken  in  applying  Eq.  (2.13).    If  one  integrates  equation  (2.10)  as 
fol 1 ows , 


J  K  +  cEkE  u^]  dr  =  /  {[v(r)  +  iLL±_U]  ^  + 
rc-e  rc-e  r2 


+  (o£k2£  -  q2)j1(qr)}  dr,  (2.17) 

it  is  apparent  that  the  product  co^(q,  k^;  r)v(r)  does  not  vanish  inside 
the  core  region  even  though  w^(q,  k  ;  r)  does.    In  fact,  for  a  potential 
containing  a  hard  core  (La68), 

v(r)Ul(q,  kE;  r)  =  6(r  -  rQ)  ^^(q,  k£;  r,+)  + 

(q2  -  aEk2)  j1  (qr)  (2.18) 

for  r  <  r  ,  where 
—  c 

Hr  kE>  rc+)  =     1im      clr  Vq'  kE;  rc  +  e)  '  (2,19) 

e  ■*  o+ 
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A  derivation  for  the  case  of  a  Dure  hard  core  potential  is  given  in 
Appendix  I. 

3.    The  Variable-Phase  And  -Amplitude  Method 
The  aim  of  the  VPA  method  is  to  take  advantage  of  the  known 
oscillatory  nature  of  the  off-shell  wave  function,  w^(q,  k£;  r).  This 
is  done  by  considering  a  solution  to  Eq.  (2.10)  which  has  the  form 
below: 

o)-|(q,  kE;  r)  =  j^qr)  +  ^(q,  k£;  r)  cos  [fi^q,  k£;  rjju^^r) 

-  a^q,  kE;  r)  sin  [6}{q,  k£;  rjjv^r),  (2.20) 

where  u^(k^r)  and  v^(k^r)  are  Riccati-Bessel  functions  (j^k^r)  and 
n^(k^r))  or  Riccati-Hankel  functions  (k^(k^r)  and  i^k^r))  (Ca67).  The 
choice  is  determined  by  the  nature  of  the  region  through  which  one  is 
propagating,  i.e.,  whether  the  region  is  classically  allowed  or  not. 
Using  equations  (2.10)  and  (2.20),  one  may  obtain  a  set  of  coupled 
equations  for  the  amplitude,  a^(q,  k£;  r),  and  the  phase,  5^(q,  k  ;  r). 
These  equations  will  now  be  derived  using  an  integral  equations  approach 
(Ro67).    The  functions  u^k^r)  and  v^(k^r)  used  in  expression  (2.20) 
satisfy  the  Riccati-Bessel  differential  equations  (Ab65) 

d 


5   +  [E  -  ^Ap-]  z,  =  0  .  (2.21) 
dr*    1  r  1 

Taking  the  regular  and  irregular  functions  u-j  and  v-|  to  be  normalized 
such  that  their  Wronskian,  Wfu^,  v,]s  is  unity,  the  Green's  function 
associated  with  (2.21)  may  be  written  as 

gl(r'  r<)  =  ¥E  VkEr<}  ^1  (kEr>)  *  (2,22) 
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where  r<  is  the  lesser  of  r  and  r',  and  ry  is  the  greater  of  r  and  r' 
(Ro67).    At  this  point,  the  discussion  will  be  restricted  to  positive 
energies,  since  this  would  not  involve  a  loss  of  generality.    It  is 
readily  seen  that  equation  (2.21)  is  just  the  corresponding  differential 
equation  to  (2.10)  in  the  absence  of  a  potential.    Thus,  o^(q,  k£;  r) 
may  be  written  as  follows: 

00 

^(q,  k£;  r)  =  N'j^r)  +     /  dr'g^r,  r')[(k|  -  q2)  j^qr')  + 

v(r')Wl(q,  k£;  r')]  ,  (2.23) 

where  N'  is  a  constant.    Since  j^qr)  satisfies  equation  (2.21),  the 
first  term  in  the  above  integral, 

(k2  -  q2)     /o°°  dr'g^r,  r'tf^qr')  , 

may  be  evaluated  exactly.  This  may  be  accomplished  using  the  relation 
(A2.2  )  given  in  Appendix  II.  Introducing  the  quantities  A^(q,  k^;  r) 
and  B^(q,  k^;  r)  defined  by 

A^q,  kE;  r)  =  N  +  k"1    J   dr'v(r'  )^  (q,  k£;  r')  n^^r')  (2.24) 

and 

B^q,  kE;  r)  =  -k'1     /[  dr,v(r'  )(*,  (q,  k£;  r' )  j^r'),  (2.25) 
where  N  =  N'  +  1,  and  using  equation  (2.22),  one  obtains 
(^(q,  k£;  r)  =  j^qr)  +  A^q,  k£;  r)  j]  (kjrj  + 

-  B^q,  kE;  r)  n]  (k^rj  .  (2.26) 
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The  functions  A^  and  B,  behave  as  constants  asymptotically,  as  can  be 
seen  from  the  boundary  conditions  (2.14)  and  (2.15).    These  constants 
may  in  turn  be  related  to  an  off-shell  phase  shift,  6^(q,  k£),  and 
amplitude,  a-j(q,  k^),  in  the  same  manner  as  is  done  for  the  on-shell 
case  [(Ca63),  (Ca67) ] .    Therefore  one  has 

A,(q,  k  ;  r)     ~     a  (q,  k  )  cos  [6  (q,  k  )]  (2.27) 

and 


B,(q,  kE;  r)     ~     a  (q,  k£)  sin  [6  (q,  kJ]  .  (2.28) 

p  -y  oo 

The  variable  phase  assumption  consists  of  extending  the  definition  of 
the  phase  and  amplitude  to  all  values  of  r  (Ca63);  thus, 

A^q,  k£;  r)  =  ^(q,  kE;  r)  cos  ^(q,  kE;  r)]  (2.29) 

and 

B^'q,  k£;  r)  =  ^(q,  k£;  r)  sin  ^(q,  k£;  r)]  .  (2.30) 

This  assumption  is  of  course  valid,  since  6-j  (q ,  k£;  r^)  and  a^(q,  k_;  r^) 
are  just  the  resulting  phase  shift  and  amplitude  for  a  potential  with 
a  cut  off  at  r  =  rM  (Ca67).    In  fact  the  justification  of  this  approach 
is  best  understood  by  considering  a  finite  boundary  value  problem  and 
extending  the  boundary  to  infinity  (Ze64).    The  desired  set  of  coupled 
first  order  differential  equations,  the  VPA  equations,  may  now  be 
derived  using  (2.29)  and  (2.30)  for  A]  and  B1 .    To  accomplish  this, 
first  consider  the  derivatives  of  A^  and      as  defined  by  equation 
(2.24)  and  (2.25).    From  these  relations  it  is  apparent  that 

VkEr)  dV  Al(q'  kE;  r)  "  VkEr)  ^  Bl(q'  *E;  r)  =  °  ;  (2,31) 


therefore,  the  derivative  of  the  off-shell  wave  function  has  no  terms 
involving  the  derivatives  of  A-|  and  B^,  i.e., 

3p  wj(q,  kE;  r)  =  ^-  J^qr)  +  A^q,  kE;  r)  J1(kEr) 


B^q,  kE;  r)  ^  n^r)  .  (2.32) 


Making  use  of  (2.29)  and  (2.30),  equations  (2.26)  and  (2.32)  may  be 
written  as, 

w-,  ■  J-,  (qr)  +  a-,  cos  (6-j)  j-,  (k£r) 

-  a1  sin  (6^)  n-,  (k£r)  (2.33) 

and 

dFui  =  "d?  JVqr)  +  al  cos  ^  dV  jl(kEr) 

-  a-,  sin  (6-,)  ^  n^r)  .  (2.34) 

Taking  the  derivative  of  (2.33)  and  comparing  it  with  (2.34)  one  obtains, 

d_a  _a  sin  [^j  jl  <kEr)  +  cos  (gl)  jjl  (kEr)  d_  .  (2.35) 
dr    1       1     cos  (6,)  f,  (kEr)  -  sin  (6,)  ^  (k£r)    dr  1 

In  a  similar  manner  taking  the  derivative  of  (2.34)  and  comparing  it  to 
the  inhomogeneous  Schrbdinger  equation  (2.10),  one  finally  obtains,  with 
the  aid  of  (2.35),  the  "variable-phase  and  -amplitude"  equations 

61(q,  kE;  r)  =  -k"1  v(r)[u>*  -  jyqrja)]/^  (2.36) 

and 

^^(q,  k£;  r)  =  -k"1  vfr)^  [sin(6] ) j]  (k£r)  +  cos(«1  )n1  (k  r)]. 

(2.37) 
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These  equations  are  analogous  to  those  for  the  on-shell  case,  where 
q  =  k£  (Ca63).    However,  for  the  on-shell  case  the  phase  and  amplitude 
equations  are  uncoupled. 

The  boundary  condition  satisfied  by  the  VPA  equations  at  the 
origin  follows  directly  from  equations  (2.14)  and  (2.33).    From  these 
expressions  it  is  apparent  that 

sin  [6,(q,  kE;  r)]     ~     0  (2.38) 

r  ->  o 

and  that  the  value  of  a-j(q,  k^;  r)  for  the  same  limit  is  arbitrary,  and 
becomes  fixed  only  after  the  asymptotic  boundary  condition  (2.15)  is 
applied.    On  the  other  hand,  if  a  hard  core  is  present 

o  (q,  k  .  r)     -  I  Vqrc]  x  (2.39) 

1        E       r  ■+  rc  cosCfi^^  (k£rc)  -  sin(6]  (kErc) 

and  6-|(q,  k£;  r  )  is  arbitrary. 

4.    Computational  Aspects  of  the  VPA  Equations 
The  procedure  used  in  calculating  w. (q,  k£;  r)  consists  of  propa- 
gating two  solutions  to  the  VPA  equations  (2.36)  and  (2.37)  such  that 
they  satisfy  the  initial  boundary  condition,  and  constructing  linear 
combinations  of  these  solutions  so  as  to  satisfy  the  asymptotic  boundary 
condition  (2.15).    It  is  possible  to  construct  two  linearly  independent 
solutions  that  satisfy  the  initial  boundary  condition  due  to  the 
inhomogeneous  character  of  equation  (2.10).    This  procedure  is  analogous 
to  that  found  in  reference  (Br75);  however,  instead  of  a  piece-wise 
breakdown  of  the  potential  and  the  propagation  of  local  solutions  as 
given  there,  numerical  integration  is  involved  here. 


An  important  consideration  is  the  singular  nature  of  the  VPA 
equations  near  the  origin.    Equations  (2.36)  and  (2.37)  contain  the 
potential  as  a  multiplicative  factor.    This  problem  can  be  overcome  by 
assuming  hard  core  boundary  conditions  or  by  using  a  simple  step  poten- 
tial within  the  first  classical  turning  point  [(Be60),  (Br75)].  The 
second  of  these  procedures  was  chosen  in  this  work.    The  length  of  the 
first  step,  r^,  and  the  value  of  the  assumed  step  potential  over  this 
first  step,  V^,  were  chosen  such  that  making  any  further  optimal  changes 
in  these  parameters  would  not  lead  to  a  significant  change  in  the  value 
of  the  calculated  matrix  elements,  w, (q,  k£;  r)  and  t^(q',  q;  E).  The 
guess  work  involved  here  can  be  eliminated  by  using  a  JWKB  approximation 
to  the  derivative  of  the  imaginary  part  of  the  wave  function  near  the 
first  classical  turning  point.    The  approximations  involved  are  given  in 
Appendix  III.    From  expression  (2.10)  it  is  apparent  that  Im[(^(q,  kc;  r)] 
satisfies  the  usual  Schrb'dinger  equation  and  boundary  condition  at  the 
origin.    It  therefore  follows  that  within  the  first  classical  turning 
point  Im[w-j]  is  decaying  exponentially,  and  is  zero  at  the  origin.  All 
that  is  necessary  then,  from  a  computational  point  of  view,  is  to  start 
the  integration  in  a  region  where  the  magnitude  of  the  derivative  of 

Im[w-j]  is  smaller  than  some  given  value.    In  our  calculations  a  value  of 
-  i  0 

10  AU  was  found  to  be  adequate.    Another  important  consideration  in 
propagating  the  VPA  equations  is  the  cost  involved  in  evaluating  a 
large  number  of  Bessel  functions.    Because  of  this,  the  coupled  equations 
were  integrated  numerically  using  an  Adams  numerical  integration  routine 
developed  by  Shampine  and  Gordon,  which.is  particularly  suited  for 
problems  for  which  derivative  evaluations  are  expensive.    These  codes 
and  a  discussion  of  their  efficiency  and  accuracy  are  found  in  reference 
(Sh75). 


The  propagation  of  two  independent  solutions  to  the  inhomogeneous 
differential  equation  (2.10)  further  presents  two  problems.    One  is  that 
of  the  growth  rate  of  the  amplitude  in  a  non-classical  region  and  the 
second  is  that  of  the  propagation  through  classically  allowed  regions 
at  negative  energies.    The  first  problem  was  solved  as  given  in  reference 
(Br75).    The  idea  was  simply  to  use  the  propagated  wave  functions  to 
construct  two  new  linearly  independent  solutions  to  Eq.  (2.10),  which 
would  have  much  smaller  amplitudes.    Therefore  we  shall  restrict  our 
attention  to  the  second  problem.    If  the  energy  is  negative,  the  modified 
Riccati-Hankel  functions      and  i-|  are  the  natural  basis  to  use  for 
u-|  and  v.j  in  expression  (2.20),.  the  VPA  form  for  the  off-shell  wave 
function  a, .    However,  if  the  potential  is  attractive  and  the  energy  is 
such  that  one  must  integrate  through  a  classically  allowed  region,  the 
wave  function  in  this  region  would  be  oscillatory  in  nature.  The 
Riccati-Hankel  functions  on  the  other  hand  have  a  non-exponential  grow- 
ing and  decaying  character;  thus,  propagation  of  the  off-shell  wave 
function  under  these  circumstances  is  very  difficult.    The  solution  then 
is  simply  to  replace  k^  and  ?1  by  the  Riccati-Bessel  functions  J   and  ry 
This  in  turn  requires  a  modification  of  the  VPA  equations  (2.36)  and 
(2.37).    However,  the  modification  is  minor  and  does  not  lead  to  any 
significant  increase  in  computational  effort.    Let  us  go  back  then  to  the 
original  inhomogeneous  Schrbdinger  equation,  Eq.  (2.10),  as  it  applies 
to  negative  energies,  and  rearrange  it  into  the  form 

kE  "  — ~z  v<r)h  =  2k^  -  (k2  +  q2)jl(qr).  (2.40) 


Treating  the  left  hand  side  of  this  equation  as  an  inhomogeneous  term, 
we  can  adopt  the  same  procedure  that  was  used  in  section    I I - 3  to 
derive  the  desired  set  of  modified  VPA  equations 

^  S^q,  kE;  r)  =  -k"1  {v(r)[<o*  -  ^(qr)]  (2.41) 

+  2k^  [to,  -  j^qr)]}/^ 

and 

^^(q,  k£;  r)  =  k"1  {v(rV,  +  2k\        -  Jjfar)]}  x 

x  {$111(6^ (kEr)  +  cos(61)n1(kEr)}.  (2.42) 

A  similar  set  of  modified  equations  could  be  used  for  positive  energies 
when  a  nonclassical  region  is  encountered  there.    This  would  require 
Riccati-Hankel  functions.    This,  however,  was  not  found  to  be  necessary 
in  the  calculations  done  in  this  work.    Stabilization  was  only  a  problem 
at  negative  energies. 

5.    The  Comparison  Potential  Method 
Another  possible  approach  in  solving  the  inhomogeneous  Schrodinger 
equation  (2.10)  for  the  off-shell  wave  function  w (q,  k^;  r),  is  to  make 
use  of  the  comparison  potential  (CP)  method.    As  was  mentioned  in  the 
introduction  to  this  chapter,  this  approach  has  already  been  presented 
in  the  literature  [(vL61),  (Br75)].    We  shall,  however,  sketch  an  outline 
of  the  method  and  present  some  improvements  to  the  computational 
algorithm  for  u,.    These  changes  are  based  on  the  work  of  R.  Gordon 
t(Go69),  (Go71),  (Ro76)],  who  has  investigated  extensively  the  use  of 


local  approximations  to  the  potential  in  order  to  simplify  scattering 
calculations.    The  outcome  of  this  pursuit  will  then  be  helpful  in 
suggesting  some  improvements  to  the  VPA  method,  or  for  that  matter  any 
method  which  implements  a  numerical  solution  for  the  off-shell  wave 
function  . 

The  basic  idea  behind  the  CP  method  is  to  approximate  the  potential 
on  a  given  interval  in  such  a  way  that  one  can  express  the  solution  to 
the  resulting  differential  equation  in  terms  of  known  functions  which 
present  no  problem  computationally.    To  meet  these  requirements  the 
potential  has  often  been  treated  locally  as  a  step  potential  or  as  a 
linearly  varying  potential  (Go71).    The  centrifugal  term,  1(1  +  l)/r  , 
may  also  be  included  in  the  potential,  and  the  total  effective  potential 
may  be  approximated  (Go71).    In  the  present  work,  however,  only  the 
potential  was  approximated,  taking  it  to  be  a  collection  of  step- 
potentials.    Thus  one  has 


where      plays  a  role  analogous  to  a£  and  v..  is  tfte  local  approximation 
to  v(r)  for  r.  <.  r  <  r^.    The  advantage  of  this  particular  choice  lies 
in  the  fact  that  one  can  use  Bessel  functions  as  tthe  local  basis  from 
which  to  construct  the  solution  to  (2.43);  and  sinree  the  centrifugal 
term  is  treated  exactly,  the  convergence  toward  lite  asymptotic  boundary 
condition  (2.15)  is  accelerated  in  the  case  of  the  higher  partial  waves 
(An76).    Locally  then,  on  the  itn  interval 


(2.43) 


and 


(2.44) 
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wj  (q,  kE;  r)  =  A]  u}  (ic.r)  -  ej  v.  (ic.r)  +  c]  j,  (qr),  (2.45) 

where 

i        ackr  -  qc 
aiS-  -  q 

The  off-shell  wave  function  is  then  obtained  from  two  propagated  solu- 
tions in  a  form  analogous  to  that  given  for  the  VPA  approach.    The  prop- 
agation here,  however,  is  accomplished  by  matching  each  wave  function 
and  its  derivative  across  every  interval.    This  can  be  done  by  either 
propagating  the  coefficients  A^j  and      as  v/as  done  in  the  reference 
(Br75),  or  by  propagating  the  wave  function  and  its  derivative  directly 
as  we  propose  to  do.    To  accomplish  this,  we  shall  modify  the  procedure 
given  by  Rosenthal  and  Gordon  (Ro76)  for  the  case  of  the  usual  radial 
Schrbdinger  equation.    From  expression  (2.45)  and  the  corresponding 
expression  for  the  derivative  of  the  wave  function  on  the  given  i^ 
interval,  it  is  easy  to  show  that 

Aj    =   Wtco^r.)  -  C]  j1  (qr.),  I^r.^/W.  (2.47) 

and 

B]    =   Wto^r.)  -  cj  ^  (q,  r.),  u}  (K.r.)]/W.  ,  (2.48) 

where  W..  is  a  constant  equal  to  the  value  of  the  Wronkian  of  u-j  and  v^ , 
the  two  local  basis  functions.    The  i  superscript  on  co-j  has  been  removed 
and  is  not  necessary  because  of  the  continuity  requirement  imposed  on 
the  wave  function  and  its  derivative.    Using  (2.47)  and  (2.48)  in  the 
expression  for  u>^,  Eq.  (2.45),  and  the  analogous  expression 
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j  (q.  kE;  r)  =  a|      u}  [*jf)J  +  b]  [fc  v]  (^r)] 


dr 


(X) 


+  cj  (qr)] 


(2.49) 


for  its  derivative;  the  desired  result 


=   C-fwCr.)  -  C]j(qr.)]    +  c]j  (qr.+1) 


(2.50) 


is  obtained,  where  we  have  introduced  the  matrix  definitions 


u(r)  = 


uJ  . 


dr  1 


(2.51) 


and 


i(r)  = 


j-j  (qr) 
^J"l  (qr) 


(2.52) 


and  the  coefficient  matrix  C ,  whose  components  are  defined  by  the 
relations 


C11^W  -^W  tiF"^i)I>/Wi 


(2.53) 


C12  E  tv(r.xl)  u  (r.)  -  u  (r,4.J  v  (r,)}/W 


i+1 


i+1' 


(2.54) 


'21 


^HF^Hl)!^^)]  -  ^W^^H^^r.)]}^.  (2.55) 


and 


(2.56) 


^9 


In  order  to  simplify  the  notation,  we  have  suppressed  V."  from  the 
argument  of  u^  and  v-j ,  as  well  as  the  "1"  subscript.    Equation  (2.49) 
is  appealing  in  the  sense  that  it  provides  a  simple  algorithm  for  propa- 
gating w-i  and  its  derivative  from  r^  to  r-j+p  but  it  is  even  more 
appealing  in  the  sense  that  one  need  not  store  the  coefficients  A^ ,  Bj 
and  Cj  for  all  the  intervals  needed  in  the  calculation.    Only  one 
problem  remains  to  be  dealt  with.    In  the  CP  method  the  t-matrix  is 
calculated  replacing  (2.13)  by  the  expression 

t^q',  q;  E)  =  £ttt     ?     V1  /     j^q'rjw  (q,  k£;  r)  .  (2.57) 

i  r.j 

If  one  had  propagated  the  coefficients  A^j ,  B^J  and  Cy  on  each  interval 
one  would  have  a  set  of  integrals  over  products  of  two  Riccati-functions 
as  is  indicated  by  Eq.  (2.45),  the  expression  for  (J,  and  Eq.  (2.57), 
the  expression  for  t^ .    These  integrals  are  well  known  and  are  easy  to 
generate  (Me61).    We,  however,  have  been  able  to  derive  expressions  for 
the  integrals  needed  in  Eq.  (2.57)  directly  in  terms  of  the  off-shell 
wave  function  oi,  and  its  derivative.    A  derivation  of  these  quantities 
is  given  in  Appendix  II. 

One  of  the  major  criticisms  of  the  VPA  approach  and  other  numerical 
procedures  for  obtaining  w,,  is  that  t^q',  q;  E)  must  be  obtained 
through  numerical  quadratures.    In  practice  this  is  problematic  when 
q  and  q'  are  larger,  because  of  the  oscillatory  nature  of  the  functions 
involved.    Generally,  when  one  propagates  w-j  numerically,  it  is  deter- 
mined at  smaller  step  sizes  than  are  involved  in  the  CP  approach;  thus, 
one  could  take  advantage  of  the  integral  expressions  mentioned  above  by 
implementing  a  breakdown  of  the  potential  into  step-potentials  only  in 
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the  evaluation  of  t^q',  q;  E).    This  would  then  remove  any  problems  in 
this  regard. 

6.    The  Ba tenia n  Method 
Very  little  has  been  done  in  the  literature  with  regard  to  non- 
partial  wave  techniques  in  the  evaluation  of  scattering  amplitudes  or 
two- body  t-matrix  elements.    Only  within  the  Eikonal  approximation 
[(6159),  (Su69),  (Ch73)]  and  approximations  based  on  the  Born  expansion 
[(Ho68),  (Ra72)]  has  this  not  been  the  case.    The  attempts  that  have 
been  made  to  solve  the  Lippmann-Schwinger  equation  for  either  the  full 
wave  function  or  the  t-matrix  have  involved  numerical  quadratures.  For 
example,  Walters  (Wa71)  has  considered  the  Lippmann-Schwinger  equation 
for  the  full  wave  function  and  converted  it  to  matrix  form  by  use  of  a 
two  dimensional  numerical  quadrature  for  the  case  of  a  screened  Coulomb 
potential.    The  amplitudes  obtained  from  this  procedure  however,  though 
in  good  agreement  with  exact  results,  showed  poor  convergence  when  the 
number  of  quadrature  points  was  increased.    Another  approach  that  has 
been  taken  is  that  of  Rosenthal  and  Kouri  (Ro73),  where  the  Lippmann- 
Schwinger  equation  for  the  scattering  amplitude  was  reduced  from  its 
two-dimensional  form  to  a  one-dimensional  integral  equation,  which  was 
then  solved  numerically.    Using  the  work  of  Walters  (Wa71)  for  comparison, 
these  authors  obtained  results  within  5  to  10  per  cent  agreement.  In 
contrast  to  these  methods,  the  approach  we  shall  snow  pursue  does  not 
involve  numerical  quadratures,  at  least  in  the  case  of  certain  model 
potentials.    This  procedure,  known  as  the  Batemam  method,  was  originally 
implemented  by  Akhmadkhodzhaev,  Belyaev  and  Wrzeoonko  (Ak70)  in  the 
computation  of  partial  wave  t-matrix  elements.    Ttey  have,  however, 
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extended  this  work  to  the  calculation  of  the  full  t-matrix  elements  for 
the  case  of  a  Gaussian  potential  (Be74a),  with  significant  success.  Lim 
and  Gianini  (Li 78)  have  gone  further  and  applied  it  to  the  case  of  a 
Yukawa  potential,  obtaining  encouraging  results. 

Formally  the  Baternan  method  is  equivalent  to  using  the  Schwinger 
variational  principle  [ (Li 78) ,  (Ad75)].    We  shall  thus  present  the 
method  in  this  framework.    Let  us  begin  then  by  considering  the  general 
identity  (Le69) 

[A(z)3  =  At(z)  +  At(z*)+  -  At(z*)"  A^z}"1  At(z),  (2.58) 

where  A  is  an  operator  that  is  L-conjugation  invariant  (A(z)  =  A(z*)+) 
and  At  is  its  trial  value.    Similar  operator  identities  have  been  used 
as  a  basis  for  studying  upper  and  lower  bounds  in  the  theory  of  bound 
state  calculations  (L865).    Taking  A(z)  to  be  equal  to  t(z)  and  using 
equation  (2.3),  which  relates  the  t-operator  to  the  wave  operator  W(z) 
in  the  variational  expression  (2.58),  one  obtains 

[t(z)]  =  Wt(z)  +  Wt(z*)+V  -  Wt(z*)+    V  -  VGQ(z)V  Wt(z),  (2.59) 

where  the  trial  operator  Wt(z)  can  in  general  be  chosen  as 

Wt(z)  =  |f><C|  (2.60) 

Here  |f>  represents  a  set  of  n  basis  functions  fffj^^S' ' ' '  '"V ' 
whose  expansion  coefficients  {Cj.Cg,. . .  »Cn}  are  dtetermined  variational  ly 
through  expression  (2.59).    In  the  usual  variational  procedures  (S172), 
(Ad74),  (Pa74),  (Ku78),  the  set  of  functions  {f}  is  usually  chosen  to 
be  square  integrable;  here  however,  {f]  will  be  chosen  to  be  a  collec- 
tion  of  plane  v/ave  states,  i.e.,  the  set  denoted  by  {k}  or  {ki  .kg*. . .  »Kn}. 
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Thus  letting  j  f  >  be  equal  to  !  k  >  and  maki ng  use  of  equations  (2.59) 
and  (2.60),  one  can  obtain  by  taking  the  proper  variations  (Go50)  the 
expression 

t(q',  q;  z)  *  I+(q')C(q,  z)  ,  (2.51) 
where  £  satisfies  the  relation 

J(z)C(q,  z)  =  1(5)  (2.62) 
and  where  the  components  of  the  matrices  J_  and     are  defined  as 

In-(q)  ■  <£iMq>  (2.63) 

and 

•     J^z)  =  <^|V  -  VG0(z)Y|t,>  .  (2.64) 

Another  equivalent  expression  to  (2.61)  would  be 

t(q',  q;  z)  =  I_+(q*'  )J(z  )'hft)    ,  (2.65) 

but  we  have  chosen  to  use  (2.61)  and  (2.62),  in  order  to  solve  a 
simultaneous  set  of  linear  equations  instead  of  the  matrix  inversion 
implied  by  this  last  expression. 

The  integrals  needed  in  this  procedure,  as  shown  by  equations 
(2.63)  and  (2.64),  are  those  required  for  the  second  Born  approximation 
[(Da51),  (Le56),  (Ho68)].    It  is  readily  found  however  that  integrals 
of  the  type  <£,-  |VG0(z)V|tj>  are  not  trivial  to  evaluate,  and  in  general 
must  be  done  numerically  (Ho68).    For  the  particular  potentials  we  shall 
be  interested  in,  the  Morse  and  exponential  potentials,  one  can  obtain 
the  desired  integrals  from  various  references  [(Da51),  (Le56),  (Li 78 ) ] . 
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In  conclusion  we  should  point  out  the  weaknesses  of  the  Bateman 
method.    From  Eq.  (2.65)  it  is  apparent  that  we  have  obtained  a 
separable  approximation  to  the  t-matrix.    Furthermore,  using  Eq.  (2.1) 
it  is  easy  to  show  that  this  could  be  viewed  as  a  direct  consequence 
of  the  separable  approximation 

V  =  V|J><ljVj|>  1  <£|V  (2.66) 

to  the  potential  [(Ad74),  (Ad75)].    This  is  an  important  point,  since 
Osborn  [(0s73a),  (0s73b)]  has  shown  mathematically  the  non-ccnvergence 
of  t-matrix  elements  based  on  separable  approximations  to  local  poten- 
tials.   This  problem  exists  irrespective  of  the  approach  used  in  con- 
structing the  given  separable  approximation.    Fortunately,  the  diffi- 
culties arise  for  large  values  of  the  momenta  q'  and  q  in  Eq.  (2.63), 
values  at  which  the  t-matrix  elements  are  small  and  have  no  significant 
effect  in  three-body  calculations  (S173).    Finally,  there  remains  the 
problem  of  the  inapplicability  of  the  Bateman  approach  to  singular 
potentials,  a  question  open  to  further  investigation. 


CHAPTER  III 
NUMERICAL  RESULTS  FOR  THE  TWO-BODY  t-MATRIX 

Having  developed  the  formalism  of  the  Variable-Phase  and  -Amplitude 
method  in  the  previous  chapter,  we  will  now  apply  it  to  various  model 
potentials  of  physical  interest.    This  will  serve  to  illustrate  various 
known  analytical  properties  of  the  two-body  t-matrix  and  to  gauge  the 
reliability  of  the  computational  procedure. 

We  will  concentrate  on  calculations  for  the  lowest  1Zg  (attractive) 
and  3EU (repulsive)  potentials  of  H2  because  of  the  role  these  play  in 
the  collision  dynamics  of  H  +  H2>  and  because  their  appearance  in 
dynamical  studies  can  be"  justified  within  a  diatomics-in-molecules  treat- 
ment (Mi 77b).    A  comparison  of  some  of  our  results  will  be  given  with 
others  obtained  within  a  variational  procedure,  developed  to  solve  the 
three-atom  problem  at  low  (thermal)  energies  (Ku78). 

A  study  will  also  be  made  on  the  reliability  of  the  Bateman  method 
in  the  computation  of  t-matrix  elements  for  potentials  of  interest  in 
atomic  and  molecular  scattering  theory. 

1.    Properties  of  the  Two-Body  t-Matrix 
A  brief  summary  of  the  various  properties  of  the  off-shell  t-matrix 
elements,  t-|(q',  q;  E),  will  now  be  given.    This  will  prove  valuable  in 
determining  the  reliability  of  any  given  computational  procedure  used  to 
obtain  t-| . 

The  most  important  properties  satisfied  by  t-j(q',  q;  E)  are  those 
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of  symmetry, 

t^q',  q;  E)  =  t-j (q ,  q!;  E)  ,  (3.1) 

and  of  off-shell  unitarity, 

Imft^q',  q;  E)]  =  -  mk£t*  (k£J  q';  E)t](kE,  q;  E)  .  (3.2) 

The  symmetry  relation  is  valid  at  all  energies,  but  the  unitarity  rela- 
tion is  valid  only  for  positive  energies.    These  relations  are  well 
established  in  the  literature  and  will  not  be  derived  (Wa67).  It 
suffices  to  say  that  in  this  case  they  are  a  direct  consequence  of 
assuming  a  real  local  potential. 

If  one  considers  the  differential  equation  satisfied  bycj^q,  k£;  r), 
Eq.  (2.10),  it  is  apparent  that  oo-j  reduces  to  the  on-shell  wave  function 
when  q  =  k£.    Similarly,  the  off-shell  t-matrix  elements  are  defined 
such  that  there  is  a  smooth  transition  to  the  on-shell  behavior  as  the 
momenta  involved  approach  the  on-shell  quantities.    This  is  exemplified 
by  the  relations  (Wa67) 

"(+) 

w,(q,  kE;  r)  j-i(qr)  -  Trmqt1(kE,  q;  E)h,    (k£r)  (3.3) 

1  r  ->  oo 

and 

ImtMq,  kF;  E)] 

tan[61(kE)]  =   1~   (3  4l 

1    E  Rett^q,  kE;  E)]  [6^> 

which  correspond  to  well  known  on-shell  relations  for  q  =  k£. 

Another  important  aspect  of  the  t-matrix  elements  is  their  behavior 
as  a  function  of  momenta  and  energy.    In  this  respect,  there  are  various 
limits  of  interest.    If  the  energy  E  and  initial  momentum  q  are  fixed, 
then 

lim    t^q',  q;  E)  =  0  (3.5) 
q  1  ->  oo 


and  therefore  because  of  symmetry 

lim     t^q'  ,  q;  E)  -  0  .  (3.6) 
q 1  oo 

These  relations  are  apparent  if  one  considers  the  oscillatory  nature  of 
^(q,  k£;  r)  and  j^q'r)  in  expression  (2.13)  as  q  and  q'  approach 
infinity.    A  clearer  understanding  of  relations  (3.5)  and  (3.6)  could 
be  obtained  from  a  semiclassical  analysis  of  the  two-body  t-matrix  as 
a  function  of  energy  and  momenta.    Such  an  analysis  has  been  carried 
out  by  Korsch  and  Mohlenkamp  [(Ko76),  (Ko77a)]  for  the  case  of  purely 
repulsive  potentials.    These  authors  have  classified  the  various  regions 
in  the  q-q'  momentum  plane  in  which  t-|  corresponds  to  classically  allowed 
or  forbidden  processes.    Classically  forbidden  processes  would,  for 
example,  be  associated  with  complex-valued  classical  paths  or  trajec- 
tories. 

As  a  function  of  energy,  the  t-matrix  must  also  show  continuity  as 
the  energy  goes  from  positive  to  negative,  and  consequently 

lim     Mq',  q;  E)  =    lim     Re[t,(q\  q;  E)]  (3.7) 
E  +  0-  E  ■*  0+ 

and 

lim       ImiMq',  q;  E)]  =  0    .  (3.8) 
E  -  0+ 

For  negative  energies,  the  t-matrix  elements  are  real  and  have  poles  at 
energies  corresponding  to  the  bound  states  of  the  system.    This  can  be 
established  by  considering  the  expansion  of  t(q',  q;  z)  in  terms  of  the 
eigenf unctions  of  the  full  Hamiltonian  of  the  system  (Ne66).    The  result 
is  as  follows: 
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t(5'.  3;  z)  .  i  (en1  -  V)(z-eq'  4     (q1 )  *nlm  (?) 

nlm         U  -  en1  J 
+  t(q\  q;  eq+) 

+  j  dq"[(£qll  -  sq  -  ie)"1  +  (z  -  e^"1] 

x  t(q',  q";  eq..+)  t*(q",  q;  eq»)    ,  (3.9) 

where  ^nim(q)  is  a  bound  state  function  in  the  momentum  representation 
having  the  eigenvalue  of  en-j.    The  quantities,  eq,  are  eigenvalues  to 
the  scattering  states  of  the  system. 

2.  Numerical  Calculations  Using  the  VPA  Equations 

We  will  now  present  various  results  for  the  H-  system  for  the 
1   +  ? 

lowest   Eg    and    Zu  potentials  [(Mi77c),  (0178)].    Some  of  these  results 
have  already  been  reported  in  the  literature  [(Be78),  (Ku78),  (Ko77a)]. 
A  Morse  potential 

V(r)  =  D[l  -  e"a(r-ro)]2  -  D  (3.10) 

was  used  to  represent  the  *Eg+  interaction  potential,  and  a  Hulthe"n 
potential 

V(P)  -  Ae"r/a[l  -  e-^0]"1  (3.11) 

3 

was  used  for  the  repulsive    E    interaction  potentials.    The  parameters 

used  were  D  =  4.786  eV,  a  =  2.1123  A"1,  rQ  =  0.7411  A,  A  =  20.11  eV  and 

o 

a  =  0.4984  A  (M177c).    This  particular  study  will  prove  interesting, 
since  it  serves  to  draw  a  contrast  between  the  amaTytic  properties  of 
t-|  and  w-j  obtained  from  an  attractive  potential,  and  those  obtained  from 
a  strictly  repulsive  potential. 


In  figures  (3- la)  through  (3-lc),  we  have  plotted  various  components 
of  the  positive  energy  off-shell  wave  function  w    for  the  case  of  a  Morse 
potential.    If  one  recalls  the  VPA  expression  for  w.  given  by  Eq.  (2.26), 
it  is  clear  that  one  can  consider  oo-j  (q ,  k^;  r)  to  have  two  components: 
j-j  (qr)  and  y^{q,  k^;  r),  where  x-i  is  defined  by 

X-,{q,  k£;  r)  =  ^(q,  k£;  r)  -  j^qr)  ,  (3.12) 

and  is  often  called  the  wave  defect  (Sr68).    The  term  j. (qr)  in  Eq. 
(3.12)  is  present  in  equations  (2.26)  and  (3.12)  because  of  the  inhomo- 
geneous  character  of  the  differential  equation  obeyed  by  w-j  [see  Eq. 
(2.10)].    Thus,  we  see  from  Eq.  (3.12)  that  for  positive  energies 

Re^l  =  j^qr)  +  Retx-,1  (3.13) 

and 

Im[u}]  =  Im[Xl]  .  (3.14) 

Equation  (3.13)  implies  that  the  real  part  of  w,  is  highly  influenced 
by  the  driving  term  j-j(qr).    This  can  be  seen  in  figures  (3-la)  and 
(3-lc),  which  show  how  the  function  [2/u]1/2  Re[w0(K,  KE;  R)/K]  oscil- 
lates about  the  driving  term  [2/tt]  1^2[ j'q(KR)/K]  .    The  data  used  to 
obtain  these  figures  were  !<=  0.23  a.u.  and  E  =  0.01  a.u.    If  ona.  were 
to  increase  the  value  of  either  K  or  E,  the  off-shell  wave  function 
would  become  much  more  oscillatory.    Similar  results  were  also  obtained 
for  the  Hulthe'n  potential.    We  note,  however,  that  an  attractive  poten- 
tial pulls  the  nodes  of  the  off-shell  wave  function  in,  whereas  a  repul- 
sive potential  pushes  them  out.    For  negative  energies,  the  off-shell 
wave  function  becomes  real.    An  illustration  of  to,  at  negative  energies 
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is  given  by  figures  (3-2),  (3-3)  and  (3-4)  for  the  Morse  potential 

case.    Since  the  t-matrix  elements  have  poles  at  the  bound  states  of 

the  system  [see  Fig.  (3-5)],  we  see  from  the  relation  connecting  t-| 

and  a),  [Eq.  (2.13)]  that  the  amplitude  of  to,  must  grow  as  the  energy 

1/2 

approaches  a  bound  state  eigenvalue.    In  Fig.  (3-2),  [2/-]  [cog/K] 
is  plotted  versus  R  for  K  =  4.475  a.u.  and  an  energy  value  of  E  = 
-0.1446  a.u.    The  large  magnitude  of  the  amplitude  displayed  by  this 
function,  when  compared  to  the  result  at  E  =  -0.153  a.u.  given  in  Fig. 
(3-3),  is  a  clear  indication  that  the  energy  is  quite  near  to  an  eigen- 
value of  the  potential.    A  similar  effect  can  be  seen  at  higher  partial 

1/2 

waves,  as  is  shown  in  Fig.  (3-4)  where  [2/ir]      fu^(K,  K£;  R)/K]  is 
plotted  versus  R  for  K  =  4.5  a.u.  and  E  =  -0.153  a.u.    Note  that  the 
effect  of  the  driving  term,  j^(KR),  is  always  present  in  all  of  the 
figures  (3.2)  through  (3.4).    Within  the  interaction  region  of  the 
potential,  it  is  the  wave  defect  x-|  which  dominates;  whereas  the  driving 
term  dominates  in  the  asymptotic  region.    For  the  Hulthe"n  potential,  the 
off-shell  wave  function  does  not  display  the  large  changes  in  magnitude 
mentioned  above  for  the  Morse  potential,  since  the  Hulthe"n  potential 
used  in  this  study  does  not  support  bound  states. 

Various  studies  of  the  partial  wave  t-matrix  elements  as  a  function 
of  energy  and  momenta  were  carried  out  for  both  the  Morse  and  Hulthe'n 
potentials.    We  have  already  shown  the  effect  of  the  bound  states  on 
the  off-shell  wave  function  co-j .    Their  effect  on  the  t-matrix  elements 
is  also  quite  important.    In  Table  (3-1),  we  have  tabulated  the  Morse 
tg(q,  q;  E)  matrix  elements  for  q  =  4.0  a.u.  as  a  function  of  energy. 
Column  I  gives  the  VPA  results  and  column  II  the  corresponding  results 
obtained  by  a  variational  method  developed  by  Kuraoglu  and  Micha  (Ku78). 
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Table  (3-1)    The  t«(q,  q;  E)  matrix  elements  are 
given  for  q  =  4.0  a.u.    Column  I 
gives  the  VPA  calculations  and  column 
II  the  results  obtained  from  a  varia- 
tional procedure  (Ku78).  (Morse 
Potential ) 


Energy  (a.u.)  tn(q,  q;  E)  x  10 


I 

II 

-0.18 

-2.5971 

-2.5798 

-0.17 

-4.4003 

-4.3841 

-0.16 

-2.0273 

-2.0069 

-0.15 

-8.8552 

-8.7964 

-0.14 

9.9843 

10.0183 

-0.13 

2.0411 

2.0680 

-0.12 

1.5197 

1.5486 

-0.11 

-1.2378 

-1.2106 

-0.10 

1.0856 

1.1159 
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We  note  that  in  general  the  agreement  is  quite  good.    A  more  detailed 
picture  about  the  structure  of  the  Morse  tQ  matrix  elements  at  negative 
energies  can  be  seen  from  the  results  given  in  Table  (3-2)  and  Fig.  (3-5), 
each  of  which  describes  the  energy  range  corresponding  to  the  first  four 
poles.    The  figure  shows  that  the  t-matrix  elements  change  sign  as  one 
goes  across  a  pole,  in  agreement  with  equation  (3.9).    For  the  Hulthe'n 
potential,  the  t-matrix  elements  show  a  very  smooth  and  plain  structure 
as  can  be  seen  in  Table  (3-3),  where  tQ(q,  q;  E)  is  given  versus  energy 
for  q  =  4.5  a.u.    The  agreement  with  the  variational  results  is  again 
quite  good.    At  positive  energies,  we  have  shown  in  Fig.  (3-6)  the 
behavior  of  the  Morse  tQ(q,  q;  E)  matrix  elements  near  E  =  0  a.u. 
[q  =  4.0  a.u.].    We  see  the  continuity  one  expects;  i.e.,  both  equation 

(3.7)  and  (3.8)  are  satisfied.    In  general,  the  t-matrix  elements  of 
both  the  Morse  and  Hulthe'n  potentials  show  a  much  more  oscillatory 
structure  at  positive  energies.    Figures  (3.7a)  and  (3.7b)  illustrate 
this  for  the  diagonal  (K1  =  K)  Hulthe'n  tQ(K' ,  K;  E)  matrix  elements. 
The  large  oscillations  seen  in  these  figures  do  not  extend  indefinitely. 
Korsch  and  Mohlenkamp  (Ko77a)  have  extended  the  present  calculations  to 
higher  energies  and  have  shown  that  the  amplitude  of  the  t-matrix  even- 
tually decays  to  zero.    In  the  comparison  of  their  semiclassical  results 
with  the  present  quantum  mechanical  ones,  they  obtained  good  agreement. 
As  a  function  of  momenta,  the  t-matrix  elements  show  a  much  greater 
oscillatory  nature.    This  can  be  seen  in  Fig.  (3-8)  where  the  Morse 
tg(K,  K;  E)  matrix  elements  are  plotted  versus  K  for  E  =  -0.1378  a.u., 
and  in  Fig.  (3-9)  where  the  Morse  tQ(K';  K;  E)  matrix  elements  are 
plotted  versus  K'  for  E  =  -0.153  a.u.  and  K  =  4.5  a.u.    Both  figures 

(3.8)  and  (3.9)  show  how  the  t-matrix  goes  to  zero  asymptotically  in 


Table  (3-2)    Selected  values  of  tQ(q,  q;  E)  for 
q  =  4.475  a.u.  are  presented  about 
the  first  four  poles.  (Morse 
Potential ) 


Energy  (a.u.) 


t(K,K:E)       Location  of  Pole 


■0.165078 
■0.165075 
•0.165074 
■0.165072 


-2.6856 
-5.4287 
174.7388 
5.0667 


-0.1650741 


-0.144800 
-0.144600 
•0.144400 
■0.144200 


•1.2866 
■3.5240 
4.6930 
1.4018 


-0.144486 


■0.125600 
■0.125400 
■0.125200 
■0.125000 


•0.0290 
■0.0936 
0.2220 
0.0656 


-0.125268 


•0.107800 
-0.107600 
■0.107400 
■0.107200 


■0.0011 
-0.0026 
0.0243 
0.0020 


-0.1137421 


Table  (3-3)    Selected  values  of  tQ(q,  q;  E) 
for  q  =  4.0  a.u.  are  given. 
Column  I  gives  the  VPA  calcula- 
tions and  column  II  the  results 
obtained  from  a  variational  pro- 
cedure (Ku78).    (Hulth£n  Poten- 
tial) 


Energy  (a.u.) 

tQ(q,  q; 

E)  x  103 

I 

II 

-0.18 

6.4807 

6.4592 

-0.17 

6.2683 

6.2477 

-0.16 

6.0525 

6.0300 

-0.15 

5.8277 

5.8057 

-0.14 

5.5946 

5.5742 

-0.13 

5.3574 

5.3348 

-0.12 

5.1079 

5.0868 

-0.11 

4.8526 

4.8295 

Retn(K,K;E)X104(a.u.)    Im  tn(K;K;E)X104  (a.u.) 


o 


■*n — r 


I  1  i  1 — r~ ■ — r 


0 

o 


-r— | — 7-1 — i — [— i — r 


10 

i 


j_i  i  ,  '  i  ' 


 L_ 


o 

if)  ; 

c_  (0 

O  w 

2  ^ 


n: 


'    i    I    ■  L 


ix) 


(■n-B)  JDIXQ-^)0} 


m  ■ 
i 


o 

(/) 

c 

a> 

E 

oJ 

X 

•>— 

t- 

+J 

10 

E 

i — 

<T3 

•r— 

4-> 

u 

C 

CD 

+J 

Li 

O 

a. 

c 

c:  -c 

o 

+J 

Dl 

(a 

•i — 

Q 

<T3 

CO 


QJ 
=S 

cn 


□ 
>» 

01 

to 
c 


o 


o 

to 
+-> 

c 

0) 

E 

QJ 


B  r- 

ta 

c 
«  a) 

-  +j 
^  o 

—  a. 

e—  C 
(T3  VI) 
C  -C 
O  4-> 

enr- 
ol 3 
•r-  rc 

Q  


X5 
I 


O) 


</>  v- 
+J  +J 
c  c 

CU  CD 

a  o 
,—  o_ 

0) 

T-  S~. 

5-  O 
+->  2: 

(0  

E 


UJ 


(O  

c  o 

o  -t-> 

01 

(0  s- 

■r-  O 
Q  <4- 


co 
I 

CO 


s- 

01 


81 


Figure  (3-9)    t0(K',  K;  E)  for  E  =  -0.153  a.u. 

and  K  =  4.5  a.u.    (Morse  Potential) 
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.  84 

K  and  K'  [see  equations  (3.5)  and  (3.6)].    Comparable  behavior  was 
observed  at  positive  energies  as  well,  irrespective  of  the  attractive 
or  repulsive  nature  of  the  potential.    Table  (3-4)  shows  a  momentum 
study  for  the  Hulthe*n  potential  t0(q',  q;  E)  matrix  for  E  =  -0.15  a.u. 
and  the  1  values  0,  5  and  10.    These  results  indicate  that  the  diagonal 
q'  =  q  t-matrix  elements  are  consistently  larger  in  magnitude  than  the 
off-diagonal  q:  f  q  t-matrix  elements.    Column  I  in  this  table  gives 
the  VPA  results  and  Column  II  the  corresponding  variational  results 
(Ku78).    As  a  function  of  the  partial  wave  parameter  1,  the  t-matrix 
generally  shows  an  oscillatory  structure  for  low  partial  waves  and  a 
smoothly  decaying  structure  for  high  partial  waves.    Figures  (3-10a) 
and  (3-10b),  nicely  illustrate  this  effect  for  both  the  real  and  imagi- 
nary components  of  the  Kulthdn  ti(K* ,  K;  E)  matrix  elements  (K1  =  K.  =  4.0 
a.u.  and  E  =  0.01  a.u.).    Using  the  same  data,  *e  see  in  figures  (3-lla) 

and  (3- lib)  plots  of  the  real  and  imaginary  components  of  the  total 

->•->■  ■*■  ** 

t-matrix  [t(K\  K;  E)]  versus  momentum  transfer  £P  =  |K'  -  K| .  These 

plots  imply  that  the  t-matrix  elements  become  siraTler  as  the  momentum 
transfer  increases.    In  this  particular  case,  thiirty-five  partial  waves 
were  required  for  the  total  t-matrix  to  converge-   As  the  energy  and 
momenta  arguments  become  larger,  this  number  willi  generally  increase. 
Partial  wave  convergence  will  be  discussed  further  in  Chapter  VII, 
where  the  two-body  t-matrix  elements  are  used  in  the  three-body  atom- 
diatom  scattering  problem. 

The  t-matrix  results  we  have  reported  here  ware  all  obtained  from 
Eq.  (2.13)  using  a  numerical  quadrature,  which  in  this  case  was  a  com- 
bination Simpson's  and  Newton's  3/8  rule  (Hi 56 ) .  In  order  to  test  the 
reliability  of  this  procedure,  we  have  used  Eq.  (3B.2)  to  obtain  phase 
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Table  (3-4)    Selected  values  of  t,(q',  q;  E)  are 
given  for  E  =  -0.15  a.u.  and  the  par- 
tial wave  numbers  1  =  0,  5  and  10. 
Column  I  gives  the  VPA  calculations 
ant  column  II  the  results  obtained  from 
a  variational  procedure  (Ku78).  (Hulthgn 
Potential) 

1  q'         q  t^q',  q;  E)  x  103 

 I  II 

0       3.0       1.0  0.7899  0.7894 

2.0  4.0665  4.0693 

3.0  10.187  10.179 

4.0  2.1074  2.1160 

5.0  0.2154  0.2476 


5 

3.0 

1.0 

0.0207 

0.0207 

2.0 

1.1745 

1.1755 

3.0 

5.3773 

5.3709 

4.0 

1.0412 

1.0434 

5.0 

0.0513 

0.0547 

10 

3.0 

1.0 

0.0005 

0.0009 

2.0 

0.1174 

0.1186 

3.0 

1.4905 

1.4857 

4.0 

0.2644 

0.2668 

5.0 

0.0093 

0.0125 
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shifts  from  the  half-on-sheVI  t-matrix  elements.    Typical  results  are 
shown  in  Table  (3-5),  where  Morse  potential  phase  shifts  are  tabulated 
for  the  energies  E  =  0.01  and  E  =  0.04  a.u.    Columns  I  and  III  correspond 
to  phase  shifts  obtained  using  Eq.  (3.2),  whereas  columns  II  and  IV  list 
analogous  results  obtained  from  the  numerical  integration  of  the  on-shell 
Schrodinger  equation.    A  De  Vogelaere  integration  scheme  (Le63)  was  used 
for  the  on-shell  calculations  using  the  procedure  developed  by  Bernstein 
(Be60).    Comparably  good  results  were  obtained  for  various  other  poten- 
tial data.    Other  factors  that  demonstrate  accuracy  of  the  present  method 
are  symmetry  and  unitarity  [see  equations  (3.1)  and  (3.2)].  Generally, 
agreement  to  three  or  four  decimal  places  was  obtained.    An  exception 
to  this  occurred  at  high  values  of  momenta  where  the  t-matrix  elements 
are  quite  small  in  magnitude.    This  indicates  that  the  mesh  size,  over 
which  the  off-shell  wave  function  is  determined  and  over  which  the  quad- 
rature in  Eq.  (2.13)  is  carried  out,  may  be  too  large.    Given  the  extreme 
oscillatory  nature  of  the  functions  involved,  a  smaller  step  size  would 
be  desirable.    However,  a  compromise  on  the  step  size  is  necessary  to 
prevent  the  accumulation  of  error  from  being  too  large. 

It  is  in  this  respect  that  the  comparison  potential  method  offers 
an  advantage,  since  the  quadratures  involved  in  that  procedure  are  done 
analytically.    We  found  though,  that  the  VPA  procedure  showed  more 
stability  in  the  propagation  of  the  off-shell  wave  function  through  non- 
classical  regions.    This,  however,  warrants  further  investigation,  since 
it  is  a  function  of  the  efficiency  and  sophistication  of  the  computer 
codes  used. 


Table  (3-5)  Columns  I  and  III  are  phase  shift  values 
obtained  from  half-on-shel 1  quantitites, 
and  columns  II  and  IV  are  the  correspond- 


ing 
are 

on-shell  quantities.    All  phase 
Mod  (it).    (Morse  Potential) 

shifts 

o-|(E 

=  0.01  a.u. ) 

5-,(E  =  0.04 

a.u. ) 

1 

I 

II 

III 

IV 

0 

0.8655 

0.8647 

-1.2732 

-1.2757 

1 

-0.8463 

-0.8471 

0.1843 

0.1819 

2 

0.4431 

0.4421 

1.5293 

1.5270 

3 

-1.5482 

-1.5494 

-0.3781 

-0.3805 

4 

-0.5379 

-0.5373 

0.7473 

0.7446 

5 

0.3403 

0.3390 

-1.3752 

-1.3780 

6 

1.0825 

1.0815 

-0.4598 

-0.4627 

7 

-1.4480 

-1.4494 

0.3532 

0.3524 

8 

-0.9677 

-0.9685 

1.0698 

1.0702 

9 

-0.6154 

-0.6157 

-1.4366 

-1.4477 

10 

-0.3781 

-0.3896 

-0.8966 

-0.9154 

96 


3.    Numerical  Results  Obtained  Using  the  Batsman  Method 
From  the  results  given  in  the  last  section,  it  should  be  clear 
that  a  method  which  gives  the  total  t-matrix  t(q',  q;  E)  without  a 
partial  wave  expansion  would  be  highly  desirable.    Such  a  method  would 
hopefully  avoid  the  problems  of  computing  a  highly  oscillatory  function 
for  many  partial  waves.    To  this  end,  we  now  present  some  results 
obtained  using  the  Bateman  method  as  outlined  in  Section  (2-6). 

Computationally  all  that  is  needed  to  obtain  t(q\  q;  E)  using  the 
Bateman  method  is  the  evaluation  of  the  integrals  of  the  type  <k1-|V|kJ-> 
and  <k1-|VGQV|kj>  along  with  a  procedure  for  solving  a  set  of  complex 
simultaneous  equations  [see  Eq.  (2.61)].    For  the  Yukawa  and  exponential 
potentials,  analytic  expressions  for  the  above  integrals  may  be  found  in 
the  references  given  in  Section  (2-6)  [(Da51),  (Le56),  (L177)].  The 
integrals  for  any  other  potential  which  is  a  linear  combination  of 
exponentials  may  be  obtained  using  the  relation  (Br70) 

+00 

V(r)  =  /  dX     p(X)  e~Ar    ,  (3.15) 

-00 

where  p(X)  is  some  linear  combination  of  delta  functions.  For  the  Morse 
potential ,  one  has 

p(X)  =  D[e2ar°  6(X  -  2a)  -  2fi*r°  5(X  -  a)]  .  (3.16) 

The  set  of  complex  linear  equations  given  by  Eq.  (2.61)  was  solved 

using  the  IMSL  library  program  LEQTIC  ( Ai 75 ) .    Since  the  input  matrix 
->• 

elements  <k^ | VGQV ] kj>  have  quite  complicated  analytic  expressions,  their 
accuracy  was  verified  both  analytically  and  numerically.    In  all  cases 
the  analytic  and  numerical  results  agreed  to  tern  or  more  decimal  digits. 
In  the  present  calculations,  we  have  found  tthat  the  Bateman  method 


worked  quite  well  for  soft  potentials  such  as  the  Yukawa  potential 
studied  by  Walters  (Wa71),  and  Rosenthal  and  Kouri  (Ro73)  [V(r)  = 
-1.1825  r"1  e"r  in  units  such  that  h  =  m  =  aQ  (unit  of  length)  =  1], 
This  can  be  seen  in  Table  (3-6),  where  we  have  tabulated  the  on-shell 
t-matrix  t (q ' ,  q;  E)  for  q'  =  q  =  kE  =  1.816  a^1  as  a  function  of  the 
angle  Q^>t  between  q1  and  q.    Comparing  the  Bateman  results  to  the  VPA 
results,  we  see  that  the  Bateman  results  are  much  better  for  small 
values  of  the  momentum  transfer  |q'  -  q| .    This  was  found  to  be  true 
in  all  the  studies  that  were  made.    The  set  of  plane  wave  states  {k} 
that  was  used  in  the  above  calculations  is  given  in  Table  (3-7).  Because 
of  the  rotational  invariance  of  the  potentials  involved,  it  is  only 
necessary  to  choose  the  set  of  vectors  {kj,  k£,         kn)  to  lie  on  the 
(q'  -  q)  plane.    If  one  specifies  q  as  the  reference  vector,  it  is  then 
only  necessary  to  specify         and  kp.    The  results  given  in  Table  (3-5) 
were  quite  insensitive  to  any  variation  in  the  choice  of  {k}.  Better 
agreement  was  obtained,  however,  if  the  set  included  both  q'  and  q  as 
members  of  the  set.    A  comparison  with  the  calculations  given  by  Walters 
(Wa71),  and  Rosenthal  and  Kouri  (Ro73)  is  also  given.    The  agreement 
between  these  results  and  those  obtained  from  the  VPA  method  is  better 
than  that  obtained  from  the  Bateman  method.    The  Bateman  method,  however, 
is  much  easier  to  apply  in  the  present  case  than  either  of  these  other 
methods.    Another  observation  that  can  be  made  from  Table  (3-6)  is  that 
the  Bateman  method  yields  better  agreement  for  the  imaginary  component 
of  t  than  for  the  real  component.    The  reason  for  this  difference  is 
not  obvious  and  remains  an  open  question. 

We  did  extend  the  present  study  to  soft  exponentials  with  parameters 
and  units  comparable  to  those  shown  above  for  the  Yukawa  case.  The 
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Table  (3-7)    Basis  {k}  used  in  the  Bateman  calcu- 
lations^!' ven  in  Table  (3-6).  The 
vector  k\l)  was  chosen  to  lie  along 
the  direction  of  the  initial  momentum 
vector  3  =  tr,  all  other  lie  on  the 
(q'-q)  momentum  pjane.       is  the  angle 
between  £(n)  and  k(l). 


n 

6[k(n),k(l)]  (radius) 

|k(n)| 

1 

0 

1.816 

2 

0.785 

0.605 

3 

1.571 

1.816 

4 

2.356 

0.605 

5 

3.1416 

1.816 

6 

0.1963 

1.816 

7 

0.5890 

0.454 

8 

1.3744 

1.816 

9 

2.1598 

0.454 

10 

2.945 

1.816 

11 

0.295 

1.816 

12 

0.393 

2.421 

13 

1.178 

1.816 

14 

1.963 

2.421 

15 

2.749 

1.816 

100 


results  were  generally  good,  provided  one  again  concentrates  on  small 
momentum  transfers.    We  do  not  report  any  off -shell  studies  here, 
although,  a  similar  quality  in  the  results  was  obtained.    When  we 
extended  the  Bateman  method  to  potentials  applicable  to  chemical 
problems,  e.g.  exponential  potentials  that  are  very  steep  and  localized 
in  coordinate  space,  the  results  obtained  were  totally  inadequate.  The 
sensitivity  to  the  choice  of  it]  became  very  large.    In  light  of  these 
findings,  no  further  investigations  were  carried  out.    A  possible 
explanation  of  this  negative  outcome  might  lie  in  the  fact  that  for  a 
hard  potential,  the  interaction  region  is  highly  localized  leading  to 
a  corresponding  delocalized  function  in  momentum  space.    Thus,  for  such 
potentials,  one  would  expect  the  number  of  basis  function  in  the  set 
it}  to  be  very  high. 


CHAPTER  IV 
THE  THREE-BODY  PROBLEM 
In  Chapter  I  we  discussed  the  feasibility  of  treating  an  atom- 
diatom  scattering  process  as  an  effective  three-body  problem  in  such 
a  way  that  one  could  make  use  of  two-body  interaction  data.    It  was 
pointed  out  there    that  one  could  construct  a  proper  representation  of 
the  interaction  potential  for  a  reactive  triatomic  system  in  terms  of 
two-body  quantities,  only  if  one  introduced  a  set  of  spin-dependent 
atomic  pair  interactions.    Since  our  present  interest  lies  in  high 
energy  atom-diatom  collisions,  we  proposed  the  use  of  a  simplified 
impulsive  model  in  which  the  total  interaction  potential  would  be 
assumed  a  sum  of  spin  independent  two-body  potentials.    Having  already 
addressed  in  Chapter  I  the  merits  and  limitations  to  such  a  representa- 
tion of  the  potential,  we  turn  now  to  the  solution  of  the  scattering 
problem. 

The  formalism  we  shall  adopt  in  treating  the  scattering  of  an  atom- 
diatom  system  will  be  based  on  that  given  by  Faddeev  (Fa61),  who 
developed  the  first  rigorous  mathematical  treatment  of  the  three-body 
problem.    We  shall  not  attempt  to  address  the  more  mathematical  aspects 
of  the  problem.    Instead,  we  shall  restrict  ourselves  to  a  discussion 
of  difficulties  encountered  when  one  tries  to  apply  the  Lippmann-Schwinger 
equation  to  the  three-body  problem  and  to  a  formal  derivation  of  the 
various  transition  operators  of  physical  interest.    The  resulting 
expressions  for  these  transition  operators    will  then  be  used  to  construct 
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a  multiple  collision  expansion,  which,  as  was  pointed  out  in  Chapter  I, 
can  be  used  to  describe  the  scattering  of  an  atom  and  a  diatom  as  a 
sequence  of  atom-atom  encounters  [ (Ch71) , (Mi 75b) ] . 

1.    Problems  with  the  Lippmann-Schwinger  Equation 
The  main  difficulty  with  the  three-body  problem  lies  in  the 
existence  of  various  arrangement  channels.    By  this  we  mean,  that  given 
three  particles  (atoms)  A,  B  and  C,  one  can  have  any  of  four  possible 
arrangement  channels,  i.e., 

channel  (1)  corresponding  to  A  +  (B,  C), 
channel  (2)  corresponding  to  B  +  (A,  C), 
channel  (3)  corresponding  to  C  +  (A,  B), 
and  channel  (0)  corresponding  to  A  +  B  +  C, 

where  an  interacting  pair  is  denoted  by  enclosing  the  given  particles 
in  parentheses.    This  implies  that  for  reactive  scattering  there  are 
difficulties  in  implementing  the  boundary  conditions,  at  least  in  the 
usual  coordinate  space  formulations  of  the  scattering  problem  [(Ma66), 
(Li71),  (Mi 72a ) ] ,  and  for  dissociative  scattering  there  are  even  problems 
with  regard  to  the  proper  formulation  of  the  boundary  conditions  (Me76). 
For  rearrangement  scattering  the  problems  begin  with  the  choice  of 
coordinates,  since  the  coordinates  that  are  best  suited  for  describing 
the  motion  of  the  three-body  system  in  the  entrance  channel  are  not  the 
most  appropriate  ones  for  describing  the  motion  in  the  exit  channel 
(Ma66).    Going  further,  let  us  consider  an  integral  equation  for  the 
three-body  problem  analogous  to  the  Lippmann-Schwinger  equation  for  the 
two-body  transition  operator,  i.e.,  the  equation 

T(Z)  =  V  +  VGQ(Z)T(Z)  (4.1) 


103 


where  T(Z)  is  the  general  three-body  T-operator,  GQ(Z)  is  the  free- 
particle  propagator,  Z  is  an  energy  parameter  (E  +  ie)  and  V  is  the 
potential  for  the  three-body  system.    The  problem  with  the  boundary 
conditions  now  manifests  itself  in  the  form  of  the  "uniqueness"  question 
with  regard  to  the  solutions  to  Eq.  (4.1)  [(Sc74a),  (Re77)].    Since  a 
three-body  system  allows  for  the  existence  of  two-body  bound  states  at 
a  continuum  of  energies  within  the  range  of  physical  interest,  one  must 
provide  additional  constraints  which  select  those  solutions  satisfying 
the  proper  physical  boundary  conditions.    Furthermore,  because  of  the 
presence  of  these  subsystem  bound  states,  one  finds  that  even  after  the 
center-of-mass  motion  has  been  factored  out,  the  kernel  contains 
unfactorable  delta  functions  (Re77).    Each  of  these  delta  functions 
expresses  the  conservation  of  momentum  of  the  spectator  atom  in  the 
presence  of  the  interacting  atom-atom  pair.    A  very  nice  discussion  of 
all  these  problems  is  given  in  the  references  [(Kr71),  (Sc74a),  ( S i 7 1 ) , 
(Re77),  (Wa67)]. 

2.    The  Multichannel  Transition  Uperators 
It  is  clear  from  the  discussion  in  the  previous  section,  that  the 
Lippmann-Schwinger  equation,  Eq.  (4.1),  must  be  reformulated  by  speci- 
fically taking  into  account  the  existence  of  the  various  arrangement 
channels.    This  can  be  done  in  a  straightforward  manner  by  considering 
the  structure  of  the  potential  characterizing  the  three-body  system,  and 
by  introducing  various  operators  describing  the  interaction  of  two 
particles  in  the  presence  of  a  third  [(Fa61),  (Si7/1),  (Sc74a)]. 
Let  us  begin  by  considering  the  total  Hamiltinian, 

H  =  HQ  +  V,  (4.2) 
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where  H„  is  the  free-motion  Hamiltonian  for  the  three  atom  system,  and 
o 

introduce  the  assumption 


3 

V  =    E    V.  ,  (4.3) 
j=l  J 


i.e.,  that  of  pair-wise  additive  potentials.    Each  Vj  is  an  atom-atom 
potential  corresponding  to  the  asymptotic  interaction  in  channel  j,  and 
may  be  specified  using  the  procedure  outlined  in  Chapter  I.  Asymptoti- 
cally in  each  channel  j,  the  system  can  be  described  by  the  channel 
Hamiltonian 

VHo  +  Vj  •  <4-4> 

Therefore,  the  total  Hamiltonian  for  the  system  can  be  written  as 

H  =  H,  +  v(j)  (4.5) 

where 

3 

I 

k=l 


V(J)  -  V  -  Vj  -   Ujk  Vk  (4.6) 


is  the  channel  interaction  potential  and  6^  is  equal  to  (1  -  5-^)- 
Having  made  this  connection  it  is  apparent  that  iky  considering  various 
partitions  to  the  total  Hamiltonian  in  terms  of  channel  quantities, 
such  as  H.  and  V^),  one  should  be  able  to  construct  the  various  opera- 
tors  needed  in  the  reformulation  of  Eq.  (4.1).    Important  in  this  regard 
are  various  two-body  operators,  such  as  the  charms!  resolvents  and  two- 
body  transition  operators,  which  we  will  now  introduce. 

The  resolvent  relations  that  are  needed  follow  directly  from  equa- 
tions (4.2)  through  (4.6),  using  the  identities  (01668) 
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(A  -  B)-l  =  A-l  +  A-l  B(A  -  B)"l  (4.7) 

and 

(A  -  B)"1  =  A"1  +  (A  -  B)"1  BA"1  (4.8) 
The  expressions  obtained  through  Eq.  (4.7)  would  be 

G(Z)  =  Gj(Z)  +  Gj(Z)v(J)G(Z)  (4.9) 

and 

Gj(Z)  =  GQ(Z)  +  G0(Z)VjGj(Z),  (4.10) 

where  G(Z)  =  (Z  -  H)~*  and  Gj(Z)  =  (Z  -  Hj)-1.  These  equations  remain 
valid  for  j  =  0  if  one  defines  Vq  =  0,  so  that  v(°)  =  V.  The  two-body 
transition  operators  can  be  defined  as 

Tj(Z)  =  Vj  +  G0(Z)VjTj(Z),  (4.11) 

analogous  to  the  definition  given  for  t(z)  given  in  chapter  II.  At 
this  point  it  should  be  emphasized  that  T.  and  G.  are  two-body  operators 
embedded  in  a  three-body  space;  i.e.,  they  characterize  the  behavior  of 
two  interacting  particles  in  the  presence  of  a  third  spectator  particle. 
This  should  be  clear  from  the  structure  of  GQ,  which  is  the  free  propa- 
gator for  the  entire  three-body  system. 

Let  us  return  now  to  the  Lippmann-Schwinger  equation,  and  use 
Eq.  (4.2)  to  express  the  three-body  T-operators  as 

T(Z)  =    I    T(J')(Z)  ,  (4.12) 
j=l 

where  we  have  introduced  the  definition 

T(J)(Z)  =  Vj  +  VjG0(Z)T(Z)  .  (4.13) 
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This  is  a  final  channel  decomposition  in  which  each  term  corresponds  to 
a  collision  event  where  the  last  step  is  an  interaction  through  the 
potential  V,  [(Wa67),  (Mi 72b ) ] .    It  should  be  apparent  by  comparing  the 
expressions  for  T^')  and        that  there  is  a  close  relation  between  the 
two.    Indeed,  if         is  expressed  as 


T(j)(Z)  =  T  (Z)[l  +  t(J)(Z)] 


(4.14) 


where  t'j)  is  a  distortion  term,  one  is  immediately  led  to  the  relation 


•(j)U)  =    Z  ««6A(Z)T^J(Z) 


1=1 


(4.15) 


Finally,  substituting  Eq.  (4.15)  in  Eq.  (4.14),  one  obtains 

T(J)(Z)  =  T.(Z)  +  T.(Z)GJZ)  E  ^.T^tf)  ,  (4.16) 

O  J  U  j-^    J  1 

a  set  of  three  coupled  equations.    In  matrix  form  the  result  is 


T<« 

Tl 

0     Tj  Tj 

T(D 

T(2) 

T2 

+ 

T2    0  T2 

Go 

T(2) 

T(3) 

-  ■ 

_T3 

h    T3    0  _ 

T(3) 

- 

(4.17) 

Equations  analogous  to  these  were  first  derived  by  Watson  (Wa57), 
though  Faddeev  (Fa61)  was  the  first  to  present  a  rigorous  analysis  of 
their  mathematical  properties;  as  a  consequence,  they  are  known  as  the 
Faddeev-Watson  equations.    One  important  property  of  these  equations  is 
that  they  become  connected  upon  one  iteration    and  therefore  contain  no 
unfactorable  delta  functions.    From  a  practical  viewpoint,  these  equa- 
tions are  the  most  appropriate  expressions  to  use  in  describing  scatter- 
ing processes  in  which  both  the  entrance  and  exit  channels  correspond 
to  the  breakup  channel . 
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The  Faddeev-Watson  equations  do  contain  information  on  other 
channel  processes;  however,  it  will  be  necessary  to  derive  a  set  of 
transition  operators  more  suitable  for  rearrangement  processes.  In 
order  to  accomplish  this,  we  will  begin  by  considering  a  final -channel 
decomposition  for  the  total  scattering  wave  function  [(Wa67),  (Mi 72b ) ] . 
Channel  1  will  be  singled  out  as  the  initial  channel,  for  the  sake  of 
specificity.    Accordingly,  the  total  wave  function  will  be  labeled  by 
y^*),  and  it  will  be  associated  to  an  initial  wave  function  01  through 
the  relation  (Wa67) 

=  W(Z)*  (4.18) 

where  W(Z)  is  the  wave  operator  defined  as 

W(Z)  =  ie  G(Z)  .  (4.19) 

It  will  be  assumed  throughout,  that  these  are  half-on-shell  relations, 

i.e.,  E  is  equivalent  to  the  total  energy  of  the  system,  as  described 

initially  by         Using  equations  (4.1),  (4.7)  and  (4.9)  one  can  show 
that 

6(Z)  =  GQ(Z)  +  G0(Z)T(Z)GQ(Z),  (4.20) 

which  along  with  Eq.  (4.19)  implies  that 

W(Z)  =  WQ(Z)  +  G0(Z)T(Z)WQ(Z)  (4.21) 

where  WQ(Z)  =  ie  G0(Z).    The  final -channel  decomposition  for  W  now 
follows  immediately  from  Eq.  (4.16)  and  Eq.  (4.21),  and  is 

W(Z)  =  W0(Z)  +  Z  y(J')(Z),  (4.22) 
J 


108 

where 

Y(J)(Z)  =  G0(Z)T(J)(Z)W0(Z)  .  (4.23) 
Introducing  the  Faddeev-Watson  equations  for  T^J)  into  Eq.  (4.23)  yields 

Y<j>(Z)  =  Yj(Z)  +  I  5jiG0(Z)Tj(Z)Y(i)(Z),  (4.24) 
having  defined  Y-(Z)  as 

Yj(Z)  =  G0(Z)Tj(Z)W0(Z)  =  Wj(Z)  -  W0(Z).  (4.25) 

Applying  W(Z),  as  given  by  Eq.  (4.22)  to  $p  yields  the  desired  final- 
channel  decomposition  for  the  wave  function, 

^1(+)  =  i  T(j)  (4.26) 
J 

where 

T  U)  =  Y(j)  $j    .  (4.27) 

Here  we  have  made  use  of  Eq.  (4.25)  and  the  half-on-shell  relations 

WQ(Z)*1  =  1e  S0(Z)*1  =  0  (4.28) 

and 

WjCZ)*!  =  ie  GjCZ)*!  =  fij^j  .  (4.29) 
The  Faddeev  equations  for  the  total  wave  function, 

Ti(j)  =  6ji*i +  5  Vo(z)Tj(z)Ti(i)  (4*30) 

are  now  easily  obtained  from  equations  (4.24)  and  (4.27).  The  corre- 
sponding matrix  form  of  these  equations  is 
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"*  (1)" 
1 

$1 

*  (2) 

0 

+  Go 

Tj(3) 

0 

0     T  Tj 


T2    0  T2 


T3    T3  0 


T  (1) 

^  (?) 

1 

(4.31) 


At  this  point  we  can  now  introduce  the  formal  half-on-shel 1  definition 

T^Z^tf)  =  V(1^W(Z)  (4.32) 

for  the  transition  operator  T...,  which  describes  processes  originating 
in  channel  1  and  ending  in  channel  i  (i  =  0,  1,  2  or  3).    Using  the 
definition  of  V^',  Eq.  (4.32)  may  be  written  as 


J 


V(J) 


(4.33) 


which  leads  to  the  relation 


■1 


HI      U  0      1     j=1    ij  J  o  jl  1 


(4.34) 


This  implies  the  set  of  coupled  equations 


Tu-7ii  V1*^  VjVji 


(4.35) 


which  are  known  in  the  literature  as  the  Alt,  Grassberger  and  Sandhas 
equations,  or  simply,  as  the  AGS  equations  (AT 67).    These  operators  are 
not  unique  however,  since  one  could  always  add  a  quantity  to  each  T. 
which    does  not  contribute  to  the  transition  amplitudes  in  the  on-shell 
limit  (Lovelace  [(Lo64a),  (Lo64b)]).    Another  possibility  would  be  to 
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choose  a  different  partition  of  the  various  operators  involved,  allow- 
ing one  to  construct  equations  which  put  more  emphasis  on  selective 
channels,  which  have  a  greater  physical  consequence  (Hahn  and  Watson  (Hs71), 
Kouri  and  Levin  (Ko75)  and  Tobornan  (T074)). 

We  have  now  specified  all  the  multichannel  transition  operators 
for  the  given  atom-diatom  system.    One  striking  aspect  of  the  resulting 
relations,  as  given  by  Eq.  (4.35),  is  the  similarity  between  the  expres- 
sions for  the  various  transition  operators.    In  particular  if  one  com- 
pares the  breakup  transition  operator,       ,  to  the  remaining  operators, 
Tjj  (j  =  1  to  3),  one  discovers  the  relation  (Ko77b) 


T01  ■  <!  +  TjGo>Tjl  -  SjlV1  <4-36' 


or  more  concisely 


T01  =  w/Tjl    .  (4.37) 


The  second  term  on  the  right  hand  side  of  Eq.  (4.36),  5jj60"  '  wil1  not 
contribute  to  the  on-shell  scattering  amplitudes,  i.e., 

MQ1  ■  <*0(Z)|TQ1(Z)|*1(Z)>  (4.38a) 

=  <»jC(Z)|TJ1(Z)|»1(Z)>    ,  (4.38b) 


where 


$.c  =  Wd$0    .  (4.39) 


In  this  last  expression,  <±>-c  is  a  wave  function  (describing  the  inter- 
acting  pair  in  channel  j  as  a  continuum  state.    The  superscript  "c"  has 
been  introduced  to  distinguish  <f>.c  from  $i  which  describes  the  inter- 
acting  pair  in  channel  j  as  a  bound  state.    What  this  means  physically 
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is  that  one  can  describe  the  breakup  process  as  one  would  a  rearrange- 
ment process,  but  the  interacting  pair  in  the  exit  channel  must  be 
assumed  to  lie  in  a  continuum  state. 


In  practical  applications,  the  exact  treatment  of  multichannel 
scattering  processes  via  the  coupled  AGS  equations  for     ^  is  generally 
not  possible.    This  is  particularly  true  at  high  energies,  where  the 
number  of  accessible  states  in  each  given  arrangement  channel  becomes 
quite  large.    Fortunately,  the  form  of  the  AGS  equations  suggests  that 
one  might  be  able  to  avoid  these  difficulties  by  considering  their 
iterated  form.    We  will  now  discuss  the  implications  of  this  approach 
as  it  pertains  to  atom-diatom  collisions. 

If  one  now  iterates  Eq.  (4.35),  one  obtains  the  series 


The  physical  meaning  of  this  series  becomes  clear  if  one  keeps  in  mind 
the  nature  of  the  operators  involved.    Consider  for  example  the  first 
term  in  Eq.  (4.40).    From  the  discussion  in  the  previous  section  it  is 
obvious  that  this  term  does  not  contribute  to  either  inelastic  or  break- 
up processes.    Therefore,  it  plays  a  role  only  in  the  case  of  rearrange- 
ment scattering.    Since  this  term  does  not  contain  any  interaction  terms, 
it  is  often  called  the  "spectator  stripping"  term  [(Mi72b),  (Mi72c)]. 
The  second  term  in  the  series  contains  only  the  operators  T-,  which  are 
two-body  transition  operators  in  a  three-body  space.    Hence,  these  terms 


3.    The  Multiple-Collision  Expansion 


-1 


+  5  6ijfijiTj  + 

j 


(4.40) 
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describe  processes  in  which  the  incoming   atom  interacts  with  only  one 
of  the  target  atoms  in  the  diatomic,  while  the  third  atom  acts  as  a 
spectator.    The  higher  order  terms  in  Eq.  (4.40)  would  then  correspond 
to  double  and  higher  order  encounters.    The  role  of  the  CS-^}  in  all 
these  terms  is  quite  important,  for  they  exclude  the  possibility  of 
terms  leading  to  disconnected  diagrams  and  rule  out  terms  which  do  not 
contribute  to  the  given  process  on  physical  grounds  (Ch71).    Because  of 
this  conceptual  picture,  this  series,  Eq.  (4.40),  is  known  as  a  multiple- 
collision  expansion  [(Ch71),  (Mi 75b) ] . 

Each  term  in  the  multiple-collision  expansion  has  now  been  asso- 
ciated with  a  direct  mechanism  for  an  atom-diatom  scattering  event. 
What  one  hopes  is  that  only  a  limited  number  of  these  terms  give  a  signi- 
ficant contribution  to  the  scattering  amplitude.    The  success  of  results 
obtained  by  truncating  the  expansion  would  then  be  indicative  of  the 
types  of  mechanisms  involved  in  a  given  physical  system.    For  reactive 
scattering,  there  already  exists  ample  evidence  to  corroborate  the 
existence  of  a  stripping  mechanism  in  certain  atom-diatom  reactions, 
e.g.,  Minturn,  Datz  and  Becker  (Mi 67 )  have  shown  that  a  classical 
spectator  stripping  model  properly  describes  the  location  of  peaks  in 
the  product  distribution  of  various  alkali-atom-halogen-molecule  reactions 
[Cs  +  Br2>  K  +  (I2  or  Br2)].    Within  the  context  of  the  multiple  colli- 
sion expansion,  similar  success  was  obtained  by  Yuan  and  Micha  [(Yu76a), 
(Yu76b)],  who  used  the  spectator  stripping  model  to  study  the  systems 
[Ar+,  (H2,  D2,  HD)]  and  [K  +  (Ig,  Br2,  IBr)].    It  should  be  pointed  out 
that  the  expression  for  the  transition  amplitudes, 

MSfr  =  <*i |G0~1(Z)|*1>  ,  (4.41) 
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obtained  within  the  stripping  model  is  equivalent  to  the  Qppenheimer 
(0p28),  and  Brinkman,  Kramers  approximation  (Br30), 

MT  =  Mnr  =  <«iHMV  .  (4-42) 

which  was  introduced  in  the  study  of  electron  capture  processes  in 
atomic  physics  (Mc70).    The  extension  of  these  studies  to  include 
multiple-collision  effects  has  generally  been  pursued  only  within  the 
framework  of  classical  mechanics  [(Ba64),  (Su68),  (Ge69),  (Ma74), 
(Ma76)],  but  the  success  of  these  models  has  been  impressive  consider- 
ing their  simplicity.    We  defer  the  study  of  single  and  higer  order 
terms  to  the  next  chapter,  where  a  more  detailed  analysis  will  be 
given. 


CHAPTER  V 
THE  SINGLE  COLLISION  APPROXIMATION 
In  this  chapter  we  will  investigate  the  single  collision 
approximation  as  it  applies  to  inelastic  and  dissociative  atom- 
diatom  collisions.    As  can  be  seen  from  Eq.  (4.35),  the  multiple 
collision  expansion  for  the  transition  operators  of  interest  is 

Tn  =  T2  +  T3  +  T2GQT3  +  T3G0T2  +  ...  (5.1) 

for  inelastic  scattering,  and 

T01  =  V1  +  T2  +  T3  +  T2GoT3  +  T3GoT2  + 

+  W2  +  TlGoT3  +  •  •  *  <5-2> 

for  dissociative  scattering.    The  single  collision  approximation  to 
the  scattering  amplitudes  would  then  be 

Mr!   =    <$l'lT2  +  T3lV  ^5-3) 

and 

M01}    =    <$0|T2  +  '  (5'4> 

where  the  higher  order  terms  in  equations  (5.1)  and  (5.2)  have  been 
neglected.    Note  that  the  term  G0_1  in  Eq.  (5.2)  does  not  contribute  to 
the  breakup  amplitude  in  the  on-shell  limit.    In  the  following  sections 
we  will  discuss  the  practical  implementation  of  equations  (5.3)  and 
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(5.4)  as  well  as  the  factors  which  determine  their  validity. 

1.    Description  of  Channel  States 
For  each  arrangement  channel  i,  i  =  1,  2  or  3,  we  shall  introduce 
the  set  of  Jacobi  coordinates  ft.  and  r*,  along  with  the  corresponding 
conjugate  momenta  P-  and  p..    The  momenta  P^  and  pi  describe,  respec- 
tively, the  relative  momentum  and  the  diatomic  relative  momentum  in 
the  given  channel  i.    A  definition  of  all  these  coordinates  in  terms 
ordinary  laboratory  coordinates  for  the  three  atom  system  is  given  in 
Appendix  IV. 

Working  with  the  given  relative  coordinates  in  the  center  of  mass 
system,  we  can  describe  the  free  three-particle  motion  in  channel  0  by 
the  plane  wave 

<Riri|Pipi>  =  <  fyl^x  ?i|Pi>  .  (5-5) 


where 


<  R^Pf  =  (2tt)'3/2  e1'  Pi  '  Ri  (5.6) 


and 


<;ilV-<2*r3/2elPl'ri     •  (5-7) 

Note  that  any  one  of  the  channel  coordinates  can  be  used  to  describe  the 
system,  since 

Pj  •  Rj  +  Pj  •  rj  =  P1  •  R1  +  pn-  •r1    CI,  j  =  1,2,3),  (5.8) 

and  thus 

<VjifjV =  <ViiVi>-  •  (5-9) 
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Equation  (5.8)  can  be  proved  using  the  relations  in  Appendix  IV.  This 
equivalence  amongst  the  various  channel  coordinates  can  be  stated 
formally  by  the  relations 


<Rjrj|Riri>  =  6[Rj  -  RjfRj  ,ri  )]6[rj  -  r^R^r^]  (5.10) 
and 

<  Vj^V  =  &[h  '  PjCPi-Pi)^[Pj  -  PjOW1  ■  (5-n) 

->■->■->■ 

where  a  quantity  such  as  R^R^r^  denotes  that  value  of  Rj  obtained 
from  its  relation  to  R^  and  ri ,  as  given  by  Eq.  (A4.4)  in  Appendix  IV. 
Alternately,  one  could  derive  these  last  two  equations  by  making  use  of 
equations  (5.5)  through  (5.8)  and  the  identity  operators  (Re77) 

1  =  /  dRi  dr1  |  R1  ri><  ^  ri  |  (5.12) 

and 

1  =  /  dPi  dpi  |  Pi  ^  |    .  (5.13) 

The  free  motion  in  the  remaining  channels  can  be  described  by 
(Mi  75) 

*i  =  Xi(Pi)  *i  (n^m-j  )       (i  =  1.  2  or  3),  (5.14) 

->- 

where  X^Pj)  is  a  relative  motion  state  and  <f>.j  is  a  vibrotor  state 
characterized  by  the  quantum  numbers  n.- ,  jj  and  m,-..    In  the  coordinate 
space  representation  one  has 

<Rir.|$.>=  <Ri|x1(P1)><r||#i>     ,  (5.15) 
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where 

<  Ri Ixi (Pi )>  =  <  RilV  ■  x  (Pi .R)  (5-16) 

and 

<ril*i>  =  ^(n^  ji  nij_;  r.)      .  (5.17 
Similarly,  in  the  momentum  representation  one  obtains 

<Pi'       |  *j>  =  <Pi'  l^xPi'  1*^  (5.18) 


or 


<Pi'       |*i>  =  6 (Pi'  -  P-i  )*i  (ni ,  ji5  mj  ;       )  , 


where 


(5.19) 


J'i;  Pi'  }  =  /  d?i  'Pi'  IV<7il*i>  {5'20) 

is  a  diatomic  momentum  distribution  amplitude. 

2.    Inelastic  Scattering 

The  single-collision  inelastic  scattering  amplitude,  was 

given  by  Eq.  (5.3).    We  will  denote  each  term  in  this  expression, 

(s) 

<  $j'  |Tj|^>,  by  Mj.j(j);  and  for  the  sake  of  specificity,  we  will 
study  the  j  =  3  term.    Working  in  the  momentum  representation,  we  see 
that 

h[S,\  (3)  =  /  dQ3'  dq3'  dQ3  dq3  <  ^  |Q3'  ^  ><  Q3'  ^  [  T3 1  Q3  q*3>  x 

x  <  Q3  q3l*1>    •  (5.21) 
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By  convention,  P.  and  p.  will  be  used  to  denote  the  momentum  charac- 

■+  -y 

terizing  an  initial  or  final  state        and  the  variables  Q^.  and  will 
be  reserved  for  the  intermediate  momenta  (Mi75b).    Using  the  definition 
of  T3,  Eq.  (4.11),  we  see  that 

^3  q3  lT3lq3  V  *  S(Q3  "  Q3>  *3  (q3  •  "3-  E  "  2H3  >•  (5'22) 

where  t^  is  a  two-body  t-matrix  describing  the  interaction  between 
atoms  A  and  B.    M3  is  the  channel  reduced  mass  defined  by  Eq.  (A4.8) 
in  Appendix  IV.    Note  that  the  energy  argument  of  t3  corresponds  to 
the  total  energy  minus  the  kinetic  energy  of  the  spectator  atom.  From 
Eq.  (5.22)  one  can  also  ascertain  the  fundamental  difference  between 
the  three-body  quantity  T"3  and  the  two-body  quantity  t^.    On  the  energy 
plane  the  pole  structure  of  t3  is  replaced  by  branch  points  in  the 
structure  of  T"3  (Lo64a).    Matrix  elements  such  as  <  Q3  q3l*1>  may  be 
evaluated  via  equations  (5.11)  through  (5.19),  yielding  the  result 

<Q3  q3l*1>  =  6(q3  -  q3(P1,Q3))<q1(Pli  Qjfof  .  (5.23) 

(s) 

Thus,  the  expression  for  M  ,    (3),  Eq.  (5.21),  leads  to 

i  >  q  2 

MJ.J  (3)  -  /  dQ3  <*1,|q1'>  t3  (q3'  .  q3;  E  -  JL  )  <q1|<K1>>  (5.24) 


where 


q3  =  q3(Q3,  Pj)  -  -Px  -          Q3  (5.25a) 

q3'  =  q3(Q3,  P{  )  =  q3  -  APj  (5.25b) 
->      ->■->->        mp   ->  ->- 

ql  =  ¥Q3'  Pl)=i^Pl+  Q3  <5'25c> 
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?i  ■51Mi,|,V'Eii  (5-25d) 

and  APj  =  Pj    -P.  is  the  momentum  transfer.    As  can  be  seen  from  Eq. 
(5.24),  m|?J  (3)  has  a  very  simple  structure.    All  that  is  needed  in 
its  evaluation  is  the  product  of  the  initial  and  final  momentum  distri- 
bution amplitudes,  along  with  the  two-body  t-matrix.    This  simplicity, 
however,  is  quite  deceiving.    The  integrand  in  Eq.  (5.24)  actually  has 
a  multi-center  nature,  as  can  be  seen  from  equations  (5.25a)  through 
(5.25d).    In  essence  the  problem  we  are  facing  is  the  familiar  one  deal- 
ing with  the  angular  momentum  analysis  of  an  atom-diatom  scattering 
process.    The  usual  approach  taken  in  solving  problems  of  this  nature 
involves  the  introduction  of  various  partial  wave  expansions  for  each 
of  the  quantities  involved  [(Ah65),  (Ba69a),  (El 69) ,  (Ha70),  (0m64)]. 
However,  this  does  not  seem  practical  for  high  energies,  which  require 
a  large  number  of  partial  waves.    Another  possibility,  would  be  that  of 
using  a  three-dimensional  quadrature  scheme.    This  is  particularly 
attractive,  since  the  integrand  of  Eq.  (5.24)  contains  the  product  of 
the  momentum  wave  functions  <  <j>^'  |q^  >  and  <  qj^«>,  which  would  limit 
the  range  of  integration  over      to  some  overlap  region  between  the  two 
amplitudes  (Mi 75b).    There  is  one  major  problem  though,  and  it  is  the 
idea  of  a  three-dimensional  quadrature  itself.  Multidimensional 
quadratures  are  generally  somewhat  unreliable  [(Ba74),  (St67)].  This 
however,  is  an  open  question  requiring  further  investigation.  Finally, 
there  remains  the  possibility  of  introducing  approximations  into  the 
evaluation  of  m|^|  (3).    In  this  regard,  we  will  now  consider  what  is 
known  in  the  literature  as  the  "peaking"  approximation  [(Ro67),  (Mc70), 
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(Ek71),  (Ek74) ] .    Within  this  approximation  one  factors  the  t-matrix, 
ty  out  of  the  integral  sign  in  Eq.  (5.24)  yielding  the  expression 

M{!;Pk)(3)  =  t^    /  dQ3  <*{  |qx'  ><\\*y>  .  (5.26) 

where  tPk  denotes  that  value  of  t,  obtained  for  some  suitably  chosen 
3  3 

value  of  Q,,  say  QPk.    The  criteria  for  selecting  oj?k  will  be  discussed 
5  3  J 

later  in  Section  VI. I.     The  assumption  one  has  introduced  here  is  that 
t,  is  a  smooth  and  slowly  varying  function  of  Q?  within  the  range  of 
integration  determined  by  <  <j>  '  Iq^  >  and  <  q^cf)^,  which  on  the  other- 
hand,  are  highly  oscillatory  functions  (Ek71).    Physically,  the  peaking 

approximation  implies  that  we  have  factored  the  transition  amplitude, 

( s)  Pk 
M:,J  (3),  into  a  term  t3   which  contains  information  on  the  dynamics  of 

the  collision  process,  and  a  term 

Fl'l  "    /  ^3  <<J)i'  l?l  ><51l*1>  (5.27) 

called  a  target  "form  factor,"  which  contains  information  characterizing 
the  size  and  shape  of  the  target,  along  with  information  on  the  initial 
and  final  momentum  distributions  of  the  particles  involved,  i.e.  the 
restoring  forces  of  the  diatomic  target.    We  see  then  that  the  role  of 
the  potential  for  the  target  diatomic,  in  this  case  V3,  is  one  of  deter- 
mining the  initial  and  final  momentum  distributions  [(Ek71),  (Mi 75b) ] . 

Consider  now  the  validity  of  the  single  collision  approximation. 
We  know  that  in  general  multiple  collision  expansions  are  applicable 
only  at  high  energies  [(Ch71),  (Ch73),  (Mi75b)].  From  a  mathematical 
point  of  view,  the  convergence  properties  of  expansions  such  as  (5.1) 
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and  (5.2)  can  only  be  ascertained  by  detailed  consideration  of  the 
potentials  involved.    In  practice,  this  turns  out  to  be  quite  difficult; 
and  so,  one  must  resort  to  computational  trial  and  error.    One  can, 
however,  put  forth  various  rules  based  on  physical  considerations  which 
give  some  indication  of  the  validity  of  the  multiple  collision  expansion 
and  in  particular  that  of  the  single  collision  approximation.    We  note 
that  equations  (5.1)  and  (5.2)  are  expansions  in  terms  of  the  quantities 
T..  and  GQ.    Because  of  the  definition  of  G0,  (Z  -  H0)_1,  it  is  apparent 
that  these  expansions  should  indeed  exhibit  better  convergence  properties 
at  higher  energies.    One  cannot  be  definitive  however,  since  the  operators 
{T.,-}  are  present  in  each  term.    In  general  these  operators  would  not  be 
expected  to  cause  difficulty  for  most  potentials  of  chemical  interest 
[(Ch71),  (Mi75b )].  Since      contains      to  all  orders,  it  remains  well 
behaved  even  for  singular  potentials.    As  to  the  validity  of  the  single 
collision  approximation,  one  can  see  from  a  classical  mechanical  picture 
that  such  a  collision  process  must  be  impulsive.    In  other  words  the 
projectile  must  strike  and  leave  the  interaction  region  fast  enough  so 
that  the  restoring  forces  of  the  target  diatomic  do  not  lead  to  sequen- 
tial collisions.    This  latter  assertion  implies  that  for  the  impulse 
approximation  to  be  valid,  the  projectile  energy  must  be  larger  than 
the  internal  excitation  energy  of  the  target.    Considering  the  expres- 
sions for  T,,,  T?1  and  Tn,  obtained  from  Eq.  (4.3S),  it  follows  that 


Z       Z    o    1    o    11       2    o   3  o  31 


T,  +  T,  Gn  T.  Gn  T„ 
3       3    o    1    o  11 


+  T3  GQ  T2  G0  Jn 


(5.28) 


Using  the  relation 
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61  -  GQ  =  G0  T1  GQ  (5.29) 

which  follows  from  the  definition  of  G,  and  T.,  and  equations  (4.7) 
and  (4.8),  we  see  that  Eq.  (5.28)  is  equivalent  to 

Tn  =  T2  +  T2(Gj  -  Go)Tn  +  T2(G3  -  Go)T31  ♦ 

T3  +  T3(G!  -  G0)Tn  +  T3(G2  -  G0)T21  ,  (5.30) 

a  formal  expression  that  implies  the  weaker  the  various  potentials  V. , 
the  smaller  the  contribution  from  the  corresponding  multiple  collision 
terms.    Other  important  considerations  regarding  multiple  collisions 
relate  to  the  size  and  geometry  of  the  target  involved  [(Me51),  (Mi78)], 
If  dg  and      are  the  diameters  of  the  atoms  in  the  diatomic;  and  RgQ, 
the  average  value  of  r,;  then,  multiple  collision  terms  are  less  impor- 
tant when 


and 


Pu  RBC  »  1  (5.31) 


(dB  +  dc)/2  <RBC  (5.32) 


(Mn 78) .    The  relative  masses  of  the  atoms  in  the  atom-diatomic  system 
are  also  very  important.    The  likelihood  of  multiple  collisions  for 
inelastic  scattering  diminishes  if 


mn  mo 

^rir  <  1  (5-33) 
m£  M 

and 
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[(Sg66),  (Se69),  ( E k 7 1 ) ,  (Ek74)].    A  further  discussion  of  the  impor- 
tance of  the  above  mass  factors  will  be  taken  up  in  chapter  VI  for  the 
case  of  col  linear  scattering. 

We  will  now  consider  some  kinematical  rules  which  follow  from 
Eq.  (5.24),  i.e.  the  definition  of  M**j  (3).    It  has  been  pointed  out 
(Mi 75 )  that  the  functions  <  <^'!q^'  >  and  <  Q^l^  are  each  localized 

within  a  given  spherical  region;  implying  that  the  degree  of  overlap  is 

(s) 

a  measure  of  the  size  of  the  transition  amplitude  (3).    A  good 

quantitative  estimate  of  the  volume  of  each  sphere  can  be  obtained  by 

taking  |q, I       and  !q  'I         equal  to  their  classical  values;  and  thus 
'^l'max  1  max 


"ll  <  1*1 'max  ~  r2mlwini]  (5-35a) 


and  jq^'l  <  1^'^  =   [2miWinf]  1/2  (5.35b) 


where  W^i  and  W  nf  are  initial  and  final  classical  vibrational  energy 
of  the  diatomic  BC.    Based  on  this  analysis  a  kinematical  construction 
can  be  made  as  illustrated  in  Fig.  (5-1).    A  concise  statement  of  the 
overlap  volume  may  be  written  as  (Fe71) 

mv(3)  =  /  dQ3  etlVl^-l^'Mefl^l^-lqJ]  .  (5.36) 

where  q^  and  q    are  given  by  equations  (5.25c)  and  (5.25d).    The  func- 
tion 9(x)  is  known  as  the  Heaviside  step  function  and  is  defined  as 
(Cu75) 

1     x  >  0 


e(x)  =  { 


(5.37) 

0     x  <  0 


Figure  (5-1)    Kinematical  construction  for  in- 
elastic scattering. 
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The  integral  in  Eq.  (5.36)  is  a  simple  exercise  in  elementary  geometry 
and  is  thus  very  easy  to  evaluate  (Se70).    In  Chapter  VII  we  will  see 
that  the  overlap  volume,  m^,  can  in  many  instances  predict  the  relative 
shape  of  the  energy  loss  spectrum  for  inelastic  scattering.  However, 
from  Eq.  (5.24)  or  Eq.  (5.26),  it  is  apparent  that  one  should  really 
worry  about  the  relative  phases  and  magnitude  of  <  <j>j'|q  '>  and 
<qj|<f>.>  in  the  overlap  region.    A  simple  extension  of  the  present  model 
would  be  to  use  the  semiclassical  momentum  distributions  inside  Eq. 
(5.36),  as  recently  proposed  (Mi 77 ) .    This  would  still  not  address  the 
problem  of  the  phases,  but  it  would  lead  to  a  quantity  which  may  be 
used  to  analyze  the  possible  mechanisms  in  an  atom-diatom  collision. 
Such  a  model  is  attractive,  since  it  would  provide  a  simple,  inexpensive 
way  for  the  theoretician  or  experimentalist  to  gain  insight  about  energy 
transfer  and  reactive  processes  in  atom-diatom  collisions.    In  any  case, 

the  analysis  that  led  to  Eq.  (5.36)  may  be  used  to  determine  the  region 

(s)  (s  Pk) 

of  integration,  if  one  were  to  evaluate  MJ,J(3)  or  Mj,J    '(3)  by  a 

multidimensional  quadrature  scheme. 

3.    Dissociative  Scattering 
Following  the  same  procedure  used  in  the  previous  section  one  can 
show  that  for  dissociative  scattering  Eq.  (5.4)  leads  to 


-  t3  (p3f,  q3.;  E      )  <  5^.  |^>  , 


(5.38) 


f 


where 


(5.39) 
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mAmc      +  mA 
q  .  +    ~AJ^   P..  "    —   P1f  (5-40) 

31  11       nlABmBC      lf       %  lf 


and 


5  =  q  -  AP,  (5.41) 
H3f        M3i  1 


5i  =  5if  "  ST  A^i   •  (5'42) 

i  it  mB(. 


Equations  (5.40)  through  (5.42)  follow  from  the  relations  in  Appendix 

IV.    One  does  not  have  an  integral  in  Eq.  (5.38)  as  was  the  case  for 
(s) 

Mj.j(3),  since  the  final  state  $0  is  characterized  by  a  product  of 
plane  waves  as  given  by  Eq.  (5.5).    Furthermore,  Eq.  (5.38)  is  simpler 
than  Eq.  (5.24),  since  it  contains  only  the  half-on-shell  t-matrix. 
This  simplification  however,  does  not  help  as  much  as  one  would  think. 
Experimentally,  Eq.  (5.38)  would  imply  the  need  for  a  coincidence  study, 
one  in  which  two  of  the  exit  particles  are  detected  simultaneously 
[(Sc74a),  (Sh77)].    This  is  something  that  has  not  been  done  so  far 
[(Sc73),  (Sc74b),  (Sc75)],  and  thus  one  must  average  IMq^|2  over 
to  obtain  a  proper  comparison  with  experiment  (Sc73).    Because  of  energy 
conservation,  only  an  average  over  the  orientation  p^  would  be  needed. 
We  will  address  only  the  simpler  problem  of  colli  near  scattering,  and 
will  postpone  further  discussion  on  this  matter  to  the  next  chapter. 

In  section  IV. 2  we  presented  an  alternate  expression  for  Mg^,  Eq. 
(4.38b).    Within  the  single  collision  approximation  this  relation  leads 
to 

M0l'C)  =  <$1  'T2  +  T3l*i  >    '  <5'43) 
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The  latter  expression  would  lead  to,  e.g., 

c  "* 

an  expression  analogous  to  Eq.  (5.24),  where  however,  <      |q    >  is  a 
continuum  momentum  wave  function  for  the  diatomic  BC  in  channel  1. 
The  consequence  of  Eq.  (5.43)  will  also  be  explored  in  the  next  chapter 
in  connection  with  colli  near  scattering. 


CHAPTER  VI 
COLL I NEAR  SCATTERING 

In  the  previous  chapter  we  proposed  the  application  of  the  single 
collision  approximation  to  the  study  of  inelastic  and  dissociative 
atom-diatom  scattering  processes.    We  will  now  consider  this  approxima- 
tion for  the  case  of  collinear  scattering,  which  offers  a  simple  and 
inexpensive  check  on  the  validity  and  limitations  of  such  an  approach. 
There  is  also  the  advantage  of  having  an  extensive  literature  from 
which  one  can  ascertain  the  relative  merit  of  the  present  many-body 
treatment  when  compared  to  other  quantum  mechanical,  semiclassical  or 
classical  treatments  of  the  collinear  scattering  problem.    An  extensive 
survey  of  collinear  scattering  may  be  found  in  the  reviews  by  Shin  (Sh76), 
Rapp  and  Kassal  (Ra69),  and  Takayanagi  [(Ta63),  (Ta65)]. 

From  a  physical  standpoint  it  is  hoped  that  a  study  of  collinear 
scattering  will  provide  useful  information  about  translational-vibra- 
tional  energy  transfer  in  atom-diatom  collisions.    The  collinear  config- 
uration is  in  general  the  most  effective  one  for  vibrational  excitation 
[(Sh76),  (Ra69)].    However,  one  must  be  careful  in  generalizing  the 
results  obtained  from  such  a  model.    Collinear  (1-D)  and  full  three- 
dimensional  (3-D)  dynamical  calculations  differ  in  two  very  important 
aspects  (Be74b):    the  amount  of  phase-space  available  in  1-D  collisions 
is  less  than  that  for  3-D  collisions,  i.e.  1-D  scattering  calculations 
suffer  from  a  "dimensionality-bias";  secondly,  steric  factors  play  a 
more  important  role  in  3D-coll isions.    The  importance  of  the  role  played 
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by  either  of  these  two  factors  will  then  determine  the  validity  of  1-D 
scattering  calculations.    Because  of  these  limitations  in  col  linear 
studies,  we  will  pursue  this  study  only  as  a  test  leading  to  an  evalua- 
tion of  the  single  collision  approximation,  the  peaking  approximation 
and  the  importance  of  multiple  collision  terms. 

1.    Formulation  of  the  Colli  near  Scattering  Problem 
For  the  sake  of  specificity,  we  will  select  channel  1  as  our 
initial  arrangement  channel,  i.e.  A  +  (B  +  C).    At  present,  a  harmonic 
oscillator  model  will  be  used  to  represent  the  diatomic  (B-C).  Later 
in  Section  VI. 4   where  dissociation  is  considered,  a  Morse  oscillator 
will  be  introduced.    The  interaction  potential,  V2,  between  the  projec- 
tile atom  and  the  ncnstruck  atom  will  be  neglected. 

In  order  to  simplify  our  analysis  of  the  1-D  scattering  problem 
and  compare  to  other  work  in  the  literature,  we  shall  adopt  the  follow- 
ing scaled  Jacobi  coordinates: 


xl  =  [(m1k)1/2/Ti]1/2  ^  (6.1) 


and 


where 


and 


Xj  =  [(n1k)1/2/ft]l/2^*1  ,  (6.2) 

c 


*l  =  XC  -  XB 


(6.3) 


Xl  =  XA  "  (mAxB  +  mCXC)/mBC  (6'4) 


are  the  Jacobi  coordinates  defined  in  Appendix  IV  as  they  apply  to  1-D 
scattering.    In  equations  (6.1)  and  (6.2)  we  have  introduced  the  diatomic 
reduced  mass  m,  =  (mDnu/mDr)  and  the  harmonic  oscillator  force  constant  k. 

1  d  L  dL 
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The  particular  form  of  equations  (6.1)  and  (6.2)  allow  us  to  work  in 
energy  units  of  hu, 

u  =  [k/m^172  .  (6.5) 

The  resulting  Schrodinger  equation  for  the  three-body  system  in  the 
center-of-mass  frame  is  (Ra69) 

{-  \  ^~~o  ~  \  ^-o  +  V,  (X, )  +  V3(X,  -  X,)}*,  =  CT,  ,  (6.6) 
2    3X:2      2    3Xl2       11         J    1       1      1  1 

where 


mAmc 


(6.7) 


Vx(Xl)  =  \  (Xj  -  xleq)2  ,  (6<8) 

and  V3(X1  -  x^  is  the  interaction  between  the  projectile  A  and  the 
struck  atom  B.    M  is  the  total  mass  of  the  system  (mA  +  mg  +  mc). 
Note  that  Vj  and  V3  are  not  the  same  as  before.    They  correspond  to  a 
modification  of  the  original  potentials  consistent  with  the  new  set  of 
coordinates  and  units.    A  pictorial  representation  of  the  coordinates 
used  is  given  in  figures  (6- la)  and  (6-lb).    From  Eq.  (6.6)  and  Fig. 
(6- lb)  it  should  be  clear  that  we  have  transformed  our  former  three- 
body  problem  into  one  where  an  oscillator  of  unit  mass,  vibrating  about 
an  equilibrium  position,  is  struck  by  a  particle  having  an  effective  mass 
m;  or  equivalently,  into  a  problem  where  a  harmonic  oscillator  having  a 
particle  of  mass  and  one  of  unit  mass,  is  struck  by  a  particle  of  mass 
m.    We  note,  as  pointed  out  by  Secrest  and  Johnwi  (Se66),  that  the  mass 
parameter  m  corresponds  to  many  different  atom-diiatom  problems    and  is 
thus  useful  in  studying  the  importance  of  mass  effects  in  atom-diatom 
coll isions. 


Figure  (6-1)    Coordinates  used  in  col  linear 
scattering. 

(a)  Jacobi  coordinates 

(b)  Scaled  Jacobi  coordinates 
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Using  the  physical  picture  given  above,  we  can  now  state  the 
single  collision  approximation  for  inelastic  scattering  as 

Mi'i  =  Mri  (3>  =  <*! •  I t3 | *1>  ,  (6.9) 

where  we  have  used  Eq.  (5.3)  and  have  dropped  the  term  Mj , .  (2),  since 
^2  is  neglected  in  the  given  model.    In  the  momentum  representation, 
Eq.  (6.9)  leads  to 

(s)  i  Qi2 

Mri  =  /  dQ3  <c?lt  Iq^  >  t3  (q3  ,  q3;  E  -  ^-  )<qi|^>  ,  (6.10) 

the  analog  of  Eq.  (5.24)  obtained  in  the  three-dimensional  case.  The 
various  channel  momentum  coordinates  present  in  Eq.  (6.10)  can  be 
obtained  from  Appendix  IV,  if  one  considers  a  three  particle  system 
with  mass  mA  =  m,  mg  =  1  and  mc  =  °°  ,  as  illustrated  in  Fig.  (6- lb) . 
In  channel  1,  Pj  and  p1  are  the  conjugate  momenta  to      and  x,, 
respectively.    The  result  for  the  coordinates  q3,  q3'  ,       and  q^  in 
Eq.  (6.10)  is  analogous  to  that  given  by  equations  (5.25a)  through 
(5.25d),  but  with  the  masses  as  given  above.    The  peaking  approximation 
to  Eq.  (6.10)  is 

MriPk>  =  ^  <n'le1APl-l|n>  (6.11) 

where  x1  is  the  relative  position  operator  for  the  B-C  pair,  n  and  n' 
are  the  initial  and  final  vibrational  quantum  numbers  characterizing 
(f>j  and  (Jy,  respectively.    The  two-body  t-matrix  t3^Pk^  is  factored  out 
from  Eq.  (6.10)  at  some  value  of  Q3*.    Once  this  choice  is  made  q3'*  and 
qj  are  obtained  from  equations  (5. 25a)  through  (5.25d).    Noting  that  the 
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Gaussian  envelopes  of  the  harmonic  oscillator  functions  <n • | q^'  >  and 
<q1|n>  in  Eq.  (6.10)  peak  when  and  q.^  equal  zero,  we  have  chosen 
Q3*  as  (Ek74) 

Q3*  =  \  [Qgtq^    =  0)  +  Q3(q3  =  0)]  ,  (6.12) 
which  is  equivalent  to 

Q3*  =  "  \  (P1  +  V'  (6'13) 

as  can  be  seen  from  equations  (5.25c)  and  (5.25d).    The  choice  of  Q3* 

is  in  general  quite  arbitrary  as  long  as  it  remains  within  the  overlap 

(s  Pk) 

region  of  the  two  momentum  amplitudes.    We  have  found  M\j       to  be 
quite  insensitive  to  the  choice  of  Q3*. 

2.    Two-Body  t-Matrix  for  the  One-Dimensional  Scattering  Problem 
In  the  present  treatment  of  collinear  scattering  we  shall  always 
maintain  a  close  parallel  to  the  3-D  collision  problem.    Because  of 
this,  we  will  deal  only  with  symmetric  potentials,  i.e.  those  such  that 
V(-x)  =  V(x).    Furthermore,  we  will  assume  the  presence  of  a  hard  core 
in  order  to  take  into  account  the  physical  impenetrability  of  the 
particles  involved. 

The  analogue  to  Eq.  (2.6),  the  differential  equation  satisfied 
by  the  off-shell  wave  function,  is 

dx2  q 

where 

<x|q>  =  (2tt)"  1/2  e1qX    .  (6.15) 
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Introducing  the  definitions 


0    if    [sgn(q)sgn(x)]  =  1 
qX  tt    if    [sgn(q)sgn(x)]  =  -1     ,  (6.16) 


r  =  | x  i  (6.17) 
1^(9)  =  [cos(e)]1  (6.18) 


and 


f^qr)  =  < 
one  can  show  that 


cos(qr)  if  1  =  0 

sin(qr)  if  1  =  1     ,  (6.19) 


1 

Z 

1=0 

and 


1  1 

elclx  =  I    i1  ru  [cos(eqx)]  f  (qr)  (6.20) 


1/2    ^  1 

d)(q,  k  ;  x)  e  <x|cj  (z)>  =  (2tt)~  '      I  i  V  [cos(eqx)]io,  (q,kr;r), 
t  1=0      1  1  L 

(6.21) 

where 

2 

[  ~j  +  (^EkE2  -  v(r))]  ^(q,  k£;  r)  =  (a^2  -  q2)f](qr)  (6.22) 

and  v(r)  =  2mV(r).    The  function  n,  takes  the  place  of  the  Legendre 
polynomial  P]  and  f]  that  of  the  Riccati-Bessel  function  j]  found  in 
the  partial  wave  expansion  for  the  three-dimensional  off-shell  wave 
function,  <f|w+(z)>.    However,  the  expansions  used  in  equations  (6.19) 
and  (6.20)  are  not  partial  wave  expansions;  they  are  an  artifact  of  the 
way  we  have  chosen  to  solve  for  the  off-shell  wave  function.  Alternative 
formulations  of  one-dimensional  Schrb'dinger  equations  in  terms  of 
parity  expansions  or  in  terms  of  phase  shifts  could  also  be  pursued 
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[(Li73),  (Eb65),  (Ka61)].    We  have  chosen  the  above  formulation  because 
of  its  close  parallel  to  the  treatment  found  in  Chapter  II.    In  using 
this  approach  however,  one  must  give  close  consideration  to  the  boundary 
conditions  placed  on  w, .    If  the  potential  does  net  contain  an.  infinitely 
repulsive  hard  core,  transmission  from  one  side  of  the  potential  barrier 
to  the  other  must  be  taken  into  consideration.    Since  we  have  assumed 
the  presence  of  a  hard  core,      must  satisfy  the 

a),     ~       0  (6.23) 
l   r  +  rc 

boundary  condition  at  the  hard  core  position  rc  and  the  asymptotic 
boundary  condition 

f  (qr)  +  est  elkr    if  E  >  0 

wl  < 
r->°° 

f  (qr)  +  est  e"kr    if  E  <  0    .  (6.24) 

Note  that  when  the  asymptotic  boundary  of  the  full  off-shell  wave 
function  u>(q,  k£;  x)  is  considered,  equation    (6.21)  takes  properly  into 
account  the  direction  of  the  incoming  wave. 

The  t-matrix  may  also  be  expanded  in  a  form  analogous  to  the 
partial  wave  expansion  of  the  three-dimensional  t-matrix.    Making  use 
of  the  definition  of  t, 

+00 

t(q\  q;  E)  =  J  dx  <q'|x>  V(x)  <x|ioQ(z)>  ,  (6.25) 

-00  ' 

and  equations  (6.15),  (6.20)  and  (6.21)  one  obtains 

t(q',  q;  E)  =  E  n.,(coseqlq)  t^q',  q;  E)  ,  (6.26) 
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where 

i  00 

t^q',  q;  E)  =  ±   /  dr  f^q'r)  V(r)a>.(q,  k£;  r)  .  (6.27) 
For  a  symmetric  hard  core  potential  it  is  easy  to  show  that 

u)(q,  kE;  x)  =  (27r)-1/2  [e^x  -  e^xc  e-ikErc  eik£H  (6.28) 

and 

^'^  E)=4^q-[^  +  kE^^-kE)eiAqrC 

-  (q  -  kE)(q'  +  kE)  e_1'Aqrc]  ,  (6.29) 

where  Aq  =  (q1  -  q).    Although  equations  (6.28)  and  (6.29)  were  obtained 
from  the  equations  outlined  in  this  section,  one  would  arrive  at 
identical  results  by  considering  a  potential  of  the  form  V(x)=A9(rc-r) 
and  taking  the  limit  as  X  approaches  infinity.    For  more  complicated 
potentials,  such  as  an  exponential  potential  including  a  hard  core,  one 
must  generally  go  to  a  numerical  solution  for  w.  in  Eq.  (6.22).  The 
procedure  we  have  chosen  is  completely  analogous  to  the  comparison 
potential  method  outlined  in  Chapter  II  and  therefore  will  not  be 
repeated. 

3.    Results  for  Inelastic  Scattering 
The  quantity  of  interest  in  col  linear  scattering  is  the  transition 
probability  P.f  for  going  from  an  initial  state  i  to  a  final  state  f. 
In  the  case  of  inelastic  scattering  we  have  the  well  known  result 
[(Th68),  (Ek71)] 
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p.  =  ^T77T~  <pi  'n'  T11  Pin>  '  (6.30) 
n  n      p'P     1    1     1  ll1  1    1  v  ' 

1  1 


which  is  approximated  by 


2  2 

p      a  ^HJL_    iM(s)|2  (6  31) 


within  the  single  collision  approximation. 

We  will  now  consider  various  results  obtained  for  a  collinear 
problem  where  the  interaction  between  the  projectile  atom  and  its  colli- 
sion partner  is  a  hard  core.    These  results  are  analogous  to  those  found 
in  references  (Ek71),  (C170),  and  (Se66);  our  objective,  however,  is  to 
show  the  degree  of  discrepancy  of  the  full  single  collision  approxima- 
tion, and  the  peaking  approximation  within  the  framework  of  the  multiple 
collision  expansion  of  the  Faddeev  equations.    In  figures  (6-2a)  through 
(6-2c)  we  have  plotted  the  transition  probability  P10  versus  the  kinetic 


energy  E,  of  the  projectile  (Ei  =  for  the  case  m  =  0.125.  The 

1  1  2Mj 

exact  results  obtained  by  using  the  coupled  channels  method  of  Shuler 
and  Zwanzig  (Sh60)  are  given  in  Fig.  (6-2a).    The  single  collision  non- 
peaking  results,  or  those  obtained  from  Eq.  (6.10),  are  given  in  Fig. 
(6-2b);  and  the  peaking  results,  or  those  obtained  from  Eq.  (6.11)  are 
given  in  Fig.  (6-2c).    A  summary  of  the  various  integrals  and  analytical 
quantities  needed  in  these  calculations  is  given  in  Appendix  V.    We  note 
that  the  nonpeaking  results  are  closer  to  the  exact  results  than  the 
peaking  ones;  however,  the  discrepancy  is  not  large,  and  the  greater 
computational  simplicity  of  the  latter  approximation  makes  it  particu- 
larly attractive.    Figure  (6-2d)  gives  nonpeaking  results  that  have 
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been  renormalized  according  to  the  relation 

P  ,    =  P  ,    /[IP,],  (6.32) 
n  n      n'n  '  Ln,  n  n 

as  was  suggested  in  reference  (CI 70).    If  the  computation  of  the  tran- 
sition probabilities  were  exact,  then  the  expression  in  the  denominator 
of  Eq.  (6.32)  would  be  unity,  because  of  unitarity.    This  procedure 
brings  the  single  collision  results  into  closer  agreement  with  the 
exact  results;  however,  it  does  break  down  near  the  threshold  energies 
for  vibrational  excitation  because  of  the  presence  of  singularities. 
The  exact  reason  for  the  presence  of  these  singularities  in  the  single 
collision  approximation  is  not  clear.    Their  presence  has  been  noted  by 
other  authors  [ (CI 70) ,  (Ek71 ) ] .    By  having  considered  the  behaviour  of 

Mri  in  Eq*  (6*31)  near  tnresnold>  lt:  was  ascertained  that  the  singu- 
larity could  not  arise  there  but  from  the  1/Pj'  term  multiplying 
lM^j|2  in  (6.31).    This  would  indicate  that  the  single  collision 
approximation  to  M^,^  is  not  good  enough,  and  multiple  collision  terms 
are  needed.    On  physical  grounds  this  is  quite  reasonable,  since  near 
threshold  the  incoming  projectile  barely  has  enough  energy  to  excite 
the  diatomic;  and  thus,  the  likelihood  of  multiple  collisions  becomes 
quite  high.    The  importance  of  these  sequential  collisions  or  chattering 
in  collinear  scattering  has  been  verified  in  classical  trajectory  cal- 
culations [(Se69),  (Ke72)].    Analogous  results  to  those  discussed  above 
for  PjQ  are  given  for  P2Q  in  figures  (6-3a)  through  (6-3d).    We  see 
that  the  comparison  of  the  single  collision  results  with  the  exact 
results  is  more  favorable  in  this  case.    In  Table  (6-1)  we  list  exact, 
single  collision  and  peaking  results  for  P^q  versus       using  m  =  0.125 
and  m  =  0.5.    These  results  show  that  for  larger  values  of  m,  the  single 
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collision  and  peaking  approximations  break  down.    Secrest  (Se69)  has 
shown  that  as  m  becomes  larger  multiple  collisions  become  very  important; 
and  thus,  one  would  expect  the  single  collision  approximation  to  yield 
incorrect  results. 

Other  calculations  using  a  potential  of  the  type 


V3(x3) 


00   if    lxl  <  rc 


[  Ae_ot' x3 1  if    |x|  >  rc  (6.33) 


were  also  carried  out.    In  order  to  verify  the  reliability  of  the  com- 
parison potential  procedure  used  in  evaluating  the  t-matrix,  calcula- 
tions were  performed  using  large  values  of  A  and  a  in  Eq.  (6.33)  and 
a  comparison  was  made  with  the  corresponding  hard  core  results.  Table 
(6-2)  gives  one  such  study  for  P^g  using  the  peaking  approximation. 
The  parameters  used  in  this  set  of  calculations  were  m  =  0.5,  A  =  20.0 
and  rc  =  2.0,  where  the  units  are  those  compatible  with  the  scaled 
Jacobi  coordinates  given  by  equations  (6.1)  and  (6.2).    We  note  that 
the  exponential  potential  results  approach  the  hard  core  values  as  a 
increases  in  magnitude.    Another  important  test  of  the  single  collision 
approximation  is  its  applicability  to  softer  potentials.    In  the  litera- 
ture there  exists  an  extensive  study  of  collinear  scattering  using  a 
soft  exponential  potential  [(Sc75a),  (CI 70) ,  (An69),  (Se66)].  Unfor- 
tunately the  potential  used  is  not  that  given  in  (6.33),  but  one  that 
is  exponential  throughout  space,  i.e.  V3(x3)  =  fe~ax3  for  all  values  of 
x3.    Therefore,  it  is  not  at  all  clear  that  one  should  compare  the 
results  using  Eq.  (6.33)  with  those  reported  "in  the  literature.  The 
results  obtained  with  Eq.  (6.33)  using  small  a  dto  indicate  that  the 
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Table  (6-2)    Convergence  study  showing  how  the  transition  probability 
P^g  for  an  exponential  potential    approaches  that  for  a 
hard  core  potential  as  a  is  increased  in  magnitude 
(m  =  0.5,  A  =  20.0  and  rc  =  2.0  in  reduced  units). 


,  El 
(Units  of  hwj 

p10 

(Hard  Core) 

(a=4.0) 

(a=1.5) 

(a=0.5) 

1.1 

2.580 

2.745 

4.007 

16.661 

1.3 

1.458 

1.515 

2.046 

7.136 

1.5 

0.989 

1.015 

1.299 

3.661 

1.8 

0.586 

0.594 

0.703 

1.385 

2.0 

0.417 

0.420 

0.475 

0.739 

2.3 

0.250 

0.251 

0.265 

0.288 

2.5 

0.178 

0.177 

0.180 

0.154 
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single  collision  approximation  is  not  valid  in  the  same  range  of 
energies  that  were  studied  using  a  hard  core  potential,  i.e.  for 
projectile  energies  less  than  10  Tiu  units.    Typical  results  obtained 
are  given  in  Table  (6-3),  where  the  parameters  used  were  m  =  1/13, 
a  ■  0.1287,  rc  =  2.0,  A  =  10.0  and  E,  =  5.918  Tia)  units.    Note  that  for 
the  hard  core  results  shown  in  figures  (6.1)  and  (6.2)  the  transition 
probabilities  obtained  from  the  single  collision  approximation  were 
quite  reasonable  for  projectile  energies  one  quanta  beyond  threshold. 
This  is  not  the  case  for  the  transition  probabilities  quoted  in  Table 
(6.3),  unitarity  is  violated  even  at  several  quanta  above  threshold. 
Another  interesting  point  is  the  discrepancy  between  the  peaking  and 
nonpeaking  results  quoted  in  Table  (6-3),  indicating  that  the  peaking 
approximation  should  not  be  used  for  such  soft  potentials.  Similar 
cautions  on  the  validity  of  the  peaking  approximation  have  also  been 
noted  elsewhere  in  the  literature  (Co68).    The  reason  for  the  poor 
quality  of  the  single  collision  results  in  this  last  case  might  lie  in 
the  fact  that  for  very  soft  potentials,  the  distance  over  which  the 
potential  falls  to  zero  is  large  and  therefore    the  concept  of  a  local- 
ized collision  becomes  blurred  and  the  assumption  of  an  impulsive 
encounter  would  not  be  valid  [(Se66),  (Se69)]. 


4.    Dissociative  Collinear  Calculations 
In  the  previous  section  a  harmonic  oscillator  was  used  to  model 
the  behaviour  of  the  diatom  in  an  atom-diatom  collision.    When  this 
potential  is  replaced  by  a  more  realistic  one  such  as  the  Morse  poten- 
tial, one  can  allow  for  the  possibility  of  collision  induced  dissocia- 
tion (CID),  i.e.  those  processes  where 
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Table  (6-3)    Transition  probabilities  for  a  soft  exponential  potential, 
using  the  parameters  m  =  1/13,  a  =  0.1287,  rc  =  2.0, 
A  =  10  and  E-,  =  5.928  in  reduced  units. 


Single  Collision  Approximation 

Peaking  Approximation 

PlO 

0.893 

2.174 

P20 

1.299 

4.063 

P30 

1.317 

2.81 

p40 

0.878 

3.107 
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A  +  (B-C)  ->  A  +  B  +  C  . 
The  presence  of  three  unbound  particles  in  the  exit  channel  is  what 
makes  CID  very  difficult  to  study  from  a  theoretical  point  of  view. 
As  was  pointed  out  in  Chapter  I,  because  of  the  nature  of  the  exit 
channel,  coupled  channel  procedures  lead  to  a  continuous  infinite  set 
of  coupled  integrodifferential  equations  [(Wo75),  (Kn77)],  and  would 
thus  be  very  difficult  to  solve.    The  approach  we  shall  take,  suggested 
in  Section  V.3     with  the  single  collision  approximation,  offers  a 
much  more  tractable  procedure.    Specifically,  the  approximations  we 
shall  use  are 

(s) 

M01    =  Vp3f  q3i;  EP3f>  <chi!n>  <6-34) 

and 

M^'C'Pk)  ■  tP3k  f  dQ3  <$1c,|q1'xq1[n>  (6.35) 

the  one  dimensional  analogs  of  equations  (5.33)  and  (5.44).  Clearly 
the  amount  of  work  necessary  within  this  latter  approach  is  comparable 
in  magnitude  to  the  amount  that  was  spent  on  the  inelastic  problem,  and 
this  is  certainly  not  true  of  the  coupled  channel  procedure  [(Kn77), 
(Mi'75),  (Sh77),  (Wo75)]. 

Using  a  Morse  potential, 

V^)  =  D[l  -  e"axl]2,  (6.36) 

for  the  target  diatomic  does  not  affect  the  validity  of  the  transforma- 
tions used  in  Section  VI. I    to  obtain  Eq.  (6.6),  the  Schrodinger  equation 
for  a  particle  of  mass  m  hitting  an  oscillator  of  unit  mass  vibrating 
about  an  equilibrium  position  [(Se66),  (Sc75a),  ((Ek72)].    Energies  are 
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still  measured  in  units  of  1W,  but 


co  in  this  case  is 


oj  =  a(2D/m. ) 


1/2 


(6.37) 


implying  that  the  potential  parameters  a  and  D  are  related  by 


a  =  (2D) 


-1/2 


(6.38) 


The  Morse  vibrational  and  continuum  functions  <q^.  |n>  and  |q,'> 
present  in  equations  (6.34)  and  (6.35)  are  well  known  analytically 
[(Ek73),  (En69),  (F171),  (St35)]  and  there  is  no  problem  in  their 
evaluation.    The  same  is  true  for  the  form  factor 


Fj.j  for  inelastic  scattering  (Ek72).  A  summary  of  all  these  analytical 
expressions  may  be  found  in  Appendix  V. 

The  construction  of  the  transition  probabilities  from        must  be 
handled  with  care.    The  relation 


analogous  to  the  expression  for  P  ,  ,  Eq.  (6.30),  is  actually  a  differ- 
ential energy  transfer  probability.    In  Eq.  (6.40)  <e1|  characterizes 

the  continuum  state,         I,  of  the  diatomic,  and  Ae  =  e,'   -  e  .  For 

l  in 

inelastic  scattering  the  differential  energy  transfer  probability  would 
be  (Ek72) 


(6.39) 


one 


(6.40) 


Pn(Ae)  =  P 
n v  ' 


dAe  .-1 


(6.41) 


n'n 


drT1"  l 
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where  Ae  =  e  ■  -  e„  .    The  unitan'ty  relation  when  dissociation  is 
n  n 

allowed  would  be  (Le69) 

N  El 

Z  P.  +  /  Pn(Ae)d£e  =  1  ,  (6.42) 
n'=0    n  n  Dn 

where  the  integral  is  over  the  energy  range  allowed  by  energy  conserva- 
tion, which  in  this  case  is  from  the  threshold  of  dissociation 
(Dn  =  D  -  en)  for  the  Morse  oscillator  to  the  kinetic  energy,  E,,  avail- 
able from  the  projectile.    If  Eq.  (6.34)  is  used  to  obtain  Pn(Ae),  one 
also  has  the  problem  that  the  exit  channel  is  characterized  by  <Pj  p^  | 
where  <p^'  |  has  the  normalization 

<Pj  |Pj>  =  6(p1l  -  pj)    ,  (6.43) 

whereas  <e^'  |  is  normalized  on  the  energy  scale, 

<e1'  |e>  =  6(e1'  -  e)    .  (6.44) 

The  transformation  between  the  two  normalizations  is  given  by  the  rela- 
tion 

6(ei  -  e)  =  [6(p1I  -        +  S[p{  +  PJ)]p(P1)  ,  (6.45) 

where 

P(Pl)  =  H^-l"1  =  ^  (6.46) 
is  the  density  of  states  (Le69).    This  implies  that 

Pn(Ae)  =  I^T^        I<P1  PllT0llPin>  +  <Pi<-Pi)lT01[Pln>|2  (6'47> 
where  one  finds  that  the  particles  of  the  diatomic  are  being  scattered 
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to  the  left  and  right.    An  analogous  result  is  obtained  from  Eq.  (6.40), 
since  the  Morse  continuum  function  <X|]ei'  >  has  two  components  which 
asymptotically  correspond  to  plane  waves  moving  to  the  left  and  right. 
This  result  is  a  direct  consequence  of  the  boundary  conditions  and  the 

requirement  that  <x^|e^'>  be  a  bounded  function  (St35). 

(s)  (s  c  Pk) 

The  calculations  using  the  expressions  for  Mqj    and  '  have 

led  to  rather  unexpected  results.    When  the  differential  energy  transfer 
probability  from  the  ground  state  was  plotted  versus  Ae  for  the  system 
H  +  D2  (m  =  0.2,  n  =  0),  Mq^  led  to  a  discontinuity  at  the  threshold 
for  dissociation,  whereas  M^'C,PI<)  did  not.    This  may  be  seen  by  com- 
paring figures  (6-4)  and  (6-5).    The  interaction  potential       was  taken 
to  be  a  hard  core  and  the  value  12.28  in  units  of  hw  was  used  for  the 
Morse  parameter  D  in  Eq.  (6.38).    Eckelt  and  Korsch  (Ek73)  have  reported 
a  similar  study  on  this  system,  using  a  theory  which  leads  to  an  expres- 
sion analogous  to  that  for  M^,c'Pk^,  but  with  loss  of  time  reversal 
invariance.    On  the  other  hand  the  plane  wave  approximation  MqV > 
Eq.  (6.34),  has  been  applied  only  to  the  three-dimensional  scattering 
problem  (Sh77),  but  not  to  the  collinear  problem.    Because  of  the 
discrepancy  obtained  in  the  present  study  between  these  two  approxima- 
tions, the  difference  between  them  must  be  clearly  understood.    We  note 
that  it  is  important  to  investigate  that  region  of  the  energy  loss 
spectrum  encompassing  high  vibrational  excitation  and  transitions  to 
the  continuum.    This  was  not  done  in  the  three-dimensional  study  men- 
tioned above.    If  one  returns  to  the  formal  expressions  [equations  (5.1) 
and  (5.2)]  used  in  deriving  (6.34)  and  (6.35),  we  see  that  they  differ 
as  follows: 

M01)  =  <$o'T2  +  T3lV  (6-48) 
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and 

M^C)  =  <*0|W1+(T2  +  T3)|V 

=  «&0|T2  +  T3  +  T^Jo  +  T1GoT3|$0>  ,  (6.49) 

in  other  words  they  differ  by  the  presence  of  the  double  collision  terms 

(s,c)  (s) 
TlGoT2  ancl  Tl^oT3*    *n  t'ie  col  1  i near  case  Mqj      differs  from  Mqj  by 

the  term  TjGQT3.    Physically,  this  means  that  it  is  important  to  con- 
sider the  interaction  between  the  two  atoms  of  the  dissociating  diatomic. 

The  rather  drastic  difference  between  the  results  obtained  from  Mq^ 
(s  c) 

and  Mqj      in  the  collinear  case  may  on  the  other  hand  not  extend  to 
the  three-dimensional  results,  since  the  "dimensionality  bias"  may  play 
an  important  role  in  exaggerating  the  importance  of  the  double  collision 
terms  present  in  Eq.  (6.38).    Nevertheless,  this  question  must  be  addressed 
in  the  study  of  the  full  three-dimensional  CID  studies.    Another  aspect 
of  the  present  calculations  must  be  emphasized.    The  formalism  we  are 
using  is  based  on  the  Faddeev  equations    and  differs  from  that  used  by 
Eckel t,  Korsch  and  Philipp  (Ek74),  who  used  the  formalism  of  Bianchi 
and  Fabella  (B164).    In  this  latter  approach,  these  authors  have  intro- 
duced a  further  approximation  leading  to  what  is  known  as  the  impulse 
approximation  in  the  literature  [(Mc70),  (Ro67)].    Within  the  present 
approach,  this  approximation  would  amount  to  replacing  m1?'C'^  by 

M0l'C'n  =  q3i;  EqV  /  dQ3  <$5'  l<tf  ><qiln>>  (6.50) 

Pk     *      *  Q-?*2 
where  tg    (p3f,  q3i;  E  -  ^—  )  in  Eq.  (6.35)  has  been  replaced  by  the 


half-on-shell  t-matrix  element  t^k  (p3f,  q^;  E   *).    The  peaking  argu- 

*      *  *  3l 

ments  p3f,  q3i  and  Q3    are  determined  by  the  procedure  given  in  Section 
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VI.  1. .  Furthermore,  Eckel t  and  Korsch  (Ek71)  also  used  the  t-matrix 
for  the  potential 

°°   x  <  x 


V 


c 

0    x  >  xr  (6.51) 


instead  of  the  symmetric  hard  core  t-matrix,  Eq.  (6.29).    The  resulting 
t-matrix  from  Eq.  (6.51)  leads  to  an  unsymmetric  t-matrix.    In  order  to 
make  the  results  obey  the  principle  of  detailed  balance,  the  above 
authors  replaced  t3^  in  Eq.  (6.49)  by 

t3(+)  =  \  (t3(+)  -  t3<"))  .  (6.52) 

This  is  not  necessary  within  the  Faddeev  framework,  since  each  term  in 
the  multiple  collision  expansion  automatically  satisfies  time  reversi- 
bility (Mi 72b ) .    We  have  used  expression  (6.52)  as  it  applies  to  the 
unsymmetrical  hard  core  and  were  able  to  reproduce  all  the  results 
quoted  by  Eckelt  and  Korsch  (Ek73).    This  provided  a  check  for  the  com- 
puter programs  that  were  developed  to  do  the  present  study.    In  compar- 
ing the  calculations  obtained  from  M^,c'Pk^  and  those  from  Mq!|'c'*' 
using  equations  (6.51)  and  (6.52),  we  found  that  the  latter  yields 
results  that  were  nearly  unitary,  whereas  the  single-collision  results 
were  always  lower  in  magnitude  but  with  the  same  qualitative  behaviour. 
The  sole  reason  for  this  was  that  different  t-matrix  expressions  were 
used  in  the  two  calculations,    hence,  we  feel  that  this  indicates  that 
one  should  be  wary    in  attaching  too  much  significance  to  how  well 
results  obtained  with  approximations  such  as  (6.34),  (6.35)  and  (6.50) 
obey  the  unitarity  relation  Eq.  (6.42).    One  interesting  feature  of 
these  calculations  as  seen  by  comparing  figures  (6-5),  (6-6)  and  (6-7), 
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which  give  plots  of  Pq,  Pp  and  P,Q  versus  A?,  is  that  the  envelope 
always  reflects  the  initial  momentum  distribution  of  the  target  diatomic. 
This  supports  the  idea  discussed  in  Section    V.2    of  the  information 
obtainable  from  either  the  simple  overlap  volume  weighed  by  the  semi- 
classical  momentum  distribution,  or  from  the  form  factor  Fq^  present 
in  the  expression  for  Mi,    *     .    Figure  (6-6)  is  also  interesting  in 
that  it  is  an  example  of  how  the  present  approach  can  be  used  to  study 
scattering  of  excited  targets,  something  which  is  quite  difficult  to 
do  with  other  coupled  channel  procedures. 


CHAPTER  VII 

THREE-DIMENSIONAL  INELASTIC  SCATTERING  RESULTS 
In  this  chapter,  we  will  apply  the  single  collision  peaking 
approximation  to  inelastic  atom-diatom  scattering  processes  for  the 
case  of  full  three-dimensional  motion.    Specifically,  we  will  investi- 
gate the  systems  (Li  +  ,  CO)  and  (Li+,  N2),  and  discuss  various  numerical 
results  in  light  of  recent  experiments  [ (B676) ,  (Ea78)].    This  investi- 
gation will  serve  to  complement  other  numerical  results  found  in  the 
literature  [(Bo74),  (Ek74),  (Ph76)],  obtained  using  the  impulse  approx- 
imation [see  Section   V.4  ]. 

Within  the  present  approach,  we  will  approximate  the  inelastic 
scattering  amplitude  by  [see  Section    V.2  ] 

(s.Pk)  _     Pk    (2)        Pk  F(3)  ,7  n 

Mri    "  t2    Fi'i    t3    Fi'i  '  i7,1) 

where  the  target  form-factors  Fj?|  and  f|?|  are  defined  by  Eq.  (5.27). 
The  two-body  t-matrix  elements  t2Pk  and  t3Pk  are  off-the-energy  shell, 
and  their  arguments  are  chosen  by  a  procedure  analogous  to  that  used  in 
the  collinear  case  in  Section  VI. 6  .  One  of  the  goals  of  this  chapter 
will  be  to  illustrate  how  the  target  form-factors  play  a  major  role  in 
characterizing  the  energy  loss  spectra  in  atom-diatom  scattering 
processes. 

1.    Practical  Implementation  of  the  Peaking  Approximation 
The  differential  cross  section  for  inelastic  scattering  processes 
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is  formally  related  to  the  scattering  amplitude  by  the  relation  (Le69) 

P  1 

da      ,  1 A      ,2-rr, 4    2    1    ...  ,2 

m  I1  ^  =         Mi    Pf  \!Vl\  .  (7.2) 

Here,  we  continue  to  work  with  the  Jacobi  coordinates  introduced  in 
Chapter  V.    Since  only  inelastic  processes  will  be  considered,  we  have 
chosen  channel  1  as  the  entrance  and  exit  arrangement  channel.  Intro- 
ducing the  peaking  approximation  as  given  by  Eq.  (7.1),  we  see  that 

&(1'H>  ■  (f  )4  M?  £        F<?>        AP!)  +       f{?)  (-  ^  ^,|*  , 

(7.3) 

where 

•  Fl'l  {4c^l)  ~~  <nl  J"l  mj1le1mBlTAPl'rl|n1j1mj  >  (7.4) 

->■ 

and      is  the  position  operator  for  the  relative  coordinate  of  the 

diatom  (B-C).    Note  that  the  initial  and  final  states  of  the  target 

diatom  are  characterized  by  the  quantum  numbers  {n,,  j,,  m-  },  with  a 

l     1     J 1  ' 

prime  denoting  the  final  exit  state.    To  simplify  our  notation,  we  will 

drop  the  "1"  subscript  in  the  above  labels  and  the  "j"  subscript  in  the 

magnetic  quantum  number  m.. 

3 

In  the  experiments  we  shall  discuss  in  the  next  section,  the  mag- 
netic quantum  numbers  are  not  analyzed.    Therefore,  the  differential 
cross  section  expression  given  by  Eq.  (7.3)  wilH  be  replaced  by 

£  (-■J-.J.ejfc)  ■  $>*  i^flj+T)      lt2Pk  p«  ♦      r{?»l2,  (7.5, 

where  we  have  averaged  over  the  initial  magnetic  quantum  numbers  and 
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and  summed  over  the  final  magnetic  quantum  numbers. 
Introducing  the  vibrator  states 

<r1|njm>  =  uftJ (r^j  YJm(P1 )  (7.6) 

and  the  expansion  (Me61) 
•>  -> 

e^APl'rl  =  4tt    S  iX  Y*y  (Ap\)  Y^CP^j^lYAP^rj)  (7.7) 
Au 

into  Eq.  (7.4),  one  can  show  that  (Ph76) 

F<?>  =  4,    I  iA  C^.    f(?),   .y  Y*    (AP.)  .  (7.8) 
1  1  xy       J  m  jm    n'j'njA    Ay       1  v  ' 

where  we  have  defined  the  radial  form-factor 


and  the  coefficient 


°] VJ.  ■  '  d?l  1v  <?1>  V?l»  Vfl'  •  MO) 


The  integral  in  Eq.  (7.10)  is  well  known  (Ed74)  and  leads  to  the  relation 


c: 

j 'm' jm 


(-l)ra,(4,)-l/2[(2j.+l)(2x+l)(2j+l)jV2(^  JY°?  ^ 

\0  0  0/  \m'v  my 


(7.11) 

The  above  expression  for  FJ*]  [Eq.  (7.8)]  and  the  corresponding  equation 
(3) 

for  Fj, j,  lead  to  considerable  simplification  of  Eq.  (7.5).    The  alge- 
braic manipulations  are  quite  straightforward,  although  tedious.  Since 
they  are  given  in  the  literature  (Ph76),  we  will  simply  state  the  result 
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^  (n'j'-nj;  ef.f  )  =  (f  )4  ^      (2j+l)  E  {(2X+1) 

1  ^ 

JN  2  JI+J    pk    (3)       ,2}  (?12) 

\0  0  0/  n'j'njA  3  n'j'njA1 

Because  of  the  symmetry  properties  of  the  3j  coefficients  in  Eq.  (7.11), 
the  sum  in  Eq.  (7.12)  is  restricted  to  those  values  of  X  for  which  the 
relations 

|J'  -  j|  <  X  <  j  +  j  (7.13) 

and 

j'  +  j  +  X  =  even  (7.14) 
are  satisfied.    For  homonuclear  diatoms,  the  (-lp+J  'phase  factor  in 
Eq.  (7.12)  implies  that  odd  Aj  transitions  are  symmetry  forbidden 
[(Ek74),  (Bb76)]. 

The  differential  cross  section  can  now  be  obtained  by  the  evaluation 
of  the  radial  form  factors  defined  by  Eq.  (7.9)  and  of  the  off-shell 
two-body  t-matrix  elements.    In  the  present  investigation,  we  will 
assume  the  radial  functions  unj  to  be  harmonic  oscillator  wave  function 
with  modified  frequency  arguments  allowing  for  rotational  distortion 
[(Pa35),  (Yu76a)].    Adopting  the  treatment  given  by  Pauling  and  Wilson 
(Pa35),  the  resulting  expression  for  the  radial  wave  function  is 

where 

Snj(?)  =  [*1/2xnj2nn!]-1/2  e-[?/xn/  Hn(c/xnj),  (7.16) 
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xnj  =  [ti/(mlCu)nj)]1/2      ,  (7.18) 


r    2     6j  j+1  Beh  1/2 
co  •  -         +   ,  (7.19) 

1  e 


and 


j(j+l)Bph 

rnj  ■  re[l  +   e   j  t  {7  20) 

3j(j+l)Beh  +Icm1re2(u,e)2 

The  function  Hn  in  Eq.  (7.16)  is  the  ntn  Hermite  polynomial.  Anharmonic 
effects  could  be  incorporated  into  Eq.  (7.18)  by  replacing  ooe  by  an 
effective  vibrational  constant  (Yu76a) 

un  =  coe  -  u3exe(n  +  1/2)  +  Weye(n  +  1/2)2    .  (7.21) 

The  corresponding  vibrational -rotational  eigenenergies  are 

(nc)"1  Wnj  =  a)n(n  +  1/2)  +  Bej(j  +  1)  -  De[j(j  +  l)]2  ,  (7.22) 

where  the  equilibrium  rotational  constant  Be,  the  centrifugal  distortion 
De  and  un  are  ail  in  units  of  cm"1  (Le70). 

It  is  clear  that  the  presence  of  the  Gaussian  term  in  Eq.  (7.16) 
limits  the  range  of  integration  used  in  evaluating  the  radial  form- 
factors.    Therefore,  the  numerical  evaluation  of  Eq.  (7.9)  is  rather 
straightforward  and  no  complications  occur.    We  have  used  a  32-point 
Gauss  integration  procedure  as  suggested  by  Philflpp  and  coworkers  (Ph76). 
The  use  of  a  higher  order  quadrature  led  to  no  significant  improvement 
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for  the  present  case  in  which  only  lower  vibrational  levels  are  con- 
sidered. 

Finally,  one  should  note  that  the  evaluation  of  the  target  form- 
factor  or  the  differential  cross  section  requires  the  evaluation  of  3j 
symbols  to  very  high  j  values,  when  a  large  amount  of  rotational  excita- 
tion is  involved.    For  the  systems  we  shall  study  in  the  following 
section,  j  values  up  to  sixty  will  be  considered.    Often  the  straight- 
forward use  of  recursion  relations  found  in  standard  texts  [(Ed74), 
(Mi 61)]  will  lead  to  computational  problems.    To  avoid  this  difficulty, 
we  have  adopted  the  algorithms  given  by  Schulten  and  Gordon  (Sc75c) 
which  are  stable  to  very  high  j  values  (j>100). 

2.    Numerical  Results 
The  results  of  various  calculations  on  the  systems  (Li+,  CO)  and 
(Li+,  N2)  will  now  be  presented.    Within  the  present  many-body  approach, 
these  systems  are  convenient  because  the  mass  ratios  are  small  and 
repulsive  forces  are  dominant  in  the  collision  energies  of  interest. 
Using  equations  (5.33)  and  (5.34),  one  obtains  the  mass  ratios  0.15 
and  0.27  for  the  (Li+,  CO)  system,  and  the  mass  ratio  0.2  for  the 
(Li+,  N2)  system.    With  these  mass  ratios,  the  single  collision  approx- 
imation should  be  applicable,  since  the  ratios  are  less  than  0.5  [see 
Section  VI. 3  ].    The  potential  energy  hypersurfaces  have  been  recently 
calculated  by  Staemmler  [(St75),  (St76)].    Using  the  information  found 
in  these  references,  Micha  et  al .  (Mi78b)  have  extracted  two-body  poten- 
tial data  corresponding  to  the  collinear  configuration  of  each  of  the 
above  systems.    For  the  interaction  between  the  projectile  and  each  of 
the  target  atoms,  they  used  a  potential  of  the  form 
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V(r)  =  Ae_br    ,  (7.23) 

and  obtained  the  parameters  A  =  1908.51  eV  and  b  =  5.949  A"1  for  (Li+-C), 
A  =  5237.46  eV  and  b  =  4.8578  A"1  for  (Li+-0),  and  the  parameters 
A  =  3079.77  eV  and  b  =  5.2582  H_1  for  (Li+-N).    Using  such  a  simple 
potential  does  not  properly  describe  the  potential  energy  surface  at 
low  energies,  where  there  is  an  attractive  component.    However,  the 
experiments  we  shall  compare  to  are  at  high  energies  (4  to  8eV)  and  the 
scattering  angles  involved  are  beyond  the  rainbow  angle  (Bb'76).  Because 
the  above  parameters  lead  to  very  repulsive  potentials,  we  will  simply 
use  a  hard  core  interaction  potential  for  each  atom-atom  pair.  The 
data  we  have  used  for  each  of  the  two  systems  are  given  in  Table  (7-1). 

We  shall  restrict  the  present  investigation  to  the  study  of  rota- 
tional excitation  from  the  ground  vibrational  level.    This,  however,  is 
by  no  means  a  limitation  of  the  present  approach.    The  experiments  to 
which  we  shall  compare  our  calculations  are  time  of  flight  studies  on 
the  scattering  of  Li+  ions  from  ground  state  N2  and  CO  molecules  (Bo76) . 
In  particular,  we  shall  compare  our  calculated  differential  cross  sec- 
tions [in  arbitrary  units]  as  a  function  of  final  rotational  quantum 
number  jf  versus  those  obtained  from  experiment-    Typical  results  ob- 
tained for  the  (Li+,  CO)  system  are  shown  in  figures  (7-1),  (7-2)  and 
(7-3).    These  figures  correspond  to  a  center  of  mass  collision  energy 
of  4.23  eV  and  center  of  mass  scattering  angles  aofi  37.1°,  43.2°  and 
49.2°,  respectively.    Initially,  the  target  diatom  is  assumed  to  be 
in  its  ground  vibrational  and  rotational  level.   We  note  that  as  the 
scattering  angle  is  increased  the  overall  envelope  of  the  plotted  dis- 
tributions is  broadened.    The  same  effect  v/as  netted  in  the  experimental 


Table  (7-1)    Parameters  used  in  peaking  calculations 
[(He50),  (Mi  78b)] 


Spectroscopic 

Constants  CQ[Xh+]  N2[X1Eg] 

coe  (cm"1)  2170.21  2359.62 

Be  (cm'1)                     1.9313  2.010 

re  (A)                          1.1281  A  1.094  X 
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Atom- Atom  Pair  re(A) 

(Li+  -  C)  1.17 
(Li+  -  0)  1.05 
(Li+  -  N)  1.05 


O  «— i 

uo 

t— 1 

co+ 

+J 

•i—  CO 

_l  II 

c 

3 

f  a 

a;  _  i- 

S- 

+J  CD 

n3 

S- 

S-  "O 

+-> 

o  c: 

•r— 

4—  r0 

_Q 

S- 

(0 

CU  CU 

JD  CO 

E  CM 

c 

3  • 

o 

C  «* 

•r- 

II 

4-> 

E  t-< 

U 

3  LU 

Ol 

+-> 

CO 

C  « 

to  o 

10 

3  II 

to 

CX  4— 

o 

C 

S- 

CJ 

(O  - 

C  . — - 

o  o 

fO 

■1— 

4->  O 

4-> 

f0  — - 

+->  II 

O'— 

J-  -r- 

CD 

H- 

1—  •> 

4- 

CO  •!- 

c  c 

"C 

•p —  * — ' 

4-  >— ' 

4- 

O 

to  E 

3  a; 

+-> 

CO  +-> 

o 

S-  CO 

cu  >> 

>  CO 

V2 


o  o 

CM 

10  +  00 
■+->•»— 

-T—   1  II 

£=«--' 
3  fQ. 
CD  _t— I 

>)J=tQ- 
S_  -P  CD 

<o 

s-  t.  -a 
+->  o  c 

v—  t-  <0 

S_  S-  > 
(O  d)  d) 

- — xi  co 

E  CM 
C    =3  • 

o  c  «3- 
■I-  II 

4->   E  r-i 

a>  +-> 
«  n  « 
<c  o 
to  3  II 
in  a  <t- 
o  c 

S-  i— 

ore)" 
i—  o  o 

ro  *r —  »> 
•r-PO 
-t->  (D  

G  -4->  II 

CD  O  ' — - 

S-  S- 

CU  •<-) 
M-  r-  - 
"4-  fO 
•r-  C  C 

■a  -i  

If.  L_, 


3  cu 

to  +-> 

s-  i/i 

cu  >> 

>  to 


CM 

I 


cu 
s- 

CD 


174 


o  <— < 
o  o 

CVJ 


• 

-M 

•r- 

'i — 

_J  II 

C 

S   1 — 1 

rs 

s_ 

+->  O 

ro 

s_ 

1_  "O 

+j 

O  C 

4-  ro 

_Q 

S-  > 

t«3 

GJ  <U 

-Q  CO 

E  C\J 

c 

rs  • 

o 

c 

II 

-M 

£  rH 

3  UJ 

CI 

(/) 

C  •> 

ro  O 

i/> 

=3  II 

in 

CT  4- 

o 

C 

s_ 

u 

n3  •» 

c  ■ — 

O  O 

ro 

•r-  « 

•i— 

4-»  O 

■M 

fO  — - 

C 

+J  II 

OJ 

O  ■ — - 

s_ 

i-  •!- 

OJ 

4- 

1  ti 

4- 

(O  -r- 

■r- 

C  C 

-o 

•l—  — ' 

4-  ■— ' 

4- 

o 

</)  E 

Z3  0) 

•4-) 

to  +J 

o 

i-  V) 

<u  >, 

D_ 

>  </> 

o 

CD 


I-  00 


CD 
CO 


Si 


-  CM 


e 

<U3 


o 


o 


lO 
O 


o 
If) 
o 


CM 

O 


o 


177 


results  (B676).    In  each  of  these  calculations,  the  peak  in  the  distri- 
bution envelope  fell  at  a  j^-value  that  was  higher  than  the  one  observed 
experimentally.    In  Fig.  (7-1)  the  peak  occurs  at      -  29,  whereas  the 
peak  in  the  corresponding  experimental  curve  falls  at      =  20.  A 
similar  result  was  obtained  by  Micha  et  al .  (Mi 7Sb )  in  their  space-time 
correlation  function  approach.    In  their  studies,  they  also  assume  the 
single  collision  approximation.    As  a  possible  reason  for  this  discrep- 
ancy between  theory  and  experiment,  they  suggested  that  the  lock-in 
effect  of  neglected  long-range  anisotropy  would  lead  to  an  apparent 
increase  in  the  moment  of  inertia  of  the  target  diatomic.    Upon  treat- 
ing the  moment  of  inertia  as  a  parameter  and  increasing  its  value,  they 
were  able  to  obtain  agreement  with  experiment.    In  order  to  compare  the 
above  results  to  those  obtained  using  a  homonuclear  target,  we  have 
plotted  in  Fig.  (7-4)  the  differential  cross  section  versus  final  rota- 
tional quantum  number  of  the  (Li  ,  N2)  system,  using  the  same  energy 
and  scattering  angle  data  used  to  obtain  Fig.  (7-1).    We  note  that  the 
structure  obtained  for  the  (Li+,  N2)  system  is  much  simpler  than  that 
for  the  (Li+,  CO)  system.    Part  of  this  simplicity  is  due  to  the  fact 
that  N2  is  a  homonuclear  diatomic  and  odd  Aj  transitions  are  symmetry 
forbidden.    The  interference  pattern  in  the  (Li**  N2)  results  is  also 
less  pronounced  than  that  in  the  (Li  +  ,  CO)  resulltts. 

The  experimental  data  to  which  we  have  compared  the  above  results 
(BS76),  do     not  indicate  an  interference  structure.    This,  however,  is 
due  to  lack  of  resolution.    Other  experiments  om  (iLi  +  ,  CO)  have  recently 
been  done  which  show  various  maxima  in  the  energj  loss  spectra  for  rota- 
tional excitation  (Ea78).    In  Table  (7-2),  we  ccmpare  the  probabilities 
of  rotational  excitation  obtained  from  the  analysis  of  experimental 
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Table  (7-2)    Comparison  of  probabilities  of 
rotational  excitation  for  the 
(Li+,  CO)  system  (n^MO.O), 
nf=0,  EjM.28  eV  and  efi^p  =10° 
Column  I  gives  theoretical  results 
and  Column  II  the  ones  obtained 
from  experiment  (Ea78). 

P 

jf  I  II 


0 

0.018 

0.058 
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0.011 

0.100 

2 

0.002 

0.100 

3 

0.027 

0.040 
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0.045 

0.0 
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0.0 
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0.228 

0.0 
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0.0 
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0.0 

11 

0.102 

0.006 
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0.028 
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13 

0.015 
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0.004 
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results  (column  II)  to  those  obtained  theoretically  (column  I).  The 
center  of  mass  collision  energy  and  scattering  angle  are  4.28  eV  and 
10°  [8°  in  the  laboratory  frame].    The  rainbow  angle  is  at  5.5°  in  the 
laboratory  frame.    The  theoretical  probabilities  were  obtained  from 
the  relation 

As  noted  in  the  results  mentioned  earlier  in  this  section,  the  peak 
in  the  theoretical  probabilities  falls  at  a  higher  value  of  jx  than 
the  peak  found  in  the  experimental  results  [see  Table  (7-2)].  Because 
of  this  and  the  general  simplicity  of  the  present  model,  one  cannot  at 
this  time  state  whether  the  maxima  are  due  to  effects  coming  from  the 
potential  or  from  dynamical  effects  not  contained  in  Eq.  (7.12).  One 
should  also  note  that  the  scattering  angle  is  not  much  larger  than  the 
rainbow  angle,  and  therefore  the  experimental  data  quoted  in  Table  (7-2) 
are  not  the  most  favorable  to  our  model. 

The  results  we  have  been  discussing  have  all  been  obtained  using 
the  single  collision  peaking  approximation  as  it  arises  in  the  Faddeev 
formalism.    The  use  of  this  formalism  has  led  to  expressions  involving 
off-shell  t-matrix  elements.    On  the  other  hand,  the  impulse  approxima- 
tion [see  Section  VI. 4  ]  leads  to  half-on-shell  t-matrix  elements.  For 
the  present  case  of  a  hard  core  potential,  we  have  done  calculations 
using  off-shell,  half-on-shell  and  on-shell  t-matrix  elements  to  compute 
the  energy  loss  spectrum  given  in  Fig.  (7-1).    The  results  are  given  in 
Table  (7-3).    Column  I  corresponds  to  the  off-shell  results,  column  II 
to  the  half-on-shell  results  and  column  III  to  the  on-shell  results. 
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The  renormalized  results  were  obtained  by  scaling  the  entire  spectrum 
so  that  the  largest  peak  has  the  value  one.    We  note  that  the  differences 
in  all  cases  are  minor.    One  of  the  reasons  for  these  unexpected  small 
differences  is  that  the  magnitudes  of  the  momenta  that  occur  in  the 
present  calculations  are  all  similar.    We  caution,  however,  that  for 
softer  potentials  this  may  not  be  the  case.    As  was  pointed  out  in 
Chapter  III.  the  t-matrix  may  be  a  very  oscillatory  function  of  its 
momentum  arguments. 

It  was  argued  in  Section    V.2    that  the  momentum  overlap  volume 
give  by  Eq.  (5.36)  could  provide  a  rough  measure  of  the  probability 
amplitudes.    From  Eq.  (7.12),  it  is  clear  that  two  such  overlap  volumes 
are  involved.    Each  single  collision  term  leads  one  to  consider  two 
overlapping  spheres  as  defined  by  the  functions  <n,j'm'!r>  and  <r|njm> 
present  in  Eq.  (7.4).    Each  of  these  functions  is  normalized,  implying 
that  a  meaniful  comparison  to  the  single  collision  peaking  calculations 
can  be  obtained  only  by  normalizing  the  square  mot  of  the  volume  of 
the  sphere  determined  by  <n,j'm,jr>.    From  the  fLi+,  CO)  data  used  to 
generate  the  energy  loss  spectra  given  in  figures  (.7-1)  through  (7-3), 
we  obtained  the  rencrmalized  overlap  volumes  plotted  in  Fig.  (7-5) 
versus  jf.    From  the  three  curves  plotted  in  this  figure,  it  is  apparent 
that  the  peaks  are  shifting  to  higher  final  j-vallues  as  the  scattering 
angle  increases.    This  is  the  same  conclusion  that  was  drawn  from 
figures  (7-1)  through  (7-3).    There  is  also  a  gemeral  similarity  of 
the  shapes  of  these  curves  to  those  found  in  the  quanta 1  results. 
However,  the  momentum-overlap  model  predicts  a  greater  amount  of  rota- 
tional excitation.    The  interference  features  are  also  not  contained  in 
this  simpler  mode.    To  give  an  idea  of  the  behavior  of  the  form-factors 


c  •> 

O 

1-  o 

+J  - 

rC  O 
-t->  — - 
O  !l 
I-  — 


C  -r- 


3  • 

oo  E 
i-  0) 
<L)  4-> 

>  to 
>> 

00  oo 
0) 

E  1 — 1 
3  OO 

I—  <_)  r-H 
O 

>  -  r-~ 
+  ro 

Q--r-  || 
ro  —I  i-l 

i  f  Q. 

S_  -I-H 

>r<D 

O  4-> 

I  -o 

e  i-  c 

3  O  n3 
■»->  M- 
c  > 
tt)  S-  <u 
E  O  co 
O  -Q  CVI 
E  E  • 

<4-   C  II 
O  t-l 
E  LlJ 

+J  3 
O  +■>  « 

i—  c  o 

Q.  n3  II 

3  *t- 
<  CJC 


i 


0) 


CD 


185 


K~o  ty 


o 

CD 


3<  o 


ca  x 


O 
if) 


o  o 


9  O 


o 


O  X 


O 
CO 


O 


o  °  0 

N  co  O) 

n  ^  ^ 

O  O  X 


9 


O  X 


O  X 


O 


o 


in 

5 


o 
in 

a 


in 

a 


o  x- 


o 


o 


186 


in  the  quanta!  calculations,  we  have  plotted  the  radial  form-factor 
1 

f0jf OOjf  [A=Li+'  B=c  and  C=°J  as  a  Unction  of  jf  in  Fig.  (7-6).  This 
figure  shows  a  rather  pronounced  oscillatory  behavior  as  a  function  of 
jf,  but  with  a  rapidly  decaying  tail  for  high  jf-values.    This  decaying 
tail  leads  to  the  rapid  decrease  in  the  rotational  transition  probability 
for  high  j-values. 

An  attempt  was  made  to  extend  the  present  results  by  using  the 
exponential  potentials  given  in  the  introduction  to  this  section. 
However,  the  energies  involved  are  quite  high,  making  it  necessary  to 
use  more  than  one  hundred  partial  waves  in  the  evaluation  of  the  neces- 
sary t-matrix  elements.    Because  of  this,  the  partial  wave  t-matrix 
elements  involved  are  very  small  and  the  errors  accrued  at  the  high 
partial  waves  did  not  allow  convergence  in  the  full  t-matrix.    As  was 
pointed  out  in  Chapter  HI,  it  is  these  conditions  which  are  least  favor- 
able for  numerical  procedures  such  as  the  VPA  method. 

Generating  the  various  differential  cross  section  plots  shown  in 
this  section  required  a  rather  large  number  of  t-matrix  evaluations. 
It  is  the  evaluation  of  these  functions  which  can  be  quite  costly, 
when  numerical  methods  are  used.    Because  of  this,  plots  of  the  real 
and  imaginary  components  of  the  t£k  were  made  as  a  function  of  the  final 
rotational  quantum  number.    These  results  are  given  in  figures  (7-7a) 
and  (7-7b),  Re[t£k]  x  106  and  Im[t£k]  x  106  are  plotted  versus  jf. 

3.  Discussion 

At  hyperthermal  energies,  the  study  of  atom-diatom  scattering 
processes  using  coupled  channel  methods  becomes  prohibitively  expensive 
and  impractical.    It  was  shown  in  Chapter  IV  that  iterating  the  AGS 
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equations  led  to  a  multiple  collision  expansion,  which  in  effect  un- 
coupled the  various  transition  amplitudes.    Because  of  this,  approxima- 
tions based  on  the  multiple  collision  expansion  allow  for  the  direct 
computation  of  the  transition  probabilities  of  interest  without  the 
need  for  introducing  expansions  in  terms  of  intermediate  states.  For 
problems  such  as  collision  induced  dissociation,  we  saw  in  Chapter  VI 
that  enormous  simplifications  are  obtained  using  such  a  many-body 
approach.    We  conclude  then    that  approximations  such  as  the  ones  we 
have  pursued  in  this  dissertation  warrant  further  investigation  as  a 
viable  and  very  practical  approach  to  the  study  of  atom-diatom  colli- 
sion at  high  energies. 

From  a  practical  standpoint,  the  calculations  given  in  the  last 
section  for  inelastic  scattering  demonstrate  that  useful  qualitative 
information  is  obtainable  with  the  single  collision  peaking  approxima- 
tion.   However,  it  is  clear  that  if  more  reliable  quantitative  informa- 
tion is  desired,  one  must  input  more  realistic  two-body  t-matrix  data. 
In  light  of  the  difficulties  encountered  in  obtaining  converged  t-matrix 
elements  at  high  energies  and  large  momenta,  it  is  clear  that  an  improved 
algorithm  for  the  high  partial  waves  must  be  sought.    One  possibility  in 
this  regard  would  be  the  implementation  of  semiclassical  procedures. 
On  the  other  hand,  perhaps  the  most  practical  solution  would  be  to 
pursue  nonpartial  wave  techniques.    Even  though  the  Bateman  method  used 
in  Chapter  HI  did  not  work  for  singular  potentials,  one  could  modify  this 
procedure  using  a  two-potential  formula  (Ro67).    This  would  entail 
separating  the  hard  repulsive  core  from  the  tail  of  the  potential,  and 
using  the  Bateman  method  for  the  softer  component.    The  need  for  an 
accurate  numerical  procedure  still  remains,  however,  since  it  will  serve 
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as  a  comparison  to  more  approximate  procedures. 

The  collinear  studies  in  Chapter  VI  clearly  indicate  that  the  valid- 
ity of  the  peaking  approximation  may  be  limited  to  cases  where  the  po- 
tential is  very  repulsive  and  hard.    It  was  not  possible  to  pursue  this 
matter  in  the  case  of  three-dimensional  motion,  because  of  the  complex- 
ities arising  in  the  angular  momentum  analysis.    This  question  must  be 
answered,  however,  if  one  is  to  properly  judge  the  validity  of  the  single- 
collision  approximation.    We  pointed  out  in  Section    V.2    that  the  eval- 
uation of  the  single  collision  transition  amplitudes  would  entail  the 
introduction  of  partial  wave  expansions  for  the  various  quantities 
involved  [see  Eq.  (5.22)].    At  high  energies  this  would  not  be  practical, 
because  of  the  large  number  of  partial  waves  required  to  obtain  conver- 
gence.   In  order  to  circumvent  these  difficulties,  one  must  either  intro- 
duce decoupling  approximations,  or  make  use  of  multi-dimensional  quadra- 
ture techniques.    Given  the  simplicity  of  the  single  collision  term 
[see  Eq.  (5.24)],  the  latter  possibility  would  seem  thebest  alternative, 
at  least  in  the  case  of  repulsive  potentials.    If  the  potentials  support 
bound  states,  the  presence  of  poles  at  negative  energies  would  severely 
complicate  the  use  of  multi-dimensional  quadratures. 

Finally,  we  note  that  the  results  on  collision  induced  dissociation 
reported  in  Section  VI. 4    for  collinear  scattering  indicate  that  it  is 
essential  to  include  certain  double  collision  terms  if  reliable  results 
are  to  be  obtained.    Essentially,  it  was  demonstrated  that  the  interac- 
tion between  the  two  atoms  of  the  dissociating  dliiatomic  must  not  be 
neglected.    Recall  that  the  role  played  by  this  [potential  in  the  case 
of  inelastic  scattering  was  merely  one  of  deternrnning  the  initial  and 
final  momentum  distributions  of  the  target  diatom.    This  is  clearly  not 
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sufficient  when  dissociation  is  involved.    We  have  shown  in  Section  IV. 2 
that  adding  the  extra  double  collision  terms  [T^G0T2  and  T^GqT^  in  the 
expansion  for  Tq^]  is  equivalent  to  using  a  continuum  state  to  describe 
the  dissociating  target  diatom.    In  the  single  collision  approximation, 
a  plane  wave  was  used  for  this  purpose.    Given  the  disparity  between 
the  results  obtained  from  these  approximations,  we  feel  that  it  would  be 
very  important  to  extend  the  present  colli  near  study  to  the  case  of 
three-dimensional  motion. 


APPENDIX  I 
HARD  CORE  TWO- BODY  t-MATRIX 
We  will  now  derive  a  closed  analytical  expression  for  the  two- 
body  t-matrix  for  a  hard  core  potential  using  the  method  of  Laughlin 
and  Scott  (La68),  i.e.  the  t-matrix  elements  for  a  potential  of  the 
form 

V  =  X  6  (rc  -  r)  (Al.l) 

in  the  limit  of  X  going  to  infinity.    The  step  function  8  is  defined 
by  Eq.  (5.37).    The  off-shell  wave  function      can  be  obtained  immedi- 
ately from  the  appropriate  radial  differential  equation,  Eq.  (2.10), 
and  the  physical  boundary  conditions,  equations  (2.15)  and  (2.16).  The 
result  is 

ji(qO  ~m 

^(q.  V  r)  =  j^qr)  -  — -S   h  iJ(kErc).  (A1.2) 

.       hl  Vc> 

The  t-matrix  can  now  be  obtained  from  the  definition  of  t^q',  q;  E) 
given  by  equation  (2.13);  however,  as  was  pointed  out  in  Section  II. 2  , 
one  must  make  use  of  Eq.  (2.18).    From  this  we  obtain 

rc 

+  (q2  -  kE2)  /    j  (q'r)j ,(qr)dr].  (A1.3) 
o 
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The  integral  in  Eq.  (A1.3)  can  be  solved  analytically  using  relation 
(A2.2)  in  Appendix  II.    With  a  little  algebra  one  can  also  show  that 

d 

+   —     J  (q'r  )  5—  j,(qr  )  + 

n,2  n2      JTM  V  drc  JrM V 

-  ^(q'^c)  i}{<\rc)  ^-   In  hj1^^)       ,  (A1.4) 

where  the  j-j  and  hj1^  are  Bessel  functions  related  to  the  Ricatti- 
Bessel  functions  j-|  and  h}  '  through  the  relations  (Ab65) 

51  (z)  =  z  j^z)  (A1.5) 

and 

h(1)(z)  =  z  h(1)(z)    .  (A1.6) 

Equation  (A1.4)  has  been  reduced  to  the  form  given  by  van  Leeuwen  and 
Reiner  (vL61),  though  these  authors  used  a  derivation  different  from 
the  one  given  above.    It  should  be  pointed  out  that  there  exist  many 
errors  in  the  literature  (Bi73)  where  people  have  neglected  the  first 
term  in  equation  (A1.3).    The  reason  for  this  error  is  the  neglect  of 
Eq.  (2.18)  in  the  above  procedure  when  one  applies  Eq.  (2.13). 

It  is  interesting  to  note  that  in  obtaining  the  total  t-matrix  by 
summing  over  partial  waves,  the  first  and  second  term  in  Eq.  (A1.4)  lead 
to  the  simple  closed  expression 
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t(5-.5'E,.-(E-^i)r/itol  ♦ 


2m         C    o  2,„ 
2ir  Aq 


where 


Aq  =  |q'  -  qj  .  (A1.8) 


In  obtaining  Eq.  (A1.7),  we  have  made  use  of  Eq.  (2.12)  and  the  identi- 
ties (vL61) 


(2,)"3  /  drc  e-^'^c    d_  eiq-rc  „  ^  %  ^ 

c  2tt2  Aq 


In  order  to  clarify  the  meaning  of  the  first  term  in  expression  (A1.7), 
it  is  instructive  to  consider  the  expression  io  ,  Eq.  (A1.2),  sum  it  over 
partial  waves  to  obtain  w(q' ,  q;  E)  and  use  this  expression  to  obtain 
the  total  t-matrix.    Doing  this  one  obtains 


*    c\       ,Q2     c\     2  Jl(Aqrc) 
t(q',  q;  E)  =  E)  r  c   -L-  °_  + 

2m  c  9 

2tt  Aq 


r 


2 

/  drc  <  q1 1  rc>  ^_  E;  rc)  .  (Al.ll) 


2m  **  drc 


If  one  were  to  approximate  u>(q,  E;  rc)  by  a  plane  wave,  one  would  get 
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t(q',  q;  E)  -  -  (E  -  SiS  )r*  ^^lA    .  (M.12) 

2n  Aq 

Thus,  we  see  that  such  an  expression  is  a  plane-wave  or  Born  approxima- 
tion to  the  hard  core  t-matrix.    Other  approximations  based  on  addition 
formulas  such  as  Eq.  (ALIO)  could  be  pursued.    This  for  example  has 
been  done  by  Bethe,  Brandow  and  Petschek  (Be63)  in  applications  of  the 
Brueckner-Goldstone  theory  of  nuclear  matter;  however,  approximations 
analogous  to  those  used  by  these  authors  were  not  found  to  be  satis- 
factory in  the  present  case  of  atom-diatom  collisions  at  hyperthermia! 
energies. 


APPENDIX  II 

INTEGRALS  USED  IN  "COMPARISON  POTENTIAL"  METHOD 
In  Chapter  II  we  saw  the  need  for  integrals  of  the  product  of  two 
Ricatti-Bessel  functions,  and  in  particular,  for  the  comparison  poten- 
tial method,  integrals  such  as  those  found  in  Eq.  (2.57).    Consider  now 
the  integral  of  the  product  u-jV-j,  where      and  v,  satisfy  the  Ricatti- 
Bessel  differential  equation  (2.21),  i.e. 


drc  1  2  1 


0 


(A2.1a) 


and 


(A2.1b) 


It  is  well  known  (Me61)  that  the  desired  integral 

r=r, 


r  -  '2 

*      -  -      W[u,,  v,] 
/    dr  u  v    =       1  1 


11  <«?-,!> 


r=r< 


if  Ej  ^  E2;  on  the  otherhand  if  E   =  E?  (Fa71) 


(A2.2) 


'  UlVl 


dr 


(A2.3) 
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The  derivation  of  equations  (A2.2)  and  (A2.3)  is  straight  forward  and 
is  usually  given  in  standard  texts  on  quantum  mechanics  (Me61)  or  in 
the  theory  of  special  functions  [ (Fa71 ) ,  (6r66)].    If  we  now  replace 
one  of  the  Ricatti-Bessel  function,  say  v^ ,  in  the  left  hand  side  of 
Eq.  (A2.2)  by  a^1  ,  where  u.^1   satisfies  Eq.  (2.43)  in  the  interval 
[fj.  r^],  one  can  actually  follow  a  derivation  analogous  to  that  used 
to  obtain  expression  (A2.2),  and  the  result  is 


r.2  -  ,•  Wtw-i1',  iii] 
J    u.juy  dr  =   


r 


2       (q2-ak2)      ^  * 

+   L±—     /    J"  (qr)u  dr.  (A2.4) 


(qf-a.K.2)  rl 


In  order  to  clarify  this,  we  will  now  derive  Eq.  (A2.4).  The  starting 
point  is  the  identity 


1     dr2    1        1         r2  1 
-  V        u,1  +  (a^2  -  HI  +  1)  ^  .  (aEkE2  .  q2} ^  (qp)]J  =  0 

(A2.5) 

which  is  obviously  true  since  u,  and        satisfy  Eq.  (A2.1a)  and  Eq. 
(2.43)  respectively.    Rearranging  and  using  the  identity 

d  "2  2 

3FWK'U1]  -  «!1  ~2U1  "       72  ^  '  (A2.6) 

dr  1  dr  1 

one  obtains  the  desired  relation,  Eq.  (A2.4).    An  expression  analogous 
to  Eq.  (A2.3)  is  also  obtainable  for  the  case  in  which  q2  =  o^k?  in 
Eq.  (A2.4).    Here  one  would  use  the  identities 
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dr 


2  1 


u,  +  (q 


2  1(1  +  1) 
1  "  ~ 


V 


ru1[^?a)1i+  (q2-  HL+-U)  (apkF2-  q2)31(qr)]}  =  0  (A2.7) 

dr2    1         1  r2 


d  r  ,  d  4  . ,  d  "  ,  id", 
d?  [r(  drW  )(  dFul)  "  ^  dF  u^ 


K^ui)(^Lsi)  +  r(|:Gi)(^_w..).wi(  jL^,  (A2.8) 


and 


d    _ ,  ?    1(1  +  1)  .        -        „  o  -  .. 
dr  r2       )  ^i1  u-jl  =  2q^u1Lo11  + 

(ql2"  l0r2  ^  }  EnM  B7  ^  )  +  mj  f  lr  "])  "  ]  •  (A2-9) 

These  relations  can  easily  be  verified  using  equations  (A2.1a)  and  (2.43) 
The  final  answer  is 


r2  A 

+  (q2  -  aEkE2)  /    j1(qr)(  S-  ^  )dr}  .  (A2.10) 


Finally,  the  integral  on  the  right  hand  side  of  Eq.  (A2.10)  can  be 
obtained  from  the  relation 
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i  ™i<  i?s,)dr  -  ««i2- 
i 


r2  A  A 

2q2  /    U]v  dr)  (A2.ll) 


for      ^  q2.    Equation  (A2.ll)  may  be  proven  using  the  identity 


[r(  ^  V  "  V1  ^2  "1  +  (ql2  "        2  =  °'  (A2'12) 


which  follows  from  equations  (A2.1a)  and  (A2.1b). 


APPENDIX  III 
JWKB  STARTING  PROCEDURE 
In  propagating  the  numerical  solution  to  the  radial  Schrodinger 
equation 

2 

2-=-  *(r)  +  [k2  -  U(r)  -  }0  +  l)]^(r)  =  0,  (A3.1) 

one  generally  starts  at  some  rQ  within  the  first  turning  point  (Be60) 
such  that  ^  is  small.    Trial  calculations  are  then  performed  to  test 
the  sensitivity  of  the  final  wave  function  to  the  starting  point.  If 
the  potential  is  very  singular  near  the  origin,  this  sensitivity  is 
indeed  small  (Be60);  however,  if  the  potential  is  not  singular  or  if 
the  singularity  is  not  very  strong,  this  may  not  be  the  case.  The 
latter  observation  was  found  to  be  true  for  some  of  the  calculations 
reported  in  Chapter  III,  e.g.  the  off-shell  wave  function  studies  done 
for  the  Morse  and  Hulthe"n  potentials.    In  order  to  remove  this  arbi- 
trariness we  have  adopted  a  JWKB  procedure  proposed  by  Bell  urn  (Be75a). 

It  is  well  known  that  the  JWB  solution  to  Eq.  (A3.1)  near  the 
first  turning  point,  r.,  has  the  form  [(Mo58),  (La37)] 

♦<r)  »  Hfj-pj-HA  J1/3  (w(r))  +  B  J_1/3  (w(r))]    ,        (A3. 2) 

where  A  and  B  are  constants  specified  by  the  boundary  conditions, 

q(r)  =  [k2  -  U(r)  -  0  *  I/2)2  ^2  ,  (A3-3) 

r 
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r 

w(r)  =  /    q(r')  dr'       ,  (A3. 4) 

r. 


,th 

Eq.  (A3. 2)  one  finds  that  (Mo58) 


and  Jv  is  a  vtn  order  Bessel  function.    Using  the  asymptotic  form  of 


<Kr)  *  \  tq(r)]"1/2  e"|w(r)|  (A3.5) 


and 


-  i  rafrn"172  f1  d1n'^r)l  ,  d]w(r)l  -|w(r)| 

dr     "     2  lq^rjJ        {2  dr  dr        J  e  (A3. 6) 

for  r<r^..  Near  r  it  is  clear  that  one  can  use  a  linear  approximation 
to  the  effective  potential, 


Uef  ,  U(r)  +  Ifl-lil  .  k*  ♦         (r  -  rt)  .  (A3.7) 
r 

leading  to  the  relation 


Mil  „  I  iq(r)|1/2  e"{f|q(r)|3  (rt  "  r^Uef  "  k2)> 
dr     "  2  |qin|       6    6  (A3. 8) 

for  r  <  rt.    Equation  (A3. 8)  can  now  be  used  as  a  good  estimate  to  the 
magnitude  of  ^tidl)  to  the  left  of  the  first  turning  point.    In  the 
computations  reported  in  Chapter  III,  a  value  of  less  than  10"10 

was  found  to  be  suitable. 


APPENDIX  IV 


COORDINATES  FOR  THE  THREE-BODY  PROBLEM 


Consider  three  distinguishable  particles  (atoms)  A,  B  and  C.  We 

shall  associate  with  each  atom  a,  a  equal  to  A,  B  or  C,  the  atomic  mass 

itl,,  the  position  coordinate  r   and  the  momentum  coordinate  p  .    In  each 
or         r  a  ra 

arrangement  channel  i,  i  equal  to  1,  2  or  3,  of  the  three  particle 
system  one  is  concerned  with  a  free  atom  a  and  an  interacting  pair 
(By);  thus,  it  proves  convenient  to  work  with  relative  Jacobi  coordinates. 
These  relative  position  coordinates  shall  be  defined  according  to 
Fig.  (A4-1),  where  e.g.,  Rj  is  used  to  denote  the  position  of  particle 
A  relative  to  the  center  of  mass  of  the  (BC)  pair  with  the  relative 
position  coordinate  f-..    In  order  to  give  a  general  description  of  these 
coordinates,  it  is  useful  to  introduce  the  labels  (ijk)  and  (a8y),  which 
indicate  a  cyclic  permutation  of  the  channel  indices  (123)  and  the  par- 
ticle indices  (ABC),  respectively.    Furthermore,  we  shall  associate  the 
index  i  with  the  index  a,  j  with  3  and  k  with  y.    The  reason  for  this 
association  becomes  clear  if  one  considers  the  fact,  e.g.,  that  in 
channel  1  particle  A  is  free  and  particles  B  and  C  are  bound.  Ihe 
desired  Jacobi  coordinates  can  then  be  defined  as  follows, 


±      *  -*■ 


(A4.1) 


ri  =  ry  "  r8 


and 


mBrg  +  y 


(A4.2) 
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Figure  (A4-1)    Internal  and  relative  position 
vectors  r^  and  R-,  i  =  1,  2,  3, 
respectively  for  channels  1,  2 
and  3. 
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where  m,    =  m   +  m  .    Here  we  note  that  the  definition  for  r.  is  given 

f$Y      3      Y  i 

by  Eq.  (A4.1)  can  be  associated  with  the  mixed  index  set  (iyB)  which 

is  equal  to  (1CB),  (2AB)  or  (3BA).  One  can  further  show,  using  equations 

(A4.1)  and  (A4.2),  that  the  various  channel  relative  coordinates  satisfy 
the  relations 

-  -  53L  f .  -  (ij)  1-  (A4.3) 

1        maY     J  3 

and 


1  may  By     j  By 


where  (ij)  is  equal  to  +1  if  i  and  j  are  in  cyclic  order  or  -1  otherwise. 
M  is  the  total  mass  of  the  system.    In  each  case  the  transformation 
Jacobian  is  equal  to  unity. 

->-  -± 

At  this  point,  one  can  now  introduce  the  momenta  p..  and     ,  con- 
jugate to      and  R.  respectively.    These  relative  momenta  can  be  defined 

in  terms  of  the  p   as  follows, 

a 


and 


Pi  "S^Pv  -  SMa 
1     m3Y   Y       BY  B 


Qualitatively  then  P.  describes  the  relative  motion  of  atom  a  with 
respect  to  the  pair  (By)  having  the  relative  momentum  p^ ,  in  the  given 
channel  i.    We  can  further  introduce  the  channel  reduced  masses  m.  and 
M.  defined  by  the  relations, 

mRmY 

mi  ■  4r  (M-7) 
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and 


M . 


-  m 


(A4.8) 


to  cast  equations  (A4.5)  and  (A4.6)  into  the  form 


and 


5.  =  m,A  -  %  (A4.9) 

1  1    "ly  % 


j*.  =  „A     W]  .  (A4.10) 


The  relation  between  the  various  channel  momenta  can  now  be  written 
simply  as 

Pi  =  -  l&h  -  <U>  rn^rn—  Pj  (A4'U) 
1         m6y    J  mcxym6Y 


and 


APPENDIX  V 

ANALYTICAL  QUANTITIES  USED  IN  THE  COLLINEAR  SCATTERING  MODEL 
The  computation  of  the  peaking  transition  probabilities  for  vibra- 
tional excitation  and  collision  induced  dissociation  reported  in  Chapter 
VI  require  the  evaluation  of  the  form  factors: 

Fri  =  <n'  |e"iAx  |n>  (A5.1) 

and 

FQ1    ■  <e'|e-iAPx|n>    .  ,  (A5.2) 

For  the  case  of  inelastic  scattering,  F^,^  may  be  found  in  the  liter- 
ature for  both  the  harmonic  oscillator  and  the  Morse  potentials  [(Sh60), 
(Ra69),  (En69),  (Ek71),  (Ek72)].    In  the  case  of  a  harmonic  oscillator 
potential  [(Sh60),  (Ek71)]  one  can  easily  show  that 

k  -k 

mx  -    k  n, 


where 


Xkf0(-1}    k!(n--°k)!(n-k)!    -  (A5.3) 


n0  =  (AP)2/2  (A5.4) 


and  kmx  is  the  smaller  of  n  and  n'.  The  corresponding  relation  for  the 
Morse  oscillator  [(En69),  (Ek72)]  is  given  by 


210 


211 


(n-n'+l-if) 

F, =  A  ,  r  s  ,+  s    +  i  —  )  — —  rrr— !1 

11      n'n'n      n        a   '  (2s 


where 


n1  'n' 


x  3F2  ("n'  V  +  sn  +        •  n'  "  n  +  1T;  2sn  + 

AD 

-  n  +  i^-;  1)  ,  (A5.5) 


Vn=^  («)"lAP/a  ■  (A5-6) 
An     =  {a(2sn+l)n[n!r(2sn)]-1}1/2  ,  (A5.7) 


Sn     =  2D  -  n  -  1/2  ,  (A5.8) 

<«>„ « W  •  <«-9> 

3F2  is  a  hypergeometric  function  (Ab65)  defined  by  the  relation 

_  ,  n    (-n)  (a)  (b)  k 

3F2(-n,  a,  b;  c,  d;  z)  (c^(d^        fr    ,  (A5.10) 

and  the  quantities  a  and  D  are  the  Morse  potential  parameters  defined  by 
Eq.  (6.36). 

The  form  factor  FQ1  for  dissociative  scattering  in  the  case  of  a 
Morse  oscillator  can  be  derived  following  the  analysis  given  by  Strachan 
(St35).    Using  the  notation  given  by  Eckelt  and  Korsch  (Ek73),  one  obtains 


Fm  =  A  ,  {eT 

01     e  r(l+n-s+sn)r(l-s-sn-,f ) 


,i«(-s)  r(-s+sn+i^)r(i-2s)r(n+i-i^) 
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aP  AP  AP 

x  -F,  (-n,-s+s  s+sn+i^f;  l+2s  ,-n+i^:  1) 

5  c  n    a         rid  n  a 

+  [S  ->  -s]  ,  (A5.ll) 

where 

A  ,A 

A£,n  =  -el^(4D)"iAP/a,  (A5.12) 

A  ,  =  (8tt2c  '  )"1/4,  (A5.13) 

s  =  i(2e')1/2/a  (A5.14) 

6(s)  =  arg[r(-2s)/r(Js-s-D)]  =  -6(-s)  (A5.15) 

and  [s  ->  -s]  indicates  that  a  term  identical  to  the  first  one  must  be 
included  with  s  replaced  by  -s.    We  note  that  the  above  expression  for 
Fqj  requires  the  evaluation  of  complex  gamma  functions.    These  quantities 
may  be  evaluated  using  the  Stirling's  formula  (Fr74) 

lnr(z)  =  \  In (2*)  +  (z  -  h)  In  (z)  -  z 

+  rjr  -  +  — i-r  -  .  .  .  .  (A5.16) 

x"     360z  1260zD 

This  relation  is  valid  for  the  condition 

|arg(z)|  <  tt  (A5.17) 

and  will  yield  ten  digit  accuracy  if  |z|_>15.    The  recursion  relation 

r(z  +  1)  =  zr(z)  (A5.18) 

can  be  used  along  with  Eq.  (A5.16)  for  | z | <15 . 
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